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Heat Losses in Stranded Armature Conductors 
BY WALDO V. LYON 


Associate, A. I. E. E. 
Assistant Professor of Electrical Engineering, Massachusetts Institute of Technology. 


In the present paper, which is a continuation of one presented at the last annual convention, the author 5 
the method of complex hyperbolic functions to the solution of the problem of heat losses in stranded conductors embedded 
in rectangular slots. In the preceding paper the discussion was confined to solid conductors and to those having an 


infinite number of strands. 
thickness. 
appreciable thickness. 
4n the first paper. 


HE first step in the solution of the problem of 
heat losses in stranded armature conductors is 
to obtain an expression showing the relation 

between the currents in adjacent strands. If two 
strands are adjacent at any point, they are adjacent 
throughout their length. If the conductor is turned 
over in the end connection, the strand b, which is above 
a in one coil side, is below a in the next coil side, Fig. 1. 
The difference in pressure per half turn between these 
two strands is that between the lowest element of b 
and the highest element of a. If the conductor is 
turned over in the end connection, in the next half 
turn, this difference is between the lowest element of a 
and the highest element of b. The difference in pressure 
in each of these is cases given in equations (6d) and (6r). 
By the proper substitution, these equations can be 
written in terms of the currents in the adjacent strands 
and the current in the slot below them. See equations 
(Ta) and (7b). These two equations are similar in 
form. Between the points at which the strands are 
joined together, the sum of all of these half-turn differ- 
ences in pressure must be zero. This equation, (8), 
may be written in following form: 


(14094) = I ap) = 
ad ad 
(aa + 2 tanh m ) mo NT we 
ad 9 1 1p á 
m 7 
sinh ہت‎ 7 h 
m 


where < is caleulated for the embedded portion of 
the winding. Compare this with the fundamental 
equation (1) in the preceding paper which for the case 
of infinitesimal strands may be written in the follow- 
ing form: ۱ 


w dz 0 6 - e q + fx cdx) = 0; 
1+ سب‎ o 
1 


where a is also calculated for the embedded portion 
of the winding only, and l, and l: are respectively the 
length of the armature core and the length of an end 
turn. Notice that these equations are similar. The 


Presented by title only at the 10th Midwinter Convention of the 
A. I. E. E., New York, N. Y., February 15-17, 1922. 


İn the latter case, the insulation between the strands was assumed to have no appreciable 
İn the present paper, conductors are considered which have a finite number of strands separated by insulation of 
In the mathematical development which is to follow, free use is made of the results obtained 


first term in each is the vector difference between the 
currents in adjacent strands. The coefficients of 
the second terms are the same when the number of 
strands, m, Increases without limit and there is no 
insulating space, a, between them. It is further 
shown that the vector constants, To, are identical in 
the two equations. This similarity suggests that the 


BË 
d 


Fic. 1 


currents in the strands of the stranded conductor may 
be equal to those in corresponding divisions of a solid 
conductor of the same depth. In the latter case, 
the imaginary strands would of course not be sepa- 
rated. (Fig. 2.) 

Such proves to be the case provided the hyperbolic 
angular depth, 6 d, of the solid conductor has a certain 
value. This value is determined by the relation 
(Equation 10a): 


ili 


Fic. 2 
ad 0: a ad 
ine Bd _ 2m ) كك‎ + tanh pa) 
2m ad 
m lə 
+ 
sinh — h 


Were it not for this relation between the currents 
in the two cases the completion of the solution would 
be much more difficult. The reason is that the current 


199 


691121 


200 


up to any strand in the stranded conductor must be 
determined by a summation of vector quantities where- 
as in a solid conductor the current up to any point 
may be obtained by integration, a much simpler process. 
The current in any strand may also be obtained by 
integration between the proper limits. 

This equality of the currents in a stranded conductor 
and a solid one recalls the equality of currents and 
of voltages that exists in an artificial transmission line 
and a smooth line when there is a certain relation 
between their constants. 

We are now in a position to determine the expressions 
for the resistance drop in either the top or bottom 
strand of a conductor. It is only in the case of strands 


of infinitesimal thickness that the expressions are the - 


same. Compare equations (13b) and (14b) with equa- 
tion (7) in the preceding paper. They are identical 
except that M and N in the preceding paper are re- 
placed by similar functions, M’ and N’ and that there 
are terms involving 7” and S” which would disappear 
if the number of strands, m, were infinite and there 
were no insulating space, a, between them. We 
must also derive expressions for the pressure acting 
in conductors below a given one due to flux within 
the latter. These expressions are given in equations 
(15b) and (16b). These expressions are identical with 
the corresponding ones for solid or finely laminated 
conductors except that N’ replaces N, that T’ replaces 
a? d? to which it reduces if there are an infinite number 
of laminations with no space between them, and that 
there is an added term involving S’, a pure imaginary, 
which depends upon the insulating space between the 
laminations. 

Having determined these four pressure equations, 
the determination of the heating loss or of the leakage 
impedance follows precisely the same procedure that 
was used in the preceding paper. 

The author wishes to point out again a fact which 
he believes is not generally appreciated. A finely 
stranded conductor with the strands continuous 
throughout the whole coil will have a resistance ratio 
equal to that of a single finely stranded half turn of 
one half its depth, provided the coil has an even number 
of turns and the end connections are turned over on 
one side only. For conductors up to about one inch 
in depth, the alternating-current resistance is then but 
1 or 2 per cent greater than the direct-current resistance. 


Mathematical Analysis 


The description of windings considered follows. 
The embedded portions of the winding are in open 
reactangular slots; there are two coil sides per slot; 
the pitch is either full or fractional; each coil side has 
n layers or conductors, and each layer has m laminations 
or strands; the strands occupy the same relative posi- 
tions throughout the embedded portion of any half 
turn; the strands are arranged one above the other and 
have equal rectangular cross-sections; the strands are 
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separated by insulation of uniform thickness; when the 
conductors are turned over or twisted in the end con- 
nections it is done in such. à manner that any strand 
occupies the same relative position with respect to the 
bottom of the conductor of which it is a part that it 
does with respect to the top of the successive half 
turn; the coil as a whole is not turned over in the end 
connection; strands of any conductor may be joined 
at the beginning and end of a half turn, a full turn, or 
of a complete coil; the end connections may be turned 
over on neither side, on one side, or on both sides. 

The strands are designated by the subscripts q p. 
The q p strand is the p strand of the q conductor. The 
number of the strand may be counted from the bottom 
of the conductor of which it is a part, in which case the 
order of the strands is said to be direct, or from the top 
of the conductor of which it is a part, in which case the 
order of the strands is said to be reverse. Thus in a 
winding whose strands are continuous throughout a 
whole coil and whose end connnections are turned over 
on both sides the strands are in the direct order in one 
coil side and in the reverse order in the other coil side 
of the same coil. The number of the conductor is 
always counted from the bottom of the coil side of 
which it is a part. 

The armature currents are assumed to be balanced. 
In the end connections the r.m.s. current density 
throughout any strand is assumed to be constant. 
In the embedded portion of the coil the r. m. s. current 
density in any strand is variable. This vector r. m.s. 
current density in the embedded portion of the q? 
strand is: 


cosh a x 


C — —————— 
sinh 0 
m 


_m (7, ad 
Pod)” m 


— I» 24 tanh ae cosh a x 


1 
Ia, €f sinh at ) (1) 


c is the vector r. m.s. current density at points z 
centimeters from the bottom of the pq strand; d is 
the net depth? of the conductor; d/m is the depth of 
one strand; w is the width of the conductor or strand; 
Ia, is the current in the q p strand; Isa, is all of the 
current in the slot below the q p strand; 


_, 8 7*01 ۰ 


a? : 
ps 


1. See: Heat Lossesin the Conductors of Alternating- 
Current Machines presented at Annual Convention A. I. E. B., 
June 1921. (Equation (4). For solid conductors like the © 2 
strand Jo is all of the current in the slot below the conductor 
considered.) 

2. The net depth of the conductor is the depth it would have 
if the insulation between the strands were of zero thickness. 
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where f is the frequency, w, the width of the conductor, 
and s, that of the slot, and o is the resistivity of the 
material of the conductor at the working temperature. 
p is assumed the same for every strand. 

The current densities at the bottom, center and top 
of the q p strand are respectively Cops, Cape and ۰ 
Substitution in equation (1) shows that 
Since z = O, 


ad 
m m ad ad 
Ct d 06 "NET: - təyi and kə 
۱ m 
od 
; d m 2m 
Since 2 = Im? C = yd I+ ann ad (3) 
m 
Since 2 = ——, Capt = wq eo ad coth کیک‎ 
+ Ij di hz) (4) 
ad 
or O 1 7 
٩2۶ wd 7 sinh ad 
MEZ 
+ CT gg za tanh 2.6. (4a) 


It is shown in the paper already referred to that the 
flux, $4, within the embedded portion of any solid 
conductor, the وه‎ strand for example, due to current 
within the conductor, Ia, or example, and to all of 
the current in the slot below it, which in this case is 


Ioan, is 
1 | ۱ 
ay = jo "ri (15 + 25: } 


where /, is the length of the armature core. The flux, 
$a, between the q (p + 1) and the © strands is 


Arla 
5 


ad 
m 


tanh 2-2 (5) 


Ive (p+ 1), 


which may be written in a form similar to equation (5). 
a is the thickness of the insulation between the strands 
1 l p m 


əmirə” a (Ta, 5e) . aa 


m 
The proper combination of these five equations, viz., 
(2), (3), (4), (5) and (5a) will give the heat loss in any 
of the types of stranded conduetors we are considering. 
The leakage reactance due to flux within the embedded 
portion of the conductor may also be found. 


(Sa) 


RELATION BETWEEN THE CURRENTS IN ADJACENT 
STRANDS OF THE SAME CONDUCTOR 


For strands numbered in the direct order, the dif- 
ference in pressure between adjacent elements of ad- 
jacent strands is the difference in pressure between that 
in the topmost element of the pstrand and the bottom 
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element of the (p + 1) strand of the same conductor. 
For strands numbered in the reverse order this same 
difference is that between the pressures in the lowest 
element of the p strand and in the topmost element of 
the (p + 1) strand. The total difference in pressure 
between these adjacent elements must be zero between 
the points at which the strands are joined. 

For strands numbered in the direct order this dif- 
ference in pressure, Da, for a half turn is 


Dı = (lhipeeoe ون‎ +h Ino») 
-(hbtm ho ug Io 
+] Bər p cud aa) (6d) 
1 wd bq ip 1 m . 


l, and J, are respectively the lengths of the armature core 
and of the end connections for a half turn. The second 
term in each parenthesis is the resistance drop in the 
end connections. The third term in the second 
parenthesis is the drop in pressure in the q p strand 
due to the flux ¢,. The reactance drop in the end 
connections due to internal leakage flux is neglected, ?. e., 
the.current density is assumed to be uniform in each 
strand of the end connections. 


Substitute equations (2) and (4) in equation (6d) 
x ad 
m m 
D, = na? | liat + (ر‎ sinh ZÛ 
m 
— h و‎ (ç + » £2 (aa + tanh 2) 
{lex | -İan, $ coth شیک‎ 
d d 
+ Li I, کے‎ tanh — + 2 
but | Isap = Daptpo-—ls 
Thus: 
ad 
Di = —4 P (u جوت‎ th) سو‎ e 
sinh qq 
= hI ua 7ی‎ ۰ 
1 ۶ وا ون‎ + D m öm | 


For strands numbered in the reverse order the pressure 
difference, D,, for a half turn is: 


mv 
D, = ) :6م‎ + Ve Flap سح‎ icon 
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十 h p 一 一 


pon 
. aa 


- (lı p caps + سكم ما‎ Tap ) 


Substitute as before 
D, Bro} یں‎ 


۸0۳ D) 一 一 od Ed. tanh SŽ thls | 


—I I موه‎ = 


(6r) 


ad coth ad 
m m 


Ioap‏ 7 تک 


ed ) » + tanh 27-) thle | | 


Now İşa, = La و + ن +م)‎ ç + 9) 
Thus 
ad 
D, = — p 6 xu هجو )وج‎ 
m 


ad 


+ ly ed ( aa + 2 tanh 2) L 


In equation (6a) Da G+ و‎ = Do + 2 Io; 


where the first part is the current in the ilər below the 
q conductor and the second part is the current in the q 
conductor below its (p +1) strand. In equation 


(6b) Ise = Ison + Ia- 2 To where the first part 


is the current in the slot below the q conductor, the 
second part is the current in the q conductor and the 
third part, as before, is the current in the q conductor 
from the first strand to the p strand inclusive. 

With these substitutions, equations (6a) and (6b) 
become 


ed 
m 

Di = =F (^ .. 2. 

m 
- lı $2 (aa + 2tanh اف‎ 

- n S ( aa + 2tanh 27 ( 51s, (Ta) 

| ad 
_ MP x... 5 

D, = wd | (^ sinh & 4 23 (fa Ç + 1) Ia) 
m 

+h SE (aa + 2 tanh FÊ ( 7 EZ 


: HEAT LOSSES 


Transactions A. I. E .E. 


, ad 
یں ےچ‎ 2m 8) 
Between the points at which the strands are joined 
together the sum of all of these differences in pressure 


in adjacent half-turn elements, viz., 2 (Dz + D,), must 


(aa + 2 tanh $ ) 21 
1 


be zero. In general this sum may be written in the 
form 
ad 
. TT? lə 
ii — 
m 
- 24 ( وی‎ + 2 tanh Şİ ) fo 


ad ad Y? ۳ 
- £ (sa + 2tanh 27 ( Š To = 0 (8) 
where I, is the sum of the Isas and the ( — Isem 
—.1:)”s divided by 2 n, that is in general 
I, = . (Isai — İan — 14) 
We will now calculate the values of T, for the cases 
that we wish to consider. 
. Case'1. Strands joined at the beginning and end of 
each half turn. For the q conductor of the upper coil 
side 
I =(q— DI +] 0 
The first part of L, is the current in the upper eoil side 
below the q conductor and the second part is the current 
in the lower coil side which in general differs in phase 


by an angle 6. There are n conductors in each coil 


side. The current in each conductor is J amperes 
Iq = I). 

Case 2. Strands joined at the beginning and end of a 
whole turn but not turned over in the end connection. 
For the q conductor of the coil: 

Io =1⁄2((q— 1( I + 7 /0 + )- DI] 
= (q— DI + 0/2 I / 0 

Case 8. Strands joined at the beginning and end of a 
whole turn and turned over 0 the end connection. For 
the q conductor: 

Ty = 1/2{ )-1(1 +n1/0-(@-1)I-T} 
I,--I/24 
Case 4. Strands joined at the beginning and end of a 


whole coil of n turns. End connection not turned over on 
either side. 2 


I, = 37 2 ] 0-11۰1 /۵+ 6-1 ] 


Case 5. Strands joined at / 
whole coil of n turns. End connections turned over on 
one side only. The order of the strands in the first, 
or lowest, conductor of the upper coil side is direct, 
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that of the next half turn is reverse. The strands in 
the second conductor of the upper coil side are in the 
reverse order while those in the next half turn are in 
the direct order. 

In general 


I= əy 10/0-1)4(-s/68-1-141) 
-F(z/0--2-2—1) + )- 10/0 - 3 
| - 143 + ete. | I 
ly = JL | m/0- D +)- »/8- 1( 
° 2n — - 


—(0/9-1)4.(-n/6- 1) ete. | I 


(a) man even number 

I 0 = 一 I / 2 
(b) n an odd number 

I; = - 1/2 + 1/2 / 6 

Case 6. Strands joined at the beginning and end of a 

whole coil of n turns. End connections turned over on 
both sides. In the upper coil side the strands are all 
in the direct order while in the lower coil side they are 
all in the reverse order. 


1 


7( - 1 -و) -(4-1) + 0/9 2 لطم 
1 


+ 1/9 ساب 


Consider a solid conductor of the same width, və, 
and same net depth, d, but which has an angle 6 instead 
of a. The current density in this conductor is 

c = A’cosh Bx + B’sinh 6 x 

Imagine that this solid conductor is divided into m 
equal parts, z. e., strands, in the same manner that the 
actual conductor is divided. The difference between 


nI/ 0 
2 


Ba 

+ wd sinh سے‎ ۱ 

I, Sirm 8d (or ». e5)‏ —( + و"] 
2m‏ 


(See equation 2) 
But Cc (p+ رو‎ C'pe = (4 cosh ptl B d 
obanı? ga) 
m 
= ( 4’ cosh 2-12 B d 
m 
+ B' sinh FO eS 


Expand in terms of p/m 6 d and pd giving 


Ba) 
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CG + De— C'pe = (A” sinh p/m Bd 


; ., Bd 
+ B’ cosh p/m B d) 2sinh 5—- 


Thus: 
il 2 sink شك‎ 
Tyra, ga (A sinh p/m Bd 
Am ` 


+ B’ cosh p/m 8 d) 
The current in this conductor up to and including the 


part is 


Pa 
™ 


? 
“r, -{ wc'dz =w/B(A’sinh p/m 0 
1 


+ B’ cosh p/m d) — w/ 8 B’ 


We will now show that these values of current can 
be made to satisfy the relation that has been established 
between the currents in adjacent strands of the actual 
conductor. Substitute these values of (I',+ 1 — I',) 


? f 
and 2 1,” in equation (8). That is assume that the 
1 


currents in the imaginary strands of the solid conductor 
are exactly equal to the currents in the corresponding 
strands of the actual conductor 


ad Ba 
m وا‎ Yod 2sinh? 2m - 
sinh € d h m Bd 
m 2m 
( A” sinh p/m Bd + B’ cosh p/m 8 رہ‎ - 2° 
ad ad ad 
(aa + 2 tanh Fo) كك ىر‎ ) =a + 2 tanh 27.) 


w/B ( A’ sinh p/m Bd + B! cosh p/m B d — B') 


= 0 (8a) 
This equation (8a) is satisfied if: 


ad 2sinh? d 


m lə 
sinh d 0 h m 
m m 
d ad 


0 
- SE (=a + 2tanh 27. ( w/8 = 0 


wd 


(9a) 


and: 


ad | ad 
- “Gü (aa + 2tanh 27 ) I. 


a d a d P 
S ( aa + 2tanh $7) w/BB = 0 (9b) 


These conditions readily reduce to 
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ad ad. ad. 
7 gö ) St + tanh FÊ) 

sinh = ad | (10a) 

m 2 

ad 7 h 

sinh 7 

and = Bw (10b) 


If the total current in the conductor, actual or 
equivalent, is I, the vector constant A’ is determined 
as in the preceding paper. 


po ۱ I 
ac v( sinh 8 d 
By equation (3) the current in the p strand is 


sinh Bd 
I wd 2m 
"7m Ba 


— I, tanh Fe) (10c) 


( A’ cosh لتحم‎ Bd + B’ sinh P= 0 二 1/2 Bd) 


8d cosh 二 二 ^. 2 Bd 
o sinh an sinh ^ d 
- I, tanh E$ cosh 2 — 1⁄2 وم‎ 
+ I, sinh نت‎ 2 Bd (11) 
The current in the conduetor below the م‎ strand is 
۵-1 c 
> I, = f wc'dz 
1 
0 
= u/ B ( A' sinh حل‎ L gd قرو رب‎ - Bd- Bi) 
sinh 2—1 وم‎ 
| . snhfd 
_ Bd .. p—1 
I,tanh — sinh p d 
+ 1, cosh Z : Bd—I, (12) 


The total current below the q p strand of the actual 


conductor is 
p-l 
I bq + > I 2p; 
1 
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where Isa is the current in the slot below the q con- 
ductor. 

If the current in any conductor and all of that below 
it in the slot are given the copper loss in the conductor 
may be calculated as described in the preceding paper. 
We are thus able to calculate the copper loss in any 
strand o° any conductor. This method of calculation 
is far too laborious and we shall content ourselves with 
ealeulating the loss for a half turn (Case 1), a single 
turn (Cases 2 and 3) or for a single coil (Cases 4, 5 and 
6). To do this it is necessary to obtain expressions 
for the voltage drop per half turn in the top element? 
of the top strand, in the bottom element? of the bottom 
strand and the voltage drop due to flux within a single 
conductor in all conductors below it. 

The voltage drop per half turn in the topmost- ele- 


` ment of the q p strand due to resistance and leakage 


below the q (p + 1) element is, 


lı pcan + م ما‎ lə + h p — a Tae + کے و‎ aa 

The third term in this expression is the voltage drop 
due to flux within the insulation—of thickness a—that 
is immediately above the q 2 strand. It is necessary 
to include this term in order that the voltage expressions 
about to be derived will be similar in certain respects 
to the general equation (8) for the currents in the 
strands. 


By equation (4a) this reduces to 
m EN 
wd ) ۰ kun 
m 


(aa + tanh $ ) 


He outs ae om 


. The voltage drop per half turn in the topmost element 

of the top strand of the q conductor of the upper coil 
side due to resistance and leakage flux below this ele- 
ment is, 


ad 
_ m lə 
drop = م‎ T ly (mr 7 
Sinn | 


1 (13) 


This pressure is used when the strands are numbered 
in the direct order. The loss in pressure per half 
turn in the lovest element of the bottom strand of the 
q p conductor due to resistance and leakage flux belov 
this element is: 


lı م‎ Cops + كك م وا‎ Tap 


^. 3. Due to resistance and flux below it. 
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۰ ۰ : 5 ۱ m~l 
By equation (2) this — | 00 z Io 
drop = p — | —m رر پر - رر ہے‎ + 2 21 ) 
۲ 0 d 1 sinh ad 7 qp = M Loa qalir qp 
m From equation (8) we have, 
ad ad 
— رون‎ 一 一 m Un. | zt (aa + 2tanh 27 (2 م2‎ 

The loss in pressure per half turn in the lowest dd 
element of the bottom—7. e. the qı—strand of the q le 
conductor of the upper coil side due to resistance and = ( MGE + U ) (Lap — La g-o 
leakage flux below this element is, sinh “nr 

ad 
一 一 ۱ ad ad 
drop = p— h (Fart : (۶ - S (sa + 2tanh FE ) I 
E m If the (m — 1) equations for the adjacent strands are 
added together we have, 
d d d d m-19-1 
- ba (a+ tanh 27-) + Ion “ə . aa که‎ (aa + 2tanh 27 ) 2 2 lə 
(4) y ہے‎ 

This pressure is used when the strands are numbered = ( 605” ) (Iam — Ia) 
in the reverse order. The equation is written in this sinh Ron 


form so that it will be similar to equation (13). When the 

strands are continuous from one half turn to the next 

and the end connection is turned over between suc- 

cessive half turns the top strand of one half turn be- 

comes the bottom strand of the next half turn. Thus 

in equations (13) and (14) the Iam is the same as I q. 
The flux in the q conductor is, 


ga = > (dap + Qa) 


1 
By equations (5) and (5a) this is, 


ad 
2m 


1 1 m ad 
J o P wd m 


| tanh > ) ap +2 "T 


+ هيه‎ (To + Do) | 


— 
سح 


1 
— p 


i ad 
Po "dl L (aa + tanh Fr) 


ad ) 
2m 

The summation in the second term in this expression 
for the flux may be written, 


2 Isa = 06 + 0 


+ رون‎ + Ia 
+ رون‎ + Ia + Ia 
+ loa ا‎ 


Dogs ad ( aa + 2 tanh‏ ا 
m‏ 1 


tha + و‎ 


ni S3 ) aa + 2tanh 22) 7, 


Therefore the flux in the q conductor when the strands 
are numbered in the direct order is, 


ad ad 
— l H—. 


(aa + 2tanh ge ) m La 


ad 


m lə 
+) پگ‎ + 
sinh € d h 
m 


) Ge. — Io) 


— (m — 1) ® (aa + 2tanh ==) n | (15) 


If the strands are numbered in the reverse order, 
so that the first strand is at the top andthe mth strand 
is at the bottom, the current in the conductor below 


t m 
the q م‎ strand is J — 2 I+, and the 2 Isep is now: 
۱ 1 


2 Iras = La + 1 - Ia 


d a a) 
+ Isom +I- 2 Ley 
کی5 سج“‎ Q m 
“mə. 

ea ھی‎ zi > İn 


7-1 
= mIsem + (m= 1( 1-8 > Te 
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Therefore the flux in the q conductor v/hen the strands making the substitutions given in equations (17) and 
are numbered in the reverse order is (18), and remembering that the Iam in (17) is the same 


1 m ad öd as Ia in (14). We will also make the substitution 
وه‎ = jo Pod i| 45 (e a + tanh öm ) shown in equation (10a). 
Equation (13) becomes: ۱ 
428 aa + 2 tanh نے‎ ( ) Las + (m — 1) D ۲ ( + tanh ©) 
od drop = fy h| tad — ed mor پ ہی‎ 
m le tanh 2m 
-( + 7 ) 6 - La) 
sinh .* 9 ۲ 
di (“5 + tanh $) Bd 
d d | + Io & d ———54—— tanh 一 一 
+ (m 一 1) S (aa 2 FÊ) İas : tanh i ۱ 2 
The loss in pressure as derived in equations (13) and ۱ 0 
(14) and the flux within a conductor depend upon the + (Ly — Io) a a (° nh سس‎ ) 
current in the top strand and the difference between the 
currents in the top and bottom strands. The current 28 
in the p strand is determined by equation (11. + (I + (رآ‎ od. — | (13a) 
The current in the top strand numbered in direct ۱ 
order, is I am. . Equation (14) becomes: 
aa ad 
۱ ( + tanh $5) 
Is, = 2 sinh -$ Í een pe drop = يشب‎ h| Taa 2 pg coth pa 
۱ 2m sinh 6 d tanh — 
2m 
— I, tanh Ñ cosh سس‎ öd aa a 
. 5 + tanh 2m d 
+ I o Q d 8 0 tanh 
: 1/2 tanh om 
+ Tosinh E ga ) 
aa a 
۱ M (Is + L + D ad( 1 ho} 
Expand in terms of 6 d and 2 | | 
Bd | Bd +I; ad. a (14a) 
Tam = 2 sinh’ ç Hİ ) ۱ ۱ 
ic Equation (15) becomes: 
d Be d 
+ I, (tanh كك‎ 5 coth ہے‎ — 1 ) | (17) _ 1 
be = سحب‎ aq il ü + 270 
Likewise: ۲ 
fa ۱۳ tanh r) 
th — 2m 8 d 
əz HI CON 2m tanh — 
——. m B d tanh fe 2 
Bd B d 
- I, ) tanh كوك‎ coth $ — 1) —— ad) 
Thus: 
a 
(Lem — Iq) = 2sinh? Be | (I + 0 2 | (15a) 
۲ , Equation (16) becomes: 
( tanh fe coth 2.4 — 1 ) | (18) 1 5 
Pa 20 2| - (I + 2 L) 
The pressure equations which determine the resist- 
ance and reactance drops, viz. equations (13) to (16) aa + tanh ad ) 
inclusive are in terms of these currents in the top and ad 2 2m / tanh Bd 
bottom strands in addition to the conductor current tanh Ba 2 
and the current 7. We will rewrite these equations, 2m 


LYON: 
ad TE 
2m 


+Iad. M (16a) 
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In each of these four equations I, is the current in the 
slot below the conductor in question. 

The writing of these equations is much simplified 
if we let: 


a ad 
1 ( > t tanh zm ) 
d coth Bd 
2m 


1 (< + tanh 3) Bd 
: 5 tanh 5— 
m 
= ۸۷۲ < ۸۷,۲ +N, 
7 aa ad 
LT. ad 2 m ( : + tanh 27) 
= T' = ۲, + 1, 
l — Q _ 8 72 wf 
EL adaa =S' = j 2: d 


and IT (lı + L) = R (the d-c. resistance of a half 
tun) ` 

These equations may now be rewritten in simplified 
form. The resistance drop per half turn in the top 
element of the top strand of a conductor (13a) becomes: 


DN TO 一 


drop = R[ IM’ + I) سس‎ 


+ (I + 70 5/2 | (13b) 

If the strands are turned over in the end connections 
the resistance drop per half turn in the bottom element 
of the foregoing strand, which is now the lowest in 
the conductor (14a) becomes: 


I, N’ 


T 


drop - ام‎ IM’ + 


+1, 5/2 | (14b) 


The average resistance drop per half turn for the 
six cases considered is obtained from these two equa- 
tions. By definition I, has such a value that the term 
involving 7” in this average resistance drop is zero in 
every case. In order to simplify the expression we 
will let 7, represent the coeffcient of S”/2. The 
average resistance drop per half turn is: 

/ , 
1 N 7 I T ) 


av. drop =F ) IM’ + ^55 (19) 
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I, is the average sum of the (I + I,)’s in equation 
(13a) and the I; s in equation (14a) taken in the 
proper combination. 

Case 1. I. =I + و)‎ - 1 + 2/60( 7 upper coil side 
= (q+ m / 0) I upper coil side 
= q I lower coil side 

Case 2. I, — 1/2 | (I + (q^ 1+ m / 0)1] 

+ ]] + )- 1( 7[ | 
= (q + ۰/2 13 
Case 3. I, = 1/2 ) ] + (q— 1+ 0/ 6 ( 7 [ 
|. *Í[(a- DII] 
= (q — 1/2 + n/2/ 80) I 


Case 4. I, = 1/n2 (q + n/2/0)I 
-(5— + ۰/2 / 9 ( 1 


Case.5. I, 


)0( + (0/6 +4 ]و 
)1+ 1( + (60 07 +4 1) + 
ete. İT‏ + )2( + )424479 1) + 


The first term within the parenthesis is the current 
in the conductor plus the current below the first con- 
ductor of the upper coil side (13b). The second term 
is the current belov the first conductor of the lover 
coil side (14b). This is zero. The third term is the 
current below the second conductor of the upper coil 
side (14b). The fourth term is the current in the 
eonductor plus the current below the second conductor 
of the lower coil side (13b). Notice that this expression 
for I, may be written: 


1, - əy. 04/04 6+ ۰/۵ 
+ )6 + 0/0( ٩ 
+ (2n-1+n/0))1 


1 n 
I= əy 2@q-1+mn/0)1 
I, = 10/2 + n/2/ 0) I 


Case 6. I, = 1/nZ (q — 1/2 + n/2/ 0)1 


= (n/2 + n/2/ 0) I 

The pressure acting in all conductors below the q 
conductor due to flux within the latter is 7 o öç. 
Written in the simplified form this (15a) becomes: 
drop = R ] )1/2 + Io) N’ + (L,— 19) T’ + I/2 S! 

(15b) 

If the strands are turned over in the end connections, 
the pressure due to flux within the next succeeding 
half turn (16a) becomes: 
drop = R ] - (1/2 + L) N' + (L, + L, + DT’ 


+ 1/2 S' (16b) 
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Reference to the preceding paper shows that these 
pressure expressions are similar to those already derived 
for solid and finely laminated conduetors except for 
the added terms involving S”. M” replaces M, N” 
replaces N and T’ replaces o? d?. The expressions 
become identica! in the limiting case of an infinite 
number of strands with no insulation between them. 


CALCULATION OF COPFER LOSS 


The method of calculation is the same as that used 
in the preceding paper. 


Case 1. When the strands are joined at the begin- 
ning and end of a half turn the copper ‘oss in the half 
turn is symbolically: 


loss = 1. R(IM' + 1/2 N: + 1,28") 
+ Is. R[(I/2 + L,/2) N’ + I/2 8’) 
Im this ease notice that I, = L, 


The first term is the power loss due to the current in 
the conductor and the resistance drop in the topmost, 
element of the topmost strand. The second term is 
the power due to the current in the slot below this 
conductor and to the pressure acting on it produced by 
the flux within the conductor. Expanding this expres- 
sion gives: 


power =I. RIM’ +I. REN +I. RSS 
+1. 8 43 +h. RDN +L. RSS’. 


The first term is R P M,'. The sum of the second 
and fourth terms is R II,N,' cos 0,. The fifth 
term is R I N,’ and since S” is a pure imaginary the 
sum of the third and sixth terms is zero. The expres- 
sion for the power may be written: 


loss = 8 ]]۶ M,' + (I. L, cos 0, + Ip) N,'] 


The phase angle 6, is between the current, J, in 
the conductor and the total current, I», in the slot 
below it. Here the letters J and J, represent the 
numerical values of the currents. 

If this copper loss due to alternating current is 
divided by the loss due to the same amount of direct 
current, we obtain the ratio of alternating to direct- 
current resistance. 

(a) The ratio of alternating to direct-current resist- 
ance for a single half turn is: 


K = IM,” + [ (1:/D? + I/I COS 7 N,” 
(b) The average resistance ratio for a one-coil-side- 
per-slot winding having n layers is: 


K = 1/05 (M; +[(q— 1: + (a — (۱ 


= ) 1, + RN.) 


(e) The average resistance ratio for the upper coil 
side of a, two-coil-side-per-slot winding having n 
layers per coil side reduces to 
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0 is the phase angle between the currents in the 
upper and lower coil sides. 

If the strands are joined at the beginning and end 
of a whole turn the loss ratio is a little more difficult 
to caleulate inasmuch as the losses in the two half 
turns are different. There are two cases to consider, 
one in which the end connections are not turned over 
and one in which they are turned over. 

Case 2. Strands joined at the beginning and end 
of a single turn; end connections not turned over. 
The heating loss in the whole turn may. be expressed 
symbolically as: 
loss = 2I. ] )۲ ۸۲۲ + ۲0/2 ۷۲ + I, S'/2) 


“4 


+I, . 2] )1/2 + Io) N’ + L, — Io) T’ 
+ I/2 S!” 

+I”. R[(I/2 + I) N’ (15 — Io, T’ 
۱ + 7۲/2 S" 


The first, term is the power due to the current in the 
conductor and the resistance drop in the top element 
of the top strand of the turn. The second term is the 
power due to the current, T:”, below the upper half 
turn and the pressure produced by the flux vithin this 
half turn. Similarly the third term is the power due 
to the current, [,”, below the lower half turn and the 
pressure produced by the flux within this half turn. 
For the q conductor, I,’ = I(q—1+n/0) and 
L” = I (q— 1). We have also shown that in this 
case (2) I, = I ؟)‎ +n/2/ 0), and L -1(q—1 
+ n/2 / 0). ۰ 

As before 0 is the phase “ni between the currents 
in the upper and lower coil sides lying in the same 
slot. Making these substitutions and dividing by the 
direct-current resistance loss for a whole turn gives 
a resistance ratio per turn of 


m (M ,' 十 (q? — q 十 n? /4 + (2q 一 1) n/ 2 COS 0) N,’ 
+ 2/4 T,’) 


The average value of this resistance ratio for a 
whole coil of n turns is 


K =| me + مر سح‎ 


+ 4 2d 


Case 3. Strands joined at the beginning and end of 
à whole turn, end connections turned over. The heat- 
ing loss for a whole turn is now: 
los -271.R (I M' + ۲0/۵2 N' + I, S'/2) 
+I. R ( — (1/2 +I) ۸۷۲ (1: + b + D ۷ 
+I S'/2 
+ 1." RI(I/2 +I) ۸۷۲ + (L, — Ip) T’ 
+ I S”/2) 


Feb. 1922 


+I” . RI — (1/2 1 N' + (L,” + L + D T 
+ 158'/2) ۱ 

' The first term is the power due to the current in the 
turn and the resistance drop in the top strand of the 
turn in the upper coil side and the bottom strand of 
the turn in the lower coil side. Due to the turning 
over of the end connections these two half-turn strands 
are a part of the same strand. The second term 15 the 
power due to the current in the turn and the pressure 
produced by the flux within the half turn that is in 
. the lower coil side. The third term is the power due 
to the current below the half turn in the upper coil 
side and the pressure acting on this current that is 
produced by the flux within this half turn. The 
fourth term is similarly the power due to the current 
below the half turn in the lower coil side and the pres- 
sure acting on this current that is produced by the 
flux within this half turn. İn this case: 


LIL = (— 1/2 +n/2/0)I "و‎ =(q-1lt+n/6)I 
I, = (q — 1/2 +m/2 / 0) I L,” = (q— 1) I 

Making these substitutions the resistance ratio for 
a whole turn reduces to: 


n 一 1 


K = 1 0, + لحك‎ wy’ 


+ [ (2 q — 1? + + 2 (2 q — 1) n cos 0] 4 | 


The average resistance ratio for a whole coil of n 
turns 1s: 


1 n — 1 


, n— 1 , 


+ n?/2 cos 0) T. | 


The method of calculating the leakage impedance 
when the stranding is continuous throughout a whole 
coil is described in considerable detail in the preceding 
paper and need not be repeated here. With a finite 
number of strands, however, there are added terms 
involving S’ which appear on account of the insulation 
between the strands. Fortunately their effect is not 
difficult to calculate. f 

The added resistance drop in a coil of n turns is 
(equation 19)2nRI,/2 S”. Values of I, have already 
been calculated for the different cases. The added 
pressure acting in the coil side which lies in the bottom 
of the slot due to flux within the coil side above it is 
n R (n /— 0) I/2 S'. The other added terms are 


those due to the pressure acting in the eonduetors of 
the coil produced by flux within the coil itself. There 
are three cases, vz., 4, 5 and 6. 

Case 4. When the end connections are not turned 
over the resistance drop taken is that in the upper 
element of the top strand of the conductors of both 
coil sides. The added drop is then: 
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Case 5. When the end connections are turned over 
on one side only the resistance drop taken is that in 
the top element of the top strand of the first half turn 
plus that in the bottom element of the bottom strand 
Of the next half turn plus that in the bottom element 
of the bottom strand of the next half turn plus that in 
the top element of the top strand of the next half turn 
plus, ete. İn this case it is readily shown that the 
added drop is: 


È (2q— 1) R 1/2 S” or 72/2 RIS’ 
I 


Case 6. When the end connections are turned over 
on both sides the resistance drop taken is that in the 
top elements of the top strands of one coil side and the 
bottom elements of the bottom strands of the other 
coil side. The added pressure is thus 


2 (q — 1( 81/2 5 + ZqR 1/2 S") or 7۵/2 7 ۲ ۲ 


The sum of these three eomponent added terms is 
the same in each of the three cases. Itis (n? + n? cos 6) 
RIS’. 

The average value for a single half turn is 
(n/2 + n/2cos 0) RIS’. ۱ 

We may now write the expressions for the slot leak- 
age impedance of a symmetrical pair of fraetional 
pitch slots. 

Refer to the preceding paper. 

Case 4. End connections not turned over. 


Z -R| m + (21 + wos) Nr 


2 
+ 5 “15 1 p: + (n/2 + n/2 cos 6, s 
where R. is the true resistance of a whole coil. 


Case 5. End connections turned over on one side 
only. Even number of conductors per coil side. 


10m — 1 


19 ( + سا 7 


+ "?/2 cos 0) T’ + (n/2 + n/2 cos 0) s | 


Case 5. Odd number of conductors per coil side. 


10 22 — 1 


7 -R| a” + ) - + 8/2 دہ‎ ) ۳ 


+ (n/2 + n/2 cos 0) s| 


Case 6. End connections turned over on both sides. 


xn, 772 1‏ 5-1و ; m‏ 
72( تست + M‏ اھ و 
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x ; . Bd _ 0.0245 / 89° 38.5” 
+ n?/2 cos 0) T' + (n/2 + n/2 cos 0) S | sinh Dn N L97— 0.827 11 105 
Since S” is pure imaginary, the terms involving it = y 0.01164 / 45”1.5” 
do not appear in the expressions for the alternating- = 0.07628 + j 0.07634 
current resistance. In the preceding paper a? d? ۱ 
which is now replaced by T” was pure imaginary. Let ge =gt+jk 


T', however, has both real and imaginary parts, and 
thus adds to the value of the alternating-current 
resistance. 

There follows a numerical caleulation of the heat 
losses in a specified winding. The pitch of the coils is 
one and the dimensions of the slot and the conductors 


are: 
Width of slot (S) .................... = 2.54 cm. 
Width of conductor (w) .............. = 1.60 " 
Length of armature core (I) ........... = 72.4 ë 
Length of end turn (L) .............. = 80.0 s 
Depth of strand (d/m) .............. = 0.254 “ 
Thiekness of insulation between strands 
(人 = 0.0381 “ 
Number of strands per conductor (m) .. = 7.0 
Number of conductors per coil side (n) . .= 2.0 
Frequency (f) ...................... = 60 cycles 
Average temperature of winding ...... = 100 deg.cent 
p = 2260 c. g. s. ohms at 100 deg. cent. 
a = 27 | 2907 /45° = 1.15 / 45° 
ps 
ad -= 0.292 / 45° 
mo C (一 一 
ad R 
2m 0.146 / 45 
全 = 0.0219 / 45° 
tanh 2% — 0.146 / ۰ 
2m 1 
ad 
m 
= 1.00 / — 0.824 
sinh & d 0202 
m 
T ad aa a d ) 
4 m? 2m ` ( 2 2m 
= 0.0245 / 89° 8.5“ | 
From equation (10a) ` ۱ 
ad { aa ad 
un B4 _ Xm ) و‎ + tanh 25) 
ده‎ fe تھا‎ e 
m 2 
+ "EB 
sinh 24. ^ 


4. From Kennelly’s “Tables of Complex Hyperbolic and 


Cireular Functions." 


sinh (g + 7 k) = sinh g cos k + j cosh g sin k 
= A +78 = C /tanh-i B/A 
That is sinh g 609 =A 
| eoshgsink = B 
Solving these two equations for g and k gives: 


-ü-C)-tvü-Cy-r4A 
2 


sinhg — ۷ 
and sink -MV 


Substitution in these solutions shows that: 


1 + C2— V/(1 + C*)?— 4 B° 
2 


sinh g = 0.0766 g = 0.0767 
snk = 0.0761 k = 0.0762 radians 
Therefore BE — 0.108 / 44° 507 
tanh pi is readily computed in this case from 
sinh 8d 
tanh B4. 2m 


2m J., BA 
)م۸۷‎ + sinh2 Be 
0.1079 / 57 
` V1 + 0.01164 / 3' 
= 0.1079 / 44° 21.5” | 


ad 
= + tanh ہچ‎ 


tanh Bd 
2m 


a d = 8.18 / 45? 19' 


The values of coth 8 d and 2 tanh — 3 q can be com- 


puted from the formulas: 

coth Bd 
_ sinh 2 mg cosh 2m g — ?sin 2m k eos 2 m k 
— sinh? 2 m g cos? 2 m k + cosh? 2 m g sin?2 m k 


2tanh £2 


sinh m و‎ cosh m g + j sin m k cos m k 
cosh? m g cos? m k + sinh? m g sin? m k 
From which 
coth Bd = 0.8835 / — 11° 20” 


2 tanh -EÉ = 1.472 ۰ 
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The complex quantities M”, N” and 7” may nov be 
calculated 
M" = 1.335 / 33°59 1.11 + 70.746 
N’ = 2.22 / 79726” = 0.408 +72.19 
T’ = 2.28 /89° 38.5’ = 0.014 + j 2.28 
If the number of strands, m, is inereased without 
limit while the depth of the conductor and the relative 
amount of insulation between the strands are unchanged 
we have: 


— 


ey Bd. Bd . ad a/d +1 
sinh oe سل‎ EO VUL 


mn l5/li + 1 
: od . ad 
since tanh-5- "m 
ad 
and — 1ہ‎ 
sinh 


M" ù Bdcoth Bd 
N’ ù 842 tanh کے‎ 


For the winding we are considering, 
Bd -0.789ad =1.51/ 45” 
. For this value of 8d; M,’ = 1.11 and N,’ = 0.418. 
In this case, since T’ is pure imaginary, T,’ = 0. 
A further calculation has been made for conductors 


of the same net cross-section, but consisting of three 
strands instead of seven. In this case: 


M" = 1.295 /38۰ 9۶ M, = 1.084 
N’ 22.04 /80°52’ N,’ = 0.824 
T' =2.086/87°52’ T,’ =0.0776 


The following is a table of the ratios of alternating 
to direct-current resistance for this winding for each 
of the six strand arrangements. 


Strand Resistance Ratio 
Arrangement m =a m =7 m =3 
Oase 1 3.17 3.15 2.70 
v 2 2.76 2.75 2.46 
8 1.42 1.47 1.66 
“ 4 2.66 2.66 2.40 
“ 5 , 1.01 1.08 1.41 
* 6 1,42 1.47 1.66 
Discussion 


W. V. Lyon (by letter): A series of curves has been plotted 
which, though they may not be of partieular value themselves, 
indicate a line of investigation that should prove of eonsiderable 
importanee. Curves similar to them have been discussed by 
Rogowski and others. 

The curves marked “T” are for the upper bar of a laminated 
bar winding, with two bars per slot. 

The curves marked “II” are the average values for the slot. 
Laminated conductors with laminations joined at the beginning 
and end of each turn. Two turns per coil, two coil sides per 
slot. 
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The curves marked “III” are the average values for the slot. 
The winding is exactly like the preceding case but the end con- 
nections are turned over. | | 

The curves marked “IV” are the average values for the slot. 
Laminated coil with laminations joined at the beginning and 
end of the coil. Two turns per coil, two coil sides per slot. 
End connections turned over. 

All of these curves are plotted for full-pitch, finely laminated 
windings. The length of the armature core is assumed equal 
to the length of one end connection. The ordinates of these 
curves show relative values only. The abscissas, « d, are roughly 
equal to the depth of the conductor in centimeters for a fre- 
quency of 60 cycles. 
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Ra MIH CONSTANT WIDTH OF CONDUCTOR OR WATTS 
PER SQ. CM OF SIDE OF CONDUCTOR WITH CONSTANT 


WATTS PER SQ.CM.OF SIDE OF CONDUCTOR. 
WIDTH OF CONDUCTOR CONSTANT 


Fig. 1 shows how entirely inappropriate laminated bar wind- 
ings are when the conductors are deeper than about one centi- 
meter. Fig. 2 shows that with conductors of constant width 
there is a depth which makes the alternating-current resistance 
a minimum. This has been called the critical depth. Notice 
that with the better types of windings, viz. III and IV, the 
critical depth is greater. The ordinates of those curves may 
also represent the watts per square centimeter of coil side for 
conductors of constant cross section. If all of the heat developed 
in the conductors passed through the sides of the coil, which of 
course it does not, these curves would show a truly critical 
depth. 
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If the width of the conductor is kept constant the watts per 
square centimeter of coil side do not reach as distinct a mini- 
mum value. 
is no minimum. . e., critical depth. Fig. 3. ۱ 

Taking the surface through which the heat is conducted from 
the coil as the entire perimeter of the cross section of the 
coil the watts per square centimeter of insulation are plotted in 
Fig. 4. These curves were plotted for cases in which the cross 
section of the conductors was nine square centimeters and the 
thickness of the insulation was one-half a centimeter. 

It seems that these curves only emphasize the need of a 
thorough investigation of the heat conduction through the 
insulation of embedded conductors. The problem is further 
complicated in that the best shape of conductor must be deter- 
mined by considering not only the heat generated within it 
but also by the eddy current and hysteresis losses produced in 
the neighboring iron. 


WAITS FER SQ CM OF INSULATION. CON 
DUCTOR HAVING CONSTANT CROSS SECTION 
CONSTANT THICKNESS OF INSULATION 


M. S. Vallarta (by letter): Although a few resistance-ratio 
and temperature measurements have occasionally been reported, 
no thorough experimental investigation of the Field theory of 
skin-effect in embedded conductors has to our knowledge ever 
been undertaken. In order to fill this gap, such an investiga- 
tion is now in progress at the Research Division of the Electrical 
Engineering Department, Massachusetts Institute of Tech- 
nology. Its main purpose is to furnish experimental proof of 
the correctness of Field’s assumptions, since the rest of the theory 
and its applications, as developed by Rogowski, Richter, Gilman 
and very recently by Professor Lyon, in his simple method 
involving the hyperbolic functions of a complex variable. follow 
from these assumptions by a process of mathematical deduction. 
Owing to experimental difficulties it has been found necessary, 
however, to test both the assumptions and the conclusions, as 
a check on the complete theory. 

A piece of standard laminated slotted armature was kindly 
furnished for these tests by the Westinghouse Company, to 
which we acknowledge our indebtedness. This is mounted on a 
wooden framework, away from any magnetic material. So far 
as the magnetic circuit is concerned, conditions are in close 
agreement with the ideal demands of the theory. - 

For the purpose of obtaining the greatest possible resistance, 
the test coils are made of thin copper strip wound longitudinally 
in the slots, with its largest. dimension parallel to the slot side 
and with paper-insulated turns. 
skin-effeet is concerned, which is the only one considered in the 
Field theory, these coils are the exact equivalent of a solid bar 
conductor of the same dimensions. Stranded conductors with 


twisted end connections have not been thus far experimentally — 


investigated. ۱ 

To test the assumption that an element of current in the slot 
produces no field below it, an exploring coil was placed longi- 
tudinally in the slot, with its plane parallel to the slot side, 
below a current carrying conductor. This exploring coil was 
connected through a shielded two-step vacuum tube amplifying 
circuit and thermocouple to a sensitive galvanometer. No 
deflection could in any case be obtained, thus confirming the 
assumption. 
no component of field strength parallel to the slot side has been 
found. 
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In fact in the better windings, III and IV there . 


So far as one-dimensional 


No satisfactory proof of Field’s assumption of . 
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An experimental curve of cross-flux distribution was deter- 
mined by winding five equi-distant exploring coils around 
the tooth and measuring the voltage induced in them when an 
alternating current flows through the slot conductors. This 
voltage was measured by connecting the exploring coils in turn 
to a sensitive galvanometer through a shielded two-step vacuum 
tube amplifying circuit and thermocouple. To eliminate 
tooth-tip leakage, an exploring coil at the top of the slot was 


always connected in series opposition with the coil through 


which the flux was measured. 

The eurve of cross-flux distribution is easily calculated from 
equation four of Professor Lyon’s first paper. The flux linking 
an exploring coil distance z from the bottom of the upper con- 


ductor is 
d x | 
s= | ات‎ Lİ veda h İdə 
x 0 


where c = f (z) is given by the equation quoted above, 2. e., 


c =a/w pa E Soak az + İşsinh az | 
sinh ad 2 


It follows by integration, if 7; = ول‎ = I. - 
4 «ll 


o 5 


cosh a z 


d 
dı = ۱ coth a d 一 — tanh E (cosh a d 


sinh « d 


— cosh a z) + sinh a d — sinh «s | ۱ 
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For the lower conductor alone: 
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Let #19 be the total flux linking an exploring coil at the bottom 


al 
6 = س‎ 
w 
whence 

4 «lI 


as 


coth aa —‏ | = ره 
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of the upper conductor. The total flux linking an exploring 
coil distance z from the bottom of the slot is: 


وش + :$1 = $ 

From which the flux at any point along the lower conductor 
can be calculated. 

The experimental and calculated curves of cross-flux distri- 
bution are given in Fig. 5. The agreement is considered satis- 
factory, within the limits of experimental error. 

Resistance-ratio measurements have also been made. Power 
is measured by the three-voltmeter method, especial precautions 
being taken to have a sinusoidal voltage wave. Oscillographic 
records show that this condition was in every case closely 
fulfilled. The disturbing effect of uneven temperatu.e dis- 
tribution is largely eliminated by operating only after a steady 
uniform temperature, such as caused by a direct current, has 
been reached, then making all a-c. measurements within a short 
interval. To correct for iron loss, a coil of fine magnet wire was 
wound around the teeth, the excitation being distributed so as to 
correspond as closely as possible with conditions in the slot con- 
ductors. Phase difference effects of course cannot be imitated 
with this arrangement. The skin effect of such a winding is 
assumed to be negligible. 

Results of resistance-ratio measurements and computation 
data are given below. Each. power measurement given is the 
average of at least five runs. The precision of individual meas- 
urements is estimated at one per cent or better. 


COMPUTATION DATA 


Dimensions of slot: VVidth..................... 1.905 cm 
Length..................... 40.6 cm 
Depth..................... 7.62 cm 

Effective width of conductor.................... 1.418 cm 

Effective height of conductor.................... 3.58 cm 

Erequency................................. 60.0 cycles 
Temperature.............................. 65 deg. cent 
TABLET 
Measurements for two con- 
ductors in series 
i 
Correction Correction 
Measured for end for iron Net A-C. 
power Current turns loss power resistance 
watts amperes watts 1 watts watts ohms 
204.9 9.14 3.84 6.05 . 195.0 2.34 
68.0 5.21 1.22 2.21 64.6 2.38 
11.53 2.131 0.20 0.50 10.83 2.39 
3.32 1.136 0.05 0.22 . 3.05 2.38 
Average 2.37 ohms 
Embedded portion d-c. resistance..... 0.2251 ohms 
Resistance ratio.................... 10.52 
Calculated resistance ratio........... 10.26 
Difference................... .. ۳ 2.5% 
TABLE II 2 
[Measurements for lower conductor 
Correction Correction 
Measured for end for iron Net A-C. 
power Current turns loss power resistance 
watts amperes watts ! watts watts ohms 
95.8 15.16 5.25 3.14 87.4 0.382 
34.9 9.16 1.92 1.21 81.7 0.378 
11.81 5.32 0.65 0.37 10.79 0.381 
1.98 2.21 0.09 0.05 1.84 0.377 
Average........... 0.3795 ohms 
Embedded portion d-c. resistance..... 0.1148 ohms 
Resistance ratio...................., 3.30 
Calculated resistance ratio............ 3.25 
Difference.............. ere دب‎ 2.0% 


1. Length of end turns is 20 per cent of total length. Correction as- 
sumes that skin effect in air for the coils is negligible, which is borne out by 
experiment. 

2. Lower conductor alone in slot. 
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TABLE III š 
Measurements for upper conductor 
Correction |Correction | ۱ 
Measured for end for iron ə Net A-O. 
power Current turns . loss power resistance 
watts amperes watts 1 watts watts ohms 
181.5 9.36 2.15 2.36 177.0 2.04 
62.7 5.47 0.74 1.54 60.4 2.01 
10.17 2.145 0.11 0.46 9.60 1.98 
2.70 1.140 0.03 0.07 2.60 2.00 
Average........... 2.01 ohms 
Embedded portion d-c. resistance..... 0.1110 ohms 
Resistance ratio..................... 18.1 
Calculated resistance ratio........... 17.27 
Difference... ə تب پٹ‎ R au eR 4.8% 


3. Both conductors carrying the same current. 


It is seen from the above that, as far as can be judged from our 
present evidence, the complex hyperbolic method developed by 
Prof. Lyon correctly describes the eross-flux distribution and 
gives simple means of computing the resistance ratio of ideal 
bar windings with engineering accuracy. 

In connection with this investigation, the following biblio- 
graphy has been prepared: 


BIBLIOGRAPHY 
(In chronological order and without pretensions as to com- 
pleteness. ) 

1. A. B. Field. ‘‘Eddy-Currents in Large Slot-VVound 
Conductors” Trans. A. I. E. E. Vol. XIV, 1905, p. 761. 

2. Fritz Emde. "'Einseitige Stromverdrangung in Anker- 
nuten" Eleki. und Masch. XXVI, 703 & 726, 1908. 

3. W. Rogowski. Uber zusátzliehe Kupferverluste, über 
das kritische Kupferhohe einer Nut, ete". Arch. fur Elek. I1, 81, 
1913. 

4. R. Rüdenberg. “Die kritische Kupferhöhe von Leitern 
für Dynamomaschinen’’. Arch. für Elek. II, 207, 1913. 

5. L. Dreyfus. ‘‘Zusatzliche Kommutierungsverluste bei 
Gleiehstrommasehinen." Elekt. und Masch. XXXII. 281, 1914. 

6. L. Fleischmann. “Uber zisatzliche Verluste durch 
Stromverdràngung in Gleiehstrommasehinen.” Arch. für Elek. 
II, 387, 1914. 


7. R. Richter. “Uber zusätzliche Kupferverluste.” ^ Arch. 


für Elek. II, 518, 1914. 


8. W. Rogowski. “Uber zusätzliche 
Arch. für Elek. II, 526, 1914. 

9. R. Richter. “Uber zusätzliche Stromwārme. I: Einfluss 
der Leitfähigkeit, ete.” Arch. für Elek. III, 175, 1914. 

10. F. Hillebrand. “Uber zusätzliche Kupferverluste bei 
Wechselstrommaschinen.” Arch. für Elek. III, 111, 1915. 

11. L. Dreyfus. “Die Theorie der zusātzlichefn Kommu- 
tierungsverluste von Gleiehstrommasehinen." Arch. für Elek. 
III, 273, 1915. 

12. R. Richter. 
wurf von Nutenwicklungen.”’ 


Kupfervderluste.” 


“Uber zusatzliche Stromwarme. II: Ent- 
Arch. für Elek. IV, 1, 1915. 


13. L. Dreyfus. ‘‘Zusatzliche Kupferverluste durch Strom- 
verdrangung bei Einankermumformern.” Arch. fir Elek. IV, 
42. 1915. 


14. W. Rogowski. ‘‘Unterteilung und Wechselstromwider- 
stand.” Arch. fiir Elek. IV, 279, 1916. 

15. R. Richter. ‘Uber zusatzliche Stromwarme. ITI: Nuten- 
wicklungen mit unterteilten Leitern." Arch. für Elek. V, 1, 1916. 

16. L. Roebel. German Patent 277012. 

17. R. Richter. “Uber zusützliehe Stromwarme. IV: 
Unterdrüekung dureh magentische Hilfskreise." Arch. für Elek. 
V. 335, 1917. 

18. W. V. Lyon. “Current Distribution in Armature Con- 
ductors.” Bulletin 19. Res. Div. Elec. Eng. Dept. M. I. T. 1919. 

19. A. E. Gilman. ‘‘Eddy-Current Losses in Armature 
Conduetors." Trans. A. I. E. E., Vol. XXXIX, 1920, Part I, 
p. 997. 
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20. H. W. Taylor. 'Eddy-Currents in Stator Windings.” 
Jour. I. E. E., LVIII, 279, 1920. 

21. R. Pohl. “Eline einfache Theorie der zusatzliehen Ver- 
luste im Nutenkuper von Wechselstrommaschinen." Elekt. Zeit. 
XLI, 908 and 997, 1920. 

22. W. V. Lyon. “Heat Losses in the Conduetors of A-C. 
Machines." Trans. A. I. E. E., Vol. XL, 1921, p. 1361. 

23. W. V. Lyon. “Heat Losses in Stranded Armature 
Conductors.” Jour. A. I. E. E., Vol. XLI, 1922, Jan., p. 37. 


V. Karapetoff: Prof. Lyon deserves much eredit for having 
worked out in detail various expressions for heat losses in arma- 
ture eonduetors, with à non-uniform distribution of eurrent. I 
wish only to point out some possible improvements in the funda- 
mental mathematical treatment, and to indicate the next steps 
‘that might be taken in the solution of this important problem. 


1. Gilman's Work. In June 1920 Mr. R. E. Gilman presented 
before the Institute & highly mathematieal paper on the same 
‘subject: (Trans. Vol 39, p. 997), and it is to be regretted that in 
Prof. Lyon's paper, presented a year and a half later, before the 
same body, no attempt is made to correlate his work with that 
of Gilman's, and to point out identical results, discrepancies, 
if any, advantages and disadvantages of the two different methods 
of mathematieal analysis. It remains for a future investigator 
to finish this part of the work. 


2. Exponential Notation. In my discussion of Gilman's paper 
(ibid, p. 1054) the advantages of exponential notation are 
pointed out, and it is shown that several pages of long tedious 
formulas may be done away with, and that a solution of certain 
simultaneous equations is obviated. It now remains to compare 
the method of exponential notation with the use of hyperbolic 
functions of a complex variable, as used in Prof. Lyon’s paper. 
In making this comparison, it is necessary to keep in mind that 
at the present time we have only Kennelly’s tables of such func- 
tions, with steps of such magnitude that the two-directional inter- 
polation is quite tedious. The building up of hyperbolic func- 
tions of a complex variable out of those of real variables also 
takes considerable time. 


3. The fundamental differential equation. The fundamental 
differential equation of distribution of alternating-current density 
in a long conductor, subjected to a transverse magnetic flux, is 
not original either with Mr. Gilman or with Prof. Lyon. See, 
for example, A. Russell, Alternating Currents, Vol. I, index 
under “eddy currents.” This equation can now be derived much 
more directly from Heaviside’s laws of circuitation which are 
familiar to the younger generation of engineers. We want as 
wide a circle of readers for our Institute papers as possible, and 
any simplification in mathematics, any deduction directly 
from a general physical law, rather than by a special “follow me” 
method, is a step in the right direction. 


4. Some mathematical simplifications. Prof. Lyon’s paper 
being based on a well-known fundamental equation, the value of 
the contribution lies mainly in the application of the solution 
to certain specific cases. The solution must be in the simplest 
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possible form for numerical computations. An inspection of 
his formulas shows the possibility of considerable further sim- 
plification. The formulas used in the present paper are based on 
equation (3) in his first paper (TRANS., 1921, Vol. 40, p. 1361). and 
it is therefore necessary to go back to that paper in order to 
indicate an alternative treatment. 

The objections to Prof. Lyon's formulas are (a) that 
they are involved and unsymmetrical, and (b) that the 
variable z enters in the same formula in two or more places. 
What seems to be a simpler form of solution both for a general 
development and for numerical work, is indicated below; the 
equations numbers 1, 2 and 3, refer to Lyon’s first paper: 

Let the general solution of equation (2) be written in the 
form of ۱ 


c = (a/u) D Sinh (ag + 8) (10) 


This solution differs from equation (3) in the following re- 
526615 (a) the expression (a/w) is written out explicitly in order to 
simplify further transformations; (b) the variable distance, y, 
is assumed to be measured from the center of the conductor, 
in order to make the equations more symmetrical. In the origi- 
nal paper the corresponding variable, z, is measured from the 
lover edge; (e) the variable y enters in the equation only once, 
while equation (3) contains 2 twice; (d) The two constants of 
integration are D and a, that is, directly the amplitude and the 
phase, instead of the components A and B. 


To determine D and 8, we substitute expression (10) in equation 
(1), and in the formula for 7,*. After simplification we find: 


tanh م‎ = [(Z1/(11 + 210] Coth ad/2 (11) 


D = 1,/(Sinh م‎ Sinh a d/2) (12) 


v|- 


Equation (10), with the auxiliary expressions (11) and (12) leads 
to simpler formulas and computations than the formulas used in 
both of Prof. Lyon's papers. For the current density at the 
center of the conductor we have, putting y = 0, 


c = (@/w) D Sinh s (13) 
or, substituting for D its value from eq. (12), 
c = a 1,/(2 w Sinh ad/2) (14) 


which checks with Prof. Lyon's expression. 


5. Future Work. There is no reason why this problem should 
not be brought now to a final solution, at least in application 
to standard turbo-alternators. By this I mean a set of curves, 
charts, tables, ete. with which a designer could safely compute 
the armature copper loss for an assumed arrangement of con- 
ductors. At the present stage the designer would have to study 
one or two highly mathematical articles and then perform rather 
long computations before he could get a specific result. Should 
he then decide to change his winding, much of the numerical 
work would have to be repeated. 


“TRANS. A. I. E. E., 1921, p. 1369. 
Trans. A. I. E. E., 1921, p. 1372. 
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Another Harmonic Analyzer 
BY F. S. DELLENBAUGH, Jr., 


: Member A. I. E. E. 
Secretary Research Division, Electrical Engineering Department 
Massachusetts Institute of Technology. 


Review of the Subject.— Modern alternating-current engi- 
neering involves the necessity of frequently analyzing curves for 
harmonics. All of the methods of analysis have more or less serious 
objections. A method based upon the averaging of ordinates 
selected at the proper intervals has been developed into a machine 
which is reasonably simple and cheap to construct and can be 
operated without great skill. The machine depends chiefly for its 
operation upon a bar containing notches at the proper intervals and 


each bar is suitable for curves of the same plotted wave length.. 


tory at Massachusetts Institute of Technology, 

Mr. David O. Woodbury devised an extremely 
simple and useful harmonic analyzer, which is thought 
worth describing as it can be made to analyze har- 
monics of very high order without the usual attendant 
exhaustion suffered by the investigator. It is also 
rapid for harmonics of any order, simple to build, and 
while it is limited to curves having the same length of 
base for one period, it can be fitted to any size of curve 
by merely making a new bar, which must be accurately 
machined, but is not complicated or expensive. 

The theory of harmonic analysis has been gone into 
very thoroughly in many text books and will not be 
here repeated. For bibliography the reader is 
referred to the A. I. E. E. JOURNAL, Vol. XL, p. 135, 
1921. The analyzer being described utilizes the 
method of selected ordinates, originally proposed by 
Houston-Kennelly, and developed in detail by Fischer- 
Hinnen. 
number of ‘portions for each harmonic, the number of 
parts being the same as the order of the harmonic, 
In this way the sine or cosine multipliers of the usual 
schedule form become unity (or zero) and the value of 
the coefficient for the particular harmonic may be ob- 
tained by simply averaging ordinates at the different 
points given by the division of the curve into parts. 
This may be expressed as follows: 


Wi ILE working in the Research Division labora- 


For Opp HARMONICS ONLY (180. DEG. OF CURVE 
REQUIRED) 


If, starting at X = 0, we measure ordinates at 
intervals of m/n, the average of these ordinates taken 
alternately plus and minus is equal to the sum of the 
amplitudes of the nth, 3nth, 5nth............ Cosine 
components. 

If, starting at X = 7۳/2, we measure n ordinates 
at intervals of m/n the average of these ordinates 
taken alternately plus and minus is equal to the sum 
of the amplitudes, taken alternately plus and minus, 
of the nth, 3nth, 5nth........ Sine components. 


The curve must be divided into a different | 


It ts thus limited to curves plotted with the same scale of abscissas, 
but different bars may be manufactured for different scales and the 
machine thus adopted to various types of curves. A comparison 
with the times required for harmonic analysis by different methods 
shows that it compares favorably with other devices. 


CONTENTS 
Introduction (430 w.) 


Construction of Analyzer (540 w.) 
Operation of Analyzer (800 w.) 


For BOTH ODD AND EVEN HARMONICS. (860 DEG. 


OF CURVE REQUIRED) 


If, starting at X = 0, we measure 2n ordinates at 
intervals of m/n, the average of these ordinates taken 
alternately plus and minus is equal to the sum of the 
amplitudes of the nth, 3nth, 5nth,......... . Cosine 
components. E. | 

If, starting at X = 7۳/2, we measure 2n ordinates 
at intervals of m/n, the average of these ordinates 
taken alternately plus and minus is equal to the sum 
of the amplitudes, taken alternately plus and minus, 
of the nth, 3rith, 5nth........ Sine components. 


CONSTRUCTION OF ANALYZER 


The analyzer consists in principle of a bar which is 
free to move in the direction of ordinates but fixed in 
the direction of abscissas. This bar contains a series 
of notches spaced 1/8rd, 1/5th or any interval of 
r /n for a distance equal to the length of curve necessary 
to explore for analysis. By rotating the bar any series 
of notches may be brought uppermost. An index 
carrying cross lines slides upon this bar and bears a 
spring-operated dog, which drops into the notches 
sufficiently to locate the position of the index, but 
lightly enough so that it may be pushed along without 
effort. The.carriage upon which the bar is mounted 
has a magnetic clutch which engages a scale, the motion 
being imparted to the scale by the magnetic clutch, 
and a brake preventing motion of the scale when the 
circuit of the magnetic clutch is open. If the curve to 
be analyzed is now placed upon a platen under the 
Index the procedure given in the rules quoted above 
may be followed mechanically by merely moving the 
index from notch to notch, bringing the cross lines into 
alinement with the curve for each position, and properly 
operating the magnetic clutch so that the respective 
ordinates are alternately added and subtracted upon 
the scale. The scale reading will then give the sum of 
the nth, 3nth, 5nth, etc. harmonic coefficients. 

The mechanism may obviously be made either to 
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have a sliding motion in a plane or rotate around 
centers. A sketch of an early model operating upon 
. the first principle is shown in Fig. 1. The scale is a 
steel rule laid across the top. The magnetic clutch is 
merely an old bell magnet running on two round rods. 
This particular machine was arranged to take a card- 
board template of the curve, the motion being imparted 
to the clutch by a small roller following the template. 
The-clutch was operated by a contact roller traveling 


Fie. 1—Earty MopEL or WOODBURY ANALYZER 


The template under roller at right hand side of base is changed for each 
Bara 


over a specially prepared plate with alternate spaces 
of insulation and metallic conductor, which took the 
place of the notched bar. Different plates were used 
for each order of harmonic. This machine was in- 
tended to be practically automatic so that the roller was 
dragged over the curve template slowly and the 
coefficient read off the scale after each traverse. Me- 
chanical difficulties made this type impractical, so an 
improved model was made, the parts swinging around 
centers. Thisi is Shown in Figs. 2 and 3. 


Fic. 2—FRoNT View or FINAL MODEL or WOODBURY ANA- 
LYZER 


The scale upon which the coefficients are read may be seen abutting 


the left hand end of cylindrical platen. 


The curve is clamped upon the cylindrical platen 
by means of the two bands. The notching bar may be 
seen with the setting head upon the right hand end, Fig. 


2. The index with cross lines does not show plainly . 


but is carried upon the forward part of the carriage 
riding upon the notching bar and plain guide bar. 
magnetic clutch is attached to the left hand radius 
supporting arm, showing in right hand end of Fig. 8, 


The. 


and the scale is adiacent to the left hand end of the 
cylindrical platen almost in contact with clutch. The 
index for the scale is attached to the left hand side of 


platen, Fig. 2, and the brake is within the box and may 


be. adjusted by the thumb screw shown in Fig. 3 under 
right hand strap. The clutch is operated by a switch 
In the handle upon one left hand radial supporting arm, 
Fig. 2. 


Fic. 3—Rear View or FINAL MopELOF WOODBURY 
ANALYZER 
The small cylindrical drum at the right is the magnetic clutch. The 
thumb screw under right hand strap is the brake adjustment to prevent 


l SUD of scale. 


OPERATION OF ANALYZER 
To perform the analysis, a notehing bar matching 


the curve to be analyzed is placed in the machine. 


If the curves are many different scales this, of course, 


is a serious objection as they must be reduced to the 


same scale, but if oscillograms are made with a syn- 
ehronous film drive, or if the curves are plotted by 
hand from tabular data, the scale can be always made 
the same and therefore one bar is sufficient unless 
very many harmonies are desired. The development 
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Fic. 4—DEVELOPMENTS or NOTCHING Bar ror Opp HARMONICS 


Up TO 7 


of the notching for a bar to be used in analyzing a curve 
with a scale of 6 in. = 180 deg. covering odd harmonics 
from the 8rd to the 21st is shown in Fig. 4, and the bar 
may be constructed out of half inch stock. 

The curve is next adjusted so that its axis is parallel 
to the notehing bar, readily determined by sliding the 
index with the bar fixed, and the index adjusted so that 
when the dog is in the zero notch, the cross lines cor- 
respond with the zero point on the curve. The index 
is mounted upon a sleeve sliding within the carriage 
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for this purpose, which is fastened by a set serew when 
in the right place. The setting head is then turned 
to the harmonic desired and the carriage moved until 
‘the dog falls in the first notch to right of zero. The 
index is then registered with the curve, the magnetic 
clutch switch closed, the index then moved until the 
dog falls into the next notch, registered with the curve, 
and the clutch switch released. The carriage is then 
moved to the next notch and the same sequence re- 
peated. With a small amount of practise this manual 
sequence of events becomes practically automatic and 
the analysis will proceed as swiftly as it is possible to 
slide the carriage from notch to notch. The reading 
upon the scale will then be the sum of the amplitudes 
of the nth, 3nth, 5nth, etc. Cosine terms. 


0 36 72 108 144 
SINE COMPONENT, DEGREES 
Fie. 5—Tar Heavy LINE SHOWS THE TRAVEL OF ANALYZER 
INDEX. 
During the solid portions the clutch switch is closed and during the 
dotted portions the clutch is released. Any path may be followed by the 


index between the points of the curve A. B. C. etc. These points are 
located by the notching bar. 


The scale of the result depends upon the divi- 
sions made upon machine and reads direct in inches 
or centimeters. For the sine terms the carriage 
is moved to the special setting notch shown in Fig. 4, to 
the left of the zero notch, the index set screw is loosened, 
the index set to the zero of the curve, the carriage 
. moved to the zero notch, and then a similar procedure 
to that for the cosine term followed. The motion of 
the index and the closing of the clutch switch for the 
9th harmonic is indicated diagrammatically in Fig. 5. 
The results can be figured in per cent of fundamental, 
or if scale of curve is known, they can be readily 
transferred to terms of the function represented. 

The fundamental is obtained from the values of the 
harmonics and from a few specially read ordinates. 
The following equations show how this is done: 


BoTH ODD AND EVEN HARMONICS 
Y, = bo + b, + b; + b; e ə ee و‎ ' 
Hence bo = Y, — (b, + b. + bs ۷0ئ0‎ ) ۱ 
1/2 (Y, 一 Y 180) = b, + b; + bs e e o es 


Hence b, — 1/2 02 — Y 180) 一 (b; + b; تہ‎ .) 
1/2 (Yo 一 Y 270) = 0 — Gs + و0‎ — Q7 ew ee 
Hence a, = 1/2 (Ys, 一 Yzo) 一 (— Gs + و0‎ - Ga...) l 
ODD HARMONICS ONLY | 

bo == 0 

Y, = b, + b, + b; + b; əəə ə € = 0 
Hence b, = — ( b; + b; + b; bs ə ) 

Yoo =G — و0‎ + A; — G, ..... 
Hence Qı = Y, 一 ( — Gs + G; — G; 000 ) 


In order to interpret the results it is eonvenient to 
have a chart or table in which the readings may be 
inserted, as otherwise the proper signs, additions and 
subtraetions may become confused. Charts of this 
sort are shown in Tables I, II and III for both odd 
and even harmonics. 


CONCLUSIONS 


While this analyzer requires intelligent operation 
and is not fully automatic, it offers the advantages of 
simplicity, ease of construction, and rapidity of opera- 
tion. It also frees the operator from the extremely 
wearisome amount of multiplication, addition and 
subtraction required by any extensive analysis with 
the usual schedule methods, and any one reading may 
be checked in a few moments without going through a 
series of calculations. With similar curves, as for in- 
stance the magnetizing current of a transformer at 
different densities, where the general character of the 
analysis is well known, but the actual magnitudes of the 
harmonics are desired, very great speed may be 
obtained, about five minutes being required for the 
determination of the third, fifth, seventh and ninth 
harmonics. Table I shows some comparative times 
required. The accuracy is about the same as that 
obtained by schedule analysis unless a large num- 
ber of ordinates is used. In a number of test anal- 
yses of curves of known composition, or by com- 
parison with other methods of analysis an average 
maximum error of about 2 per cent was found, the 
error being expressed in per cent of the amplitude of 
the fundamental. With some types of curves slightly 
larger errors may occur. If the notching bar is care- 
fully constructed the only chance for mechanical 
errors lies in the slipping of the clutch or dragging of 
the scale, and as the adjustment of these parts is 
simple and not at all critical this can not readily occur. 
Recently in the Research Division Laboratories analy- 
sis of all odd harmonics up to and including the 41st 
were made quite easily, the time required not being 
measured, but being far less than that necessary for 
any computational method, and the personal wear and 
tear being negligible instead of extreme. In general 
it might be stated that for certain classes of work this 
analyzer is extremely convenient and can be easily 
manufactured. It still is not the ideal analyzer but 
must be used with an understanding of its limitations, 
under which conditions it is very satisfactory. 
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TABLE I TABLE II—(Continued) 


TABULATION OF ANALYZER DATA Sine Coefficients 
For Odd and Even Harmonics and Constant Term 
Cosine Coefficients Net 
2-5 یی‎ Order of | Analyzer | Divide Then | Harmonics | Sum for 
wee Harmonic] Reading | Reading | Giving | Subtract|.——-—-———— 01 
Order of | Analyzer | Divide Then Harmonics Sum for by R 
Harmonic} Reading | Reading) Giving |Subtract pee bi ا ا راخب بھرہ ا سس‎ 
by + E ds 3 一 Ce +a -Q = 
05 5 — 6 + üş = 
HEBES “aaa و وا‎ uksa ص ا ر‎ 07 7 * -Q = 
bə. 4 bs + bio oa aş 9 * + ay = 
bs 6 bs + bs = an 11 * qə 
4 Š biz x dış 13 $ +ay = 
bs 10 k: +b = 816 15 * -a = 
bs 12 : 017 17 * "ar = 
b; 14 * - + b; = 
bs 16 * Y Sum (A) -....... 
bə 18 š 5o ے‎ *If known harmonics present greater than 12th subtract bnth + 
b 20 * — Din th etc. 
x Yer 0 b = — Sum B 
bu 22 5 + ۵ = ee ١ 5 
| bız 24 * SER al = Yo 一 Sum A 
۱ Note Yoo means value of ord. at 90 deg. 

Equation of curve:—Y = asin 6 + assin 3 9 + وه‎ sn560....... 
| Sum (B) -....... + bi cos 6 + bi 6093 0 + bs 60856. ..... 
| *If known harmonics present greater than 12th subtract bs: + باه برع‎ ۰ ۱ 
| TABLE III 

Sine Coefficients ۱ COMPARISON OF TIME REQUIRED FOR HARMONIC ANALYSIS 
| BY DIFFERENT METHODS 
Net 
Order of:| Analyzer | Divide Then Harmonics | Sum for No. 
Harmonici Reading | Reading | Giving | Subtract; ل‎ a Har- | Minutes 
by + 一 Method Time monics per Authority 
solved Coeff. 
a» 4 — aş +0 — —— —— -—— 0, = سس‎ 
aş 6 — aş — aş = Steinmetz......... 10 hr. 10 60.00 D. O. Miller 
04 8 | 一 a — Schedule.......... 3 hr. 8 22.5 D. O. Miller 
as 10 7 +a = Schedule.......... 1 hr. 8 7.5 F. W. Grover 
06 12 — Schedule.......... 2.5hr 17 10.6 Dellenbaugh 
a; 14 * 一 a; = Schedule.......... 15 min 3 5.0 D. O. Miller 
as 16 * — Coradi Mch........ 13 min 10 1.3 D. O. Miller 
ag 18 * + aş Coradi Mch........ 7 min. 5 1.4 D. C. Miller 
aio 20 ۱ * = Schedule.......... 30 min 6 5.0 Dellenbaugh 
011 22 - 011 Electric Mch....... 3.5 min 6 0.6 Dellenbaugh 
a12 24 a — Woodbury Mch. 5 min. 3 1.7 Dellenbaugh 
Sum (A) = ...... Woodbury Mch.....| 4.6 min. 3 1.5 Woodbury 
: Woodbury Mch.....| 9.75min 6 1.63 Woodbury 
“If known harmonics present greater than 12th subtract — ass + Schedule.......... 7 min 3 2.3 Woodbury 
asnth etc. | Schedule.......... 22 min. 6 3.7 Woodbury 
Ye RE b = Yo — Sum B 
Yo =... bi = 1/2 ( Yo — Yiso) — Sum B 
Yigg =... . . . 3 al = 1/2 (Yoo — Yeon) 一 Sum A 
Yen =. š 2 d 
Note: Yo. Yoo, Yıso, etc. means the values of ordinates at 0 deg., 90 deg., 
180 deg., etc. 
Equation of curve then:— ۱ 
Y = by + asin 6 + 62 819 2 6 + aşsin3 6 + طاو یه‎ 4 6 . 
+ bi cos 0 + وم‎ cos 2 0 + bscos3 0 + b4 cos 4 0. 
TABLE II : 


TABULATION OF ANALYZER DATA 
For Odd Harmonics Only 
Cosine Coefficients 


Net 
Order of | Analyzer | Divide Then Harmonics 
Harmonic | Reading | Reading | Giving | Subtract | تس‎ 
by + _ 

bs 8 bç + bis 

bs 5 bs . 

b; 7 7 

bo 9 * 

bu 11 * 

bis 13 * 

bis 15 * 

bu 17 * 


Sum (B) e....... 
*If known harmonics present greater than 12th subtract bt + وط‎ o. ete. 
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Radiation from Transmission Lines” 
.BY CHARLES MANNEBACK 


Electrical Engineer, Brussels, Belgium 


Review of the Subject: — The aim of this investigation was to 
contribute something to our knowledge of traveling electromagnetic 
waves. Among the transient phenomena that occur along trans- 
mission lines, these are still little known, especially as far as their 
“attenuation” and "“distortion,” or change in shape near the ۵ 
front” are concerned. : 

It is first seen that the “classical theory” of the propagation of 
electrical disturbances along lines, as it has been chiefly developed 
by Heaviside and Poincare, does not give a correct representation of 
the facts near the wave front, because tt assumes an instantaneous 
penetration of the current in the wires. It is shown that this theory 
is an “unidimensional” one, as it considers only one space variable, 
the coordinate along the line, and, from an electromagnetic point of 
view, amounts to identifying the traveling waves with plane wave 
phenomena. 

Steinmetz”s theory of the radiation from traveling waves (TRANS. 
A. I. E. E. Feb., 1919) is then examined, and, as Carson pointed out 
(Jour. A. I. E. E., Oct. 1921), found based on a misconception of 
the propagation of the electromagnetic field near the wires of a line. 
İl is remarked that this theory amounts to propagating longitudinal 
electric waves, a conception in conflict with the basis of Maxwell's 
theory. This latter proves very easily that, along a perfect line, 1. e., 
without ohmic or leakage losses, plane electromagnetic traveling 
waves are propagated without distortion and without attenuation; 
hence that there is no radiation. Losses do not change anything 
from the radiation point of view (Mie. Ann. Phys. 2, 1900). Inthe 
theory of Steinmetz, the radiation was the controlling factor at high 
frequencies. 

Although the question could be considered as settled, in so far as 
Steinmetz raised it, it is felt that the conditions under which traveling 
waves are started must be elucidated, in order to decide whether, 
during the transient process of establishing a plane electromagnetic 
traveling wave, any loss of energy can occur by radiation in free 
space. This would reduce the radiation to a transient phenomenon 
instead of to a steady one, as assumed by Steinmetz, i. e., to an effect 
produced when the condition of plane wave is departed from, as at 
the origin or the end of a line (tend effect”). 

The problem, considered in its broad aspect, 1. e., from the 
electromagnetic point of view of Maxwell, appears to be of great 
complexity. It involves a study of what might be called time-three- 
dimensional space transient, while the most complicated transients 
considered in electrical engineering are time-one-dimensional space 
transients, as was already noticed. It was found possible, however, 
to decide whether the radiation is a factor of enginering ۵ 
in allenuating and distorting waves. 

There, first, the distribution of current along the conductors of a 
line is shown to be a very close approximation to the actual, unknown, 
distribution, with regard to the possibly existing radiation. Then, 
from this assumption, the radiation of the system is calculated at a 
corresponding approximation, and its amount found to be negligible 
compared to the heat dissipated in the line during the same time. 


"This paper is the result of an investigation submitted in 
partial fulfillment of the requirements for the degree of Doctor 
of Science from the Massachusetts Institute of Technology, 
1922. It has been undertaken under the direction of Professor 
V. Bush, to whom the author wishes to express his fullest 
appreciation. 

Presented at the Midwinter روخ‎ of the A. I. E. E., 
New York, N. Y., February 14-17, 1928. 


From this it is concluded that the effect of the radiation upon the 
attenuation and the distortion of the waves along the line must also 
be negligible, compared to the effect of the joulian losses, and in this 
result a proof ‘‘a posteriori" of the correctness of the initial assump- 
tion is seen. The first idea of that procedure, but limited to station- 
ary waves, is probably to be credited to M. Abraham (Phys. Zeit. 2-p. 
329-1901) who applied it to the calculation of the radiation from a 
single isolated wire (oscillator), and found it in practical agreement 
with a more elaborate theory, based on Maxwell’s equations. As 
in this latter application the radiation is larger than the ohmic 
dissipation, the conclusions of this investigation are even strength- 
ened. ۱ 

An application to the steady radiation from a transmission line 
oscillating freely at one of its natural frequencies shows, even at very 
high frequencies, that the power radiated is negligible compared to 
the heat dissipated. For instance, a 100-kilometer transmission 
line, No. 00 B. & S. wire, oscillating at three million cycles per 
second, wastes by radiation only 1/8600 of what is wasted by heat. 
This is in accordance with the resulis given for steady a-c. traveling 
waves along a line by Carson, (loc. cit.). 

In the case of a traveling wave suddenly started at the origin of a 
perfect line, it ts found that, when the wave has become plane, the 
amount of energy 


3 
E d 1077 I? joules 


is wasted by radiation. d is the distance in.cm. between the two 
wires of the line and I the constant current suddenly flowing in the 
line, in amperes. This energy is carried to infinity by an electro- 
magnetic field limited to a thin spherical shell of variable depth, (not 
thicker tham d), expanding at the speed of light. <A reflection of the 
wave at a free end of a line is shown to add to the preceding an amount 
of radiated energy 
4 d 1075 I? joules 

and only half that amount for a grounded line. A complete trans- 
position causes a radialion four times as large as the sudden starting 
of the wave. Compared to the joulian dissipation during the time 
the wave takes to travel once along the line, the radiated energy is of 
the order of 1/10,000th. 

The influence of the radiation on the attenuation and distortion of 
traveling waves is thus found entirely negligible in engineering 
practise. Other possible factors are suggested. 
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INTRODUCTION ` 


NE of the chief reasons at present limiting the 
development and interconnection of high-power 
and long-distance electrical transmission lines 

is the trouble due to the occurrence of abnormal volt- 
ages and abnormal currents in the line and the con- 
nected apparatus. Whenever any change occurs in 
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the steady state conditions of an electrical system, 
e.g. by connecting or diseonneeting a load, or by 
producing or interrupting a short circuit, the new 
steady state conditions are generally reached after a 
very short time, during which so-called “transient 
phenomena”” or “transients” take place. As the flow 
of effective and reactive power are not the same in the 
different parts of the system during the first and the 
second steady state, the necessary adjustment is made 
by transient electromagnetic “traveling vvaves,” 7. e. 
waves of both voltage and current, carrying energy 
from one part of the system to another. 

The values of voltage and current in these waves, 
or as we shall say the amplitude of the waves, may 
reach abnormally high values and become destructive 
to the line. The steepness of the front of the wave 
also, z. e., the high rate of change of the voltage and 
current with the time may be harmful. Therefore, a 
knowledge of the “attenuation” or diminution of 
amplitude of the wave with the time, and of the “dis- 
tortion” or change in the wave shape, is of engineering 
Importance. Our present knowledge of the subjeet is 
still far from being complete. Joulian and dielectric 
losses produce both attenuation and distortion, but 
there may be other causes. C. P. Steinmetz! considered 
the possible influence of radiation, but his results may 
be questioned. 

From this began the present investigation, aiming 
towards the determination of the radiation effect upon 
the decay of traveling waves along transmission lines. 
It was found that the radiation effect is entirely neg- 
ligible compared with the effect of joulian dissipation, 
and hence that, from an engineering point of view 
radiation can be ignored in connection with trans- 
mission lines. 


CLASSICAL THEORY OF THE PROPAGATION OF 
WAVES ALONG ELECTRIC LINES 


The classical theory of the propagation of waves 
along a line, 2. e., a system of straight parallel outgoing 


and return conductors, rests upon two assumptions: 


the displacement currents in the dielectric, parallel to 
the conductors are negligible, and the so-called “four 
distributed constants” of the line can be determined. 
They are the resistance, inductance, leakance and 
capacity per unit length of the line (R, L, G, C). These 
quantities are functions of the frequency for sinusoidal 
steady alternating current; for transient currents, they 
are somewhat indeterminate, and an average value must 
be assumed. 

The fundamental equation for the propagation of 
the voltage V or the current J is under these assump- 
tions 


0? V 0? V 
LC Ə t2 oz 
(1) 


1. C. P. Steinmetz—Proc. A. I. E. E., Feb. 1919, p. 249. 
—Transient Electric Phenomena, 1920; Sec. III, 
Chap. 8, 9; See. V. 


Journal A. I. E. E. 


where t denotes the time and z the space coordinate 
along the line. From this equation, well known con- 
clusions can be drawn, of which the following only 
need be considered here: 

1. An electric disturbance, z. e., a discontinuity of 
voltage or current, is always propagated along any 
line at the constant speed v = 1: يد‎ L C, whether there 
is resistance and leakance or not. ` 

2. The disturbance is attenuated exponentially 
with the time according to 

exp [ — 1/2 (R/L + G/C) t]. | 
3. The distortion of the disturbance depends on 
the value of the term 
1/2 (R/L — G/C) . 
In the case of transmission lines, the attenuation i 15 the 
controlling factor, and the wave has died out before any 
appreciable distortion has taken place. Thus a steep 
wave front will very nearly keep its shape during the 
propagation. 

This conclusion, however, is based on the two assump- 
tions made; the first amounts to saying that the varia- 
tions of V and I along the line will not be too rapid, 
and the second, that the penetration of the current in 
the wires is instantaneous, neither of which is true at the 
very front of the wave. From Maxwell's point of view, 
the classical theory amounts to a plane wave theory, 
z. €., the electromagnetic field around the wires lies in 
planes; this implies that the dissipation is taking place 
everywhere in the plane wave, instead of being localized 
in the wires. 'This has been shown very clearly by 
O. Heaviside in his “Electrical Papers." 

We thus consider the question of change in the wave 
shape near the front of the wave, due to the resistance 
of the line, still as an open question. 


DR. STEINMETZ'S THEORY OF THE RADIATION 
Dr. Steinmetz's theory of the radiation from trans- 
mission lines deals only with the electromagnetic 
field. The electric and the magnetic fields are 
supposed to propagate from the wires perpendic- 
ularly into space, thus implying a flow of energy 
from the wires into space, and also non-transversal 
electromagnetic waves, two conceptions in conflict 
with Maxwell's theory. It is shown?, as an immediate 
deduction of Maxwell's equations that plane electro- 
magnetic transversal waves are propagated undistorted 
and unattenuated, at the speed of light in the sur- 
rounding medium, along a pair of infinitely long straight 
perfect conductors, 2. e., along a line without joulian 
or leakage losses; thus there is no radiation. Non- 
perfect conductors, as it has been shown, introduce more 
complexity, because the waves are no more plane waves, 
but nothing is changed from the radiation point of 
view. The flow of energy in the dielectric is evidently 
parallel to the wires, with à small component directed 
2. H. Hertz, Wiedemann's Annalen, 36, p. 1, 1889. 

Die Ausbreitung der elektrischen Kraft, 1892. J. R. 


Carson, Jour. A. I. E. E., Oct. 1921, p. 789. 
3. G. Mie, Annalen der Physik, 2, p. 201, 1900. 
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into the wires if there are losses. This settles entirely 
the question, in so far as Dr. Steinmetz raised it: 
traveling plane waves along a line are not attenuated 


nor distorted by radiation, z. e., by escape of energy in, 


a direction perpendicular to the line. Plane waves 
were an implicit assumption in the criticized theory; 
however, very near the source of a disturbance; or near 
the end of the line where reflection occurs, or in general 
near any place where there is a change in the geometrical 
or electrical behavior of the line, the electromagnetic 
field, of course, cannot be plane. 

This raises a new question: If plane traveling waves 
are only a limiting state, is there not any radiation 
taking place when the electromagnetic wave is initiated? 
Or, otherwise expressed, is radiation not an “end effect” 
with respect to the line? This involves the considera- 
tion of the whole surrounding space as the seat of the 
phenomena, and we need to start from a broader point 
of view than we did. 


POSITION OF THE PROBLEM FROM MAXWELL’S 
POINT OF VIEW—OUTLINE OF AN ENGI- 
NEERING SOLUTION 


According to a classification made by Steinmetz, 
we can say that the transients considered in electrical 
engineering are either simple time transients, as in 
machinery, or simple space transients, like voltage 
and current distribution along a transmission line under 
steady a-c. conditions, or both time and space transi- 
ents, as in the case of a line under transient conditions. 
We have already noticed that the notion of space 
involved in these transients 1s unidimensional or 
depending on one space variable only (equation 1); 
this supposes implicitly that the electromagnetic field 
of these transients is in planes. We are led, therefore, 
to consider more extended classes of transients; three 
dimensional space transients, as occuring near the ends 
of a steady a-c. line or near a wireless c. w. antenna 
and time and three dimensional transients, when steady 
conditions in such systems are changing. The only 
means of studying such transients is to resort to Max- 
well’s equations of the electromagnetic field, connecting 
the electric force E and magnetic force H. From these 
equations, follows the same equation for E and H, wz.: 


eu OE 4۰7 ۸*۲ OE 


c? 0 Ë c? Ot 
0? E 02 E 02 E 2 
Om? oy ' oz (2) 


E and H are space vectors, C represents the velocity 

of light in vacuo, e and u are the numerical values of 
the dielectric constant and magnetic permeability of 
"he medium and y the conductance of the medium 
ے‎ pressed in electrostatic units. In air, € and u are 
"wanity and y zero; in conductors, e is to be taken as 
ero, From (2) is shown that a disturbance, z. e. a local 
<Aisontinuity of E or H is propagated in the medium at 
+ he speed €: V eu. Equation (2), of which (1) is a 
E> articular case, is the “damped wave equation," well 


known in mathematical physics. Although much 
work has been done on this equation even in electro- 
magnetics, no solution has been obtained that seems 
directly available for the present problem, as far as we 
know. They generally deal with steady state sinusoi- 
dal waves along wires, without considering the source 
producing them, nor how they are initiated. (Som- 
merfeld, Ann. Phys. 67, 1899. Mie, Abraham, Ann. 
Phys. 2, 1900, Hondros. Ann. Phys. 1909, ete.) What 
we want is a knowledge of the whole electromagnetic 
field proceeding from a disturbance of voltage or current 
initially localized at a point of the line; the values of the 
field near the wires will give the amplitude and shape 
of the traveling wave, the values of the field at great 
distances from the line will give its radiation. 

The mathematical difficulties involved in a general 
solution are indeed very great, but we believe that the. 
position of the problem as it was sketched is the only 
one that can lead to a full answer to the question of the 
attenuation and distortion of traveling waves along 
wires. ۱ ' ۱ 

We found it, however, possible to decide whether the 
radiation is of engineering importance in the attenua- 
tion of the waves, either traveling or stationary. By 
radiation, we mean, as usual, the constant limit S 
towards which tends the flow of energy, as defined by 
Poynting’s theorem, through a sphere enveloping the 
whole electrical system, the transmission line in the 
present case, when the radius of the sphere is increasing 
above any fixed value. The time integral of the pre- 
ceding limit S will be called the “total radiation" 
and denoted by W. Wewill show first that the classical 
theory of waves (neglecting radiation) gives a “theoret- 
ical" distribution of current along the line that is a 
very close approximation to the “true” distribution, 
2. €., the-unknown distribution taking place when the 
radiation is considered. Knowing the current in the 
line, one can calculate the corresponding radiation that 
must exist; an approximate value will be obtaiñed by 
substituting the “theoretical” current for the “true” 
current. If, now, the amount of radiated energy is 
found to be quite negligible, say one to a hundred, 
compared to the joulian energy dissipated in the wires, 
we get an a posteriori proof of the rectitude of our 
premises, that is, that the true distribution of current, 
is given at, say, one per cent by the classical theory. 
As attenuation and distortion of the waves are produced 
by the dissipation of energy in the system, it is logical 
to conclude that the radiated energy, a cause hundred 
times less in magnitude than the energy dissipated by 
joulian heat, will produce a correspondingly small 
effect in attenuating and distorting the waves, 7. e., 
the radiation will be a negligible factor from an engineer- 
ing standpoint. 

The method is one of successive approximations, 
as used in experimental physics. It has already been 
used for calculating the steady radiation of given types 
of antennas, generally without being substantiated. 
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(Pierce, El. waves, Chapter IX. F. Cutting, Proc. 
I. R. E., p. 129, 1922). In the case of a straight oscil- 
lator, Abraham found the present method in quantita- 
tive agreement with the more elaborate solution based 
on the full solution of Maxwell’s equation. (Phys. 
Zeit, 2, p. 329, 1901). It is to be noticed that in the 
ease of antennas, the radiation is the controlling factor 
and the joulian dissipation only secondary, while in the 
ease of transmission lines the conditions are reversed. 
This strengthens our conclusions even more. 

We first consider perfect conductors, and later take 
into account the influence of joulian losses. 

When any kind of disturbance of current is prop- 
agated along a perfect thin wire, Heaviside shows how 
the corresponding electromagnetic field is set up by 
means of spherical traveling waves. (El. Theory, 


` Vol. I, art. 54, etc., where what he calls the “motion 


of a charge" is entirely equivalent to the propagation 
of a “current along a wire.") The reflection of current 
and the corresponding electromagnetic field generated 
by spherical waves originating at the point of reflection 
are indicated. (in El. Theory, Vol. II, art. 393, 394.) 
This theory neglects the diameter of the wires, com- 
pared to the length of the disturbance. At that order 
of approximation, the current (or traveling wave) is 
propagated undistorted and unattenuated at the speed 
of light in the air, and reflected as usual, even with 
regard to possible radiation. | 

We thus take the distribution of current along the 
line, either a stationary or traveling sinusoidal one in 
the steady a-c. case, or a traveling rectangular one in 
the case of a sudden impulse. 


STATIONARY RADIATION FROM A TRANSMISSION LINE 


We mean the radiation that takes place under steady 
a-c. conditions by means of stationary waves, such as 
the natural oscillations of a line, without joulian losses. 

Different modes exist, according to the terminal 
conditions to which the line is subjected. We consider, 
as an example, a line free at both ends, thus capable 
of half-wave length oscillations and multiples thereof. 

We assume spherieal polar coordinates, with the 
origin O, at the middle of the length / of the wire No. 
1, and the polar axis 4- 2 along the same wire. 

A length along the wires will be z, ro is the distance 
between a point P and the origin O,, 0, 15 the angle 
between O, P and the + z axis; r and 6 are the same 
quantities referred to any point of the wires instead of 
to Oi; © is the angle (longitude) that the plane P O, z 
makes with the plane of the wires, measured from the 
direction opposite to the second wire No. 2; d is the 
distance between the wires. (Fig. 1.) 

The distribution of the current along wire 1 is from 
the preceding: 


I = + ` cos ... . sin 7 e (3) 
with — 1/2 Z 2 Z + 1⁄2, and with reversed sign for the 


eurrent along wire 2. I 15 positive, if directed along 
+ 2, n is an odd integer giving the “order” of the os- 
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cillation; c represents the velocity of light; 7 is the 
maximum current. 

We have for the electromagnetic field due to the rate 
of change of a current element J d z, at distances large 
compared to the wave length of the oscillation, the 
values first given by Hertz under an equivalent form 
(Wired. Ann. 36, p. 1, 1889; or Ausbr. Elekt. Kraft., 
p. 153, 1892). ۱ ۱ 


E =H - 1/ )کتک‎ (4) 


where 7 d z is the same as the time derivative of the 
moment of Hertz”s doublet. 7 means I, where the 
time t has been replaced by £ — r/c; the current I and 


magnetic force H are expressed in electromagnetic `: 


units; the electric force E is expressed in electrostatic 
units. Ata point P, E lies in the plane of P and the 


Fic. 1 


wire, and is perpendicular to r; at the same point, 
H is perpendicular to this plane. The positive direc- 
tions are shown on the diagram. We consider the 
whole transmission line as a large sum of such current 
elements. At the point P, the elementary H” s add 
themselves, but all the E' s do not fall in the same direc- 
tion. They will do so, however, and will add, provided 
we recede far enough from the line, because all dif- 
ferences in direction then vanish, z. e. 0 can be taken as 
a constant. We can write neglecting quantities of the 
second order, 2. e., 1/1? as compared to 1/r: 


sin 0 d ə. 
n می پوس‎ 6 
where I, and I, the currents in the wires 1 and 2 are to 
be given their proper time phase £ — r/c. 
At great distances, we have with the same approxi- 
mation 
for wire 1 r = To — zcos 0, 
for wire2 r = ro — z cos 0, + d sin 00 cos do (6) 
Still with the same approximation, and neglecting 


E = H - 1 . 


Feb. 1923 


the higher powers of d/l, which is of the order of 10-3 
to 10-5, (3) gives 


1 TNZ 


. cos" (c 1 — ro 


d sin 0, cos ho (7)‏ 7 لك 
In this expression, z cannot be neglected as compared‏ 
with ro, because these quantities figure as argument of‏ 
the function cos. Substituting in (5), performing the‏ 
time: derivation and the integration with respect to‏ 
z, we get‏ 


+ z cos 69). 


mer ` 


E =H =(-1) * a1, ane 
sin 工 一 (ct 一 70) 
cos ( ə cos 00 ) cos o 
To 2 o) 8 


(8) 
where n is an odd positive integer, and the signs agree 
with the assumed convention. This formula repre- 
sents the steady electromagnetic field from an a-c. 


6۰ ExH +Z 


transmission line of length ,ا‎ at a distance 7, from the 
origin so large that the ration ۲/۲۵ is very small. The 
electric and magnetic forces E and H lie on pure spherical 
waves, traveling at the speed c from the origin on the 
hne. We consider the intersection of the + z axis with 
a sphere as a positive pole (origin of the angles 6) of 
the sphere. Then, the electric force is directed along 
meridians and the magnetic force along parallels; 
the sign in formula (8) is in accordance with the original 
conventions. (Fig. 2). 

The main difference between the “field of radiation” 
of a transmission line and that of a single wire line 
(antenna), is that, in the first case, the forces are only 
of the order of d/l, compared to the second case, which 
is very small. 


See note 1. 
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It is seen that the electric and magnetic forces are 
perpendicular to the direction of propagation, (the 
radius of the sphere), as it must be in the case of pure 
transversal waves; and thus Poynting’s vector E x H 
is perpendicular to the sphere. Hence, we have for 
the flow of energy or energy radiated per unit of time 


C 
S - سب‎ | 2 (9) 


integrated over the area of the sphere of radius r,, 
with d Z as element of area.5 

Substituting the values (8) for E and R, and perform- 
ing the mtegrations, and. taking the time average, 
we obtain: 


nd 


S=C. m2/2 (=~ jn i? ergs/sec. (10) 


Or 
2 . 2 ; 2 

= 15. m( ne ) 72 = 15 (+22 ( I? watts (11) 
where 7 is here the maximum current in amperes and 入 
the wave length. The power radiated by an oscillating 
transmission line, open at both ends, is thus, for a 
given current, numerically proportional to the square 
of the ratio of the distance between the wires to the wave 
length of the oscillation. 

We compare the power radiated to the power dis- 


` Sipated by heat, taking as an example the line whose 


constants are 
No. 00 B & 5 
d = 183 cm. 
l = 100 km. 

Total 0-6. resist. line = 48 ohms 

Total external inductance = 0.288 henry 
We have, for the power radiated S 
forn=1 A = 2 X 107 cm. 

S = 0.5 X 107” watts per max. ampere 
f = 1500 eycles/sec. 
n = 2001 A 2 104 6۰ 
S = 0.2 watts per max. ampere 
f 二 3 x 10“ cycles/sec. 

Due to the sinusoidal distribution of current, the 
power dissipated by heat is only 48/4 watts per max. 
ampere at low frequency; it is about 60 times more at 
three million cycles. ۱ 

Thus, even at three million cycles, the power radiated 
by the line is only 1/60 or 1.7 per cent of the power 
wasted in heat at low frequency, or 1/3600 of the total 
heat wasted at three million cycles. Hence, the radia- 
tion of oscillating transmission lines is entirely negligible 
in engineering practise. 

Dr. Steinmetz’s formulas amounted to writing the 
radiated power proportional to d . //A? instead of 
0/۸۵ 

Similar results are obtained in the case of forced 
oscillations of a line or impressed a-c. and load to 
Carson”s formula. (JOUR. A. I. E. E., Oct. 1921, p. 
789.) 

Instead of considering an ideal sinusoidal current 


5. See Note 2, 


kə ə Xor 


2 mem. 


= 
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distribution along the line, as when there are no losses, 
we might introduce the actual distribution given by 
the hyperbolic theory. There would be only more 
complication in the integrations, but the order of mag- 
nitude of the result would remain the same, 7. e. pro- 
portional to 02/22: thus quite negligible. 


TRANSIENT RADIATION FROM A TRANSMISSION LINE 
We mean the radiation that accompanies the produc- 
tion of traveling waves. When a disturbance of cur- 
rent is suddenly produced at the end O of a semi- 
infinite, thin, perfectly conducting wire, it is propagated 
in first and very close approximation as seen above, 
unattenuated and undistorted, at the speed of light, 
along the wire. The wire is supposed to be the only 
conductor in the empty space. The accompanying 
electromagnetic field is propagated, at sufficient dis- 
tance from the source, by spherical waves, also at the 
speed of light. Its expression at a point P is, when r 
is great, z. e.,such asto make 1/7” negligible compared to 
1/7: 
E = H = pg ا‎ ۷ 
sin 0 
r denotes the distance P O, 6 is the angle that P O 
makes with the wire; 7 the current at the origin O, 
given as a function of time £, vrhere £ has been replaced 


(12) 


m. 
و‎ - 


c is the speed of light. J and H are ex- 
pressed in electromagnetic units; E is in electrostatic 
units. .7 is positive when flowing along the positive 
direction of the wire. The electric force E is directed 
along meridians; the magnetic force H along parallels. 

This formula is substantiated* by noticing that it 
verifies both Maxwell’s equation for r large, and that 
for 0 very small, it represents plane electromagnetic 
waves traveling at the speed c along an infinitely thin 
wire; 


E = H = لت‎ (13) 


with l = + sin 0 being the distance of a near by point 
close to the wire. When the traveling disturbance 
of current reaches the end of the wire, it can be absorbed 
by the terminal apparatus without reflection; this is 
equivalent to starting at the end a new disturbance 
— Í, of opposite sign, but of same direction of propaga- 
tion as the first, so that each cancels the other along 
the wire supposed to be extended. Ií there is partial 
or total reflection, a new superimposed disturbance of 
current K Î is supposed to be sent along the wire 
from the end towards the origin, thus propagating in 
opposite direction, the value and sign of K, which is 
the coefficient of reflection, depend on the terminal 
conditions. 

We shall consider, just for simplicity, as only dis- 
turbance of current the sudden starting or stopping 
of a constant current at the origin, 2. e., “rectangular 


6. See Note 3. 
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waves" of current. The stopping will be considered 
as the starting of a current of opposite sign, in the same 
direction of propagation. 

The case of a transmission line reduces itself to the 
preceding case. We consider a line composed of two 
parallel semi-infinite wires originating.in O, and 0, 
and at the distance d apart. The coordinates used are 
the same as before, and are shown on the diagram, Fig. 
3. A disturbance is produced at the origin, consisting 
of a suddenly flowing constant current 4- I at O, and 
— I at .وم‎ At great distances from the source, the 
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electromagnetic field of the disturbance will be given 

by the sum of two expressions like (12) where due regard 

is given to the position of the sources O,and Oşand their 
T 

influence upon the “retardation” t 一 v and also to 


the sign of the disturbances. This supposes that the 


presence of one thin wire does not interfere with the 


electromagnetic field of the other wire at great dis- 
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tances; the linearity of the electromagnetic equations 
permits the addition of the fields. 

When a sufficient time has elapsed after the disturb- 
ance started at the origin, z. e., when its spherical front 
will have reached a great distance, in the sense we gave 
to that expression, the electromagnetic field of the 
disturbance can be considered as composed of two parts. 
(See Fig. 3.) One, in the space common to both “front 
spheres” O, and O,, centered upon the sources O, and 
Oz: the other in the space between the two spherical 
surfaces, which can be considered as a very thin spher- 
ical shell of variable depth 

d sin 0 cos ۵ 
this being equal to the distance O, O: projected in the 
direction considered (6, $). The electromagnetic field 
outside both spheres is obviously zero. In the shell, 
it is given by (12), with the proper sign for J, and lasts 
at a fixed point in the space during the time 


| d/c sin 0 cos © | 
the bars meaning an absolute value. Inside both 
spheres the electromagetic field is vanishing, 7. e., of 
the order of 1/r?, except near the wires where it is very 
nearly a plane wave field. It is given there by the 
difference of two expressions like (13) 
_2f _ 7 
01 02 

pı and pz being the distances of the point close to the 
wires. This is the well-known value of the electro- 
magnetic field of a plane wave traveling along a per- 
fect two-conductor line, as we considered it first from 
the classical point of view. Thus, besides the energy 
stored in the field of the traveling wave near the wires, 
there is also a small amount of energy present in the 
indefinitely expanding thin spherical shell of which we 
spoke, and which the classical theory does not account 
for. We consider this energy as the transient radia- 
tion from a transmission line, when acted upon by a 
sudden applied e. m. f. producing the current I in the 
line, because this energy is carried to infinity and thus 
is a net waste for the system. Its amount, calculated 
in a way entirely similar to that used in the case of the 
steady radiation above, 157 

W = 3/24 ۰ 10-° ۰ I? joules (15) 
where d is the distance between wires in cm. and I 
the sudden started constant current in amperes. 

We consider as an example a rectangular pulse of 
current, z. e., a traveling disturbance of constant current, 
extending over a length 入 cm。 The total waste by 
radiation will be twice as much as given by (15), be- 
cause the same transient radiation occurs when sud- 
denly stopping the current at the origin as when sud- 
denly starting it. The total energy in Joules carried 
by the pulse along the wires, when it is far enough from 


۹ Ls inf 
the source, is given by oT 


(14) 


where L is the self 


7. See note 4. 
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induction in henries per em. of the line, I the current 
in amperes. For the line already considered, L is 
equal to 23.8 X 10-9. Thus the ratio of the radiated 
energy to the energy contained in the plane electro- 
magnetic field of the pulse is given by 


1/۸ = . ۶۸ 


a 


or, for a pulse extending over a length À one hundred 
times the distance d between the wires, only 0.1 per cent 
of the total energy of the pulse. It 18 to be noticed that 
the fundamental frequency of the pulse in this example 
is as high as 820,000 cycles. Taking as effective 
resistance the resistance at one million cycles (17.3 
ohms per km. of line, Steinmetz loc. at p. 202), it is 
found that when the pulse travels only one third of its 
own length, it wastes in heat an amount of energy 
equal to the energy radiated. This shows how neg- 
ligible is the effect of radiation. | 

When a rectangular wave of current reaches the end 
of a single wire, its stopping generates a new spherical 
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wave, expanding at the velocity of light and centered 
upon the end R; (Fig. 4) and which is tangent to the 
first spherical wave originating from O. 

Similarly, in the case of a two wire line, a spherical 
shell of electric and magnetic field is generated, cen- 
tered upon the two ends R, and Rə, and tangent to the 
first shell originating from O; and 02. (Fig. 5. Thus 
the total radiation will be the energy carried into infinite 
space by the two expanding spherical shells. As the 
two shells have the part X, X, partly in common, 
where their fields add to each other, the total energy 
radiated will not be exactly that obtained by adding 
the energies of both shells separately, but will be a little 
less. It is shown that taking this sum an error is com- 
mitted that will not exceed the ratio of the angle 


بک 


to 2, +. e., of the order of d/l,? or between 10-? 


to 10-5 or less. Hence, to that approximation, we can 
say that the starting current amd successive reflections 
will each add a given amount to the energy lost by radia- 


8. See note 6. 
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tion. “The solution of specific cases is thus directly 
given by the preceding. We give the following results 
as illustrations of the method. 

1. The sudden establishment from one end of a line 
. of a sustained continuous current I, without reflection 
at the other end, wastes 


8 d 10-? I? joules (16) 
by radiation; this result is independent of the length 
of the line or of its material; it depends only on the 
distance d between the two wires of the line. In this 
formula and in all the following, d is in centimeters and 
I in amperes. 

Taking the same 100 km. transmission line as pre- 
viously, and assuming only the d-c. resistance even 
during the transient conditions, it is seen that the 
energy wasted by radiation is only 1/10,000th of the 
energy wasted in heat during the time in which the 
line is being charged, and the steady d-c. conditions are 
reached. For a 1-km. line, the radiated energy would 
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still be only one per cent of the energy wasted in heat 
during the transient conditions. The same remark as 
to the order of magnitude of the radiation compared to 
the joulian losses applies to all following examples. 

2. A total reflection of a “rectangular wave” of cur- 
rent I at an open end of a line radiates? 


4 d 10-9 72 joules; 

at a grounded end, 
2 d 10-9 J? joules. (18) 
In general, we have a partial reflection when the amount 
of current K I, (in magnitude and sign), is reflected 
and thus the rest of the wave is absorbed in the ter- 


minal apparatus. The radiated energy there is given 
by 


(17) 


3/2 (1 — 2/8 K + K?) d 10-9 I? joules (19) 
of which the three preceding formulas are particular 
eases. (K = 0, K = — 1, K = 1) 

9. A complete transposition, z. e. a sudden rotation 
9. See note 5. 
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of 180 deg. of the plane of the wires of the line, produces 
a radiation of 

6 d 10-9 J? joules, (20) 
or four times as much as given by formula (15), be- 
cause a transposition amounts to a sudden stopping 
of the current plus a sudden starting of the same cur- 
rent with opposite sign, in the same direction of prop- 
agation or to a sudden starting of twice the same cur- 
rent in the same direction, with opposite sign. 

Hitherto, we considered a current flowing suddenly 
in the line at maximum value or a “rectangular wave" 
as we called it. Any other law for the current as a 
function of the time at the origin, as for example a 
sinusoidal law, can be assumed, with only changes in 
the calculations, but yielding the same conclusions 
as to the magnitude of the radiation effect. 

We also neglected the influence of the ground and 
supposed a transmission line isolated in space. -A 
perfectly conducting ground can be accounted for by 
taking with the line its image with respect to the 
ground; the resulting radiation is obviously less. 

If, instead of considering disturbances between the 
two wires of a line, we consider disturbances between 
line and ground, the distance d between the wires would 
have to be replaced by twice the height of the line 
above the ground, in the ease of a good conducting 
ground. This can magnify the radiation effect about 
twenty times, but keeps it still negligible compared to 
the effect of heat dissipation in the conductors. 


CONCLUSIONS’ 


The present conclusions settle the question of the 
relative unimportance of radiation in the attenuation 
and distortion of traveling waves along lines. They, 
however, leave the way open for the investigation of 
other more important causes of attenuation and dis- 
tortion, among which we suggest the transient pene- 
tration of the current into the wire at the front of the 
disturbance, as mentioned above; the dissipation of 
energy in the ground due to the fact that some lines 
of force of a two-wire transmission line reach the ground 
and cause an attenuation as in the case of wireless 
waves; the non-homogeneity of the transmission line 
for high-frequency waves, 7. e., the slight periodic varia- 
tions that occur in the constants of the line due to 
presence of poles and curved spans. The influence of 
changes in direction in the line, as far as radiation is 
concerned, can be found by the method employed here 
and is negligible. 


SUMMARIZING 


1. A criticism of the classical theory of propagation 
of traveling waves along a two-wire transmission line 
has been given showing that, as far as the shape of the 
wave near the front is concerned, this theory is still 
unsatisfactory. 

2. Dr. Steinmetz’s theory of the effect of radiation 
on traveling waves has been examined. 


Feb. 1923 


3. The position of the problem of radiation of lines 
has been sketched from Maxwell’s electromagnetic point 
of view, in connection with the accepted engineering 
theories of transients. 


4. Anengineering solution of the problem is offered, 


showing that under steady a-c. as well as under tran- ` 


sient conditions, the effect of radiation upon electromag- 
netic waves propagating along transmission lines is quite 
negligible as compared to the 7 of heat dissipation in 
the wires. 


The procedure used here in the case of steady a-c. 
conditions is to be considered an extension of a method 
due to Abraham?" it is different in the case of transients, 
and thought to be new. 


5. Other possible causes of attenuation and dis- 
tortion of traveling waves are suggested. 


NOTE 1 


The integration to be performed is obtained by sub- 
stituting the time derivative of f,— f, given by (7) 
into (5); so 


rn sin 0 mnd 
E = H مح‎ 2001010. SIN OG 
l To / 
+ 2 
TNZ . T" 
. COS م‎ cos 7 , sin a (c t — To 


— 2 


+ zcos00).d. 

But 

sn 一 六 (ct — To + z 6085 0) = sin 一 (et 

— ro). cos = = (z cos 00) + cos == (ct 
5 ro) . sin ا‎ (z cos 00) 


The second term contributes nothing to the integration, 
. wn : š : 
for sin E (z eos 0,) changes its sign with z We 


thus have only to perform the integration 


an (|^ TNZ Tn 
E m cos 一 7 一 .os 一 (z cos 00( dz 
— 1/2 
۲ LM 
= f cosx . 608 0 2 . 7 2 
° n r 
-72 
n r 
+ 
- ورد‎ cos (1 +a)z . dz 
n r 
2 


T 


dra 
+1⁄2 f 


T 


cos (l—a)x.dz 


10. Phys. Zeitschrift, 2, p. 239, 1901. 
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_ _ Sin (1 + a) 2 + sin (1 — a) 2 
1 + ۵ 1 - 0 
= E Z ہہ‎ eos aş” (n odd) 
۱ Th 
Lau 00 ` (— 1)”. . COS ( — cos 6, ) 


which, multiplied by the constant factors that we did 
not write through, gives formula (8). 
NOTE 2. 
The flow of energy S, given by (9), is 


8 | كر - 5 


integrated over the area of the sphere of radius ro; 
the values of E and H are to be taken from formula 
(8), and d > is the element of area in spherical polar 
coordinate, 7. e., 

To” sin 0, d 0, d Po 
We thus have to integrate 


sin? Z 


一 (ct 一 To) 


e 


در 


2rnd 9 
4200 


T 2T 


f f cos? ( Ty cos 0) cos? 0, sin 0, ۲۶۰ 00. ۰ 


0 =0 وب‎ =0 
But we have, dropping unnecessary indices, 


ni co ۵ . 8 و‎ =r 


0 
f cos? ( T 
0 
The double integral thus reduces to n rg” — the 
value of S becomes 
S = cË ) T nd ” (ct — ro) 


of which the time — 15 də. by replacing sin? by 
1/2. Hence formula (10) 


NOTE 3. 


The equations (12) for E and H can be substantiated 
by proving that they satisfy both Maxwell’s equations 
for r large, 7.e. when 1/72 is negligible compared to 


and 


cos 0 )sin 0d0 =1 


y sin? ۰ 


1/r. Maxvvell”s equations are, in air 
° E 

1/6 T MEE curl H 
oH 

— 1/e s curl E 


where E and H are space-vectors, E being expressed 
in electrostatic units and H in electromagnetic units. 

The eomponent of the curl of a vector, at a given 
point, in a given direction is defined as the line integral 
of that vector taken around a unit area located at the 
given point, of which the positive normal coincides 
with the given direction. 


104 . MANNEBACK: RADIATION FROM TRANSMISSION LINES 


The curl eomponents of E and H in the direetion of 
the radius (see Fig. 3) are seen to be zero or of the second 
order, (z. e., of the order of 1/77), thus negligible. The 
components of curl E along a meridian line, and of 
curl H along a parallel line are obviously zero, because 
the vector is then perpendicular to the line of integra- 
tion. 

The component of curl H in the positive direction of 
a meridian line, (7. e., towards the increasing 0), is 
given by 


curl H = 
—H(r+Ar).(r+Ar)sn0.A0-+H(r).r.sn 0. A0 


Ar.rsin 0 ۵ ۵ 


where H (r) means the scalar value of H at the distance 
r. But we have, neglecting quantitles of the second 
order, and using the value of H given by (1): 


of 
a 1 + cos 0 
sin 0 


_ dE 
= H (r) 一 l/c—T . Ar 


H (r + Ar) = H (r) + . AT 


Substituting in curl, H, at the same approximation, 
we find the first equation of Maxvvell, 

The same procedure with regard to the signs leads 
to the second equation of Maxwell. 


NOTE 4. 


Formula (15) represents the energy carried to infinity 
by the spherical shell of variable depth 


d sin 0 cos $ 
in which the value of the electromagnetic field is given, 
for — 7/2 Z و‎ Z + 7/2, by formula (12) 


0 ل 1 


E=H=1/2. uno 


where 7 is the constant value of the suddenly impressed 


= Ə T : 
current. For z 7 2/¢ i , the sign of both E and 


H is to be changed. 
At a given point on the sphere of (large) radius r, 
this electromagnetic field lasts during the time 
| d/c sin 0 cos ¢ | (c speed of light) 
the bars meaning a quantity that is positive. 


The total energy radiated W will be given by the 
time integral extended over the sphere of radius r 


. of Poyting's vector, which is directed everywhere along 


the radius outwards the sphere; 7. e. 


c ; 
Tg Ë x H | d/c sin 0 cos ¢ 
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integrated over the sphere of which the element of ara f: 

15 | ١ 
72 sın 0 d 0 0 ç 

We thus have to perform 


Di [ee s J. ۱‏ تی 


| 8/6 sin 0 cos ó | ۲۶ sin 0 d 0 d ó 


But we have 


27 


f | 605 ۵ | dQ = 4 
0 


and 


f ]1 + cos 6546 = f (1 + 2 cos 6 + cos? 0) d 0 
“0 0 

| = +0 + ۵ - ST 
The value of W is thus | 


W -3/2dT elmag C. G. S. units 


. which is equivalent to (15). 


NOTE 5. 


The total reflection “with change of sign” of a 
rectangular wave of current at the free end of a line 
can be considered as the result of the stopping of the 
incoming wave + J at the end, by starting at the end 
a new wave — I along the same direction of propagation 
as that of the incoming wave, plus starting a new 
wave— I in the opposite direction, 7. e., back into the 
line. This last wave is the wave "reflected with change 
of sign." "The sign of a wave bears no relation at all 
to the direction of propagation. 


The reflection at a free end of a line is thus equiva- 
lent to starting simultaneously two waves — I from the 
end in two opposite directions. The electromagnetic 
field of such a double disturbance along a single wire is 
easily deduced from that given in formula (12), for a 
single disturbance. It is, with regard to the signs 


(=D 14 cos 0 7 (- Í) 1- ۹ء‎ 


r sin 0 r sin 0‏ مش سیف 
(2D 1‏ | | 
m r sin 0‏ 


(— T) is the current started in both directions at the 
free end, equal in absolute value to the incoming cur- ا‎ 
rent, and where the time t, of which J is a function in : 
the general case, “has been replaced by t — r/c, r being 
taken from the free end. 

The case of a two-wire (transmission) line reduces 
to the preceding one as was seen. Thus, in the case 
of a rectangular wave of current along a line, the elec- 
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tromagnetic field due to the sudden starting at the end 
Of the line of a current (— I) in both directions is: 


TEE‏ وت وو و ےت 


This field exists in a spherical shell of variable depth, 


d sin 0 cos $, expanding at the speed of light, as it has 
been shown. It lasts at a fixed point in space during 
the time 


| d/e sin 0 cos ¢ | 


the bars meaning a positive quantity. The amount of 
radiated energy W can be found exactly as in the pre- 
ceding note. VVe have: 


T 


.. me f ane | d/c sin 0 cos ف‎ | 
000 


6 
Ar 


what reduces immediately to: 
W=4df 


which is equivalent to formula (17). 

In general, when the amount of current K I (in sign) 
is reflected into the line, the rest being absorbed in the 
terminal apparatus, we have for the electromagnetic 
field, with regard to the sign, 


elmag c. g. s. units 


E ے‎ H = = Í ۱ 1 + cos 0 
r sin 0 
(Kİ) 1- cos6 
T a s sin 0 ۱ 


which exists in a spherical shell of variable depth. 
d sin 0 cos $ 


expanding at the speed of light. 
The radiated energy is, as above, 


W = ic “idi 
f di )1 — K) + 4 + K) cos 0 
[ sin 0 ۳ 
9-0 %=0 


| 0/6 sin ۵ cos ¢ | r? sin 0d0do 
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The value for W follows 


W = 8/2 (1 — 2/8 K + K?) d 7۶ elmag. c. g. s. units 
which is equivalent to (19). 


` NOTE 6. 


The position of the point X,, limiting the length 
X, X; (Fig. 5) along which the two spherical shells 
can have parts in common is determined by the inter- 
section of a circle of radius r centered upon 01 and of 
a circle of radius r — l centered upon روط‎ r being very 
large compared to l. This allows to approximate the 
circles by parabolas of thesame curvature, z. e., of which 
the parameter p is equal to the radius of the circles; 
the axis of the first parabola will coincide with the wire 
O, R,, the axis of the second with the wire O, Re We 
so have | 


y—d-42(r—Dz 
1 < V2rx Par. 1 


with obvious axis of coordinates. Calling (ğı ə) the 
coordinates of the point X;, we deduce 


له - 72۳-56 تسس TEN‏ 


As ع‎ is very small, و‎ is the distance of X, to wire 1— 
We have thus, from the preceding equation, neglecting 
the square and higher powers of ۰: 


n/r = 1 . 


which proves that, viewed from a point of the line, the 
length X, Xə, does not cover an angle larger than 


Par. 2 


P 
—-ru-d) 


ag , provided £ is small compared to r. 

Hence, calculating the energy radiated by the two 
spherical shells in contact as if they were not in contact, 
i. e., integrating over ع‎ instead of m minus the angle 
n/r, will introduce only an error of the order of d/l, 
thus quite negligible. 
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Dissymmetrical Electrical Conducting Networks 
۱ BY A. E. KENNELLY 


Fellow, A.I. E. E., 
Professor of Electrical Engineering, Harvard University and Mass. Inst. Technology 


T is well known that a symmetrical T-connection 
| of three simple impedances a o, o b, and o q, Fig. 1, 

in which the two arms q o, and b o are electrically 
equal, and with their common terminal o connected to 
the ground or return conductor gh through the im- 
pedance o q, can be completely replaced,* at any single 
assigned alternating-current frequency, by a certain 
symmetrical r, a’ b” g” h”, Fig. 2, in which the two 
pillars a’g’ and 6’h’, have equal impedances; also 
that either of these two symmetrical systems can be 
completely replaced, at the same frequency, f by a single 
smooth line conductor a" b", Fig. 3, having uniformly 
distributed series impedance in its conductor and also 
uniformly distributed lateral leak admittance between 
the conductor and the return conductor g” h”. Such 
a line conductor may be regarded as subtending or 


2 o ¿ 


a" : Ú“ 


n 2 

" | 

168, 1,2 AND 3 — Sy MMETRICAL T, T AND SMOOTH LINE, Equiva- 

LENT AND ÍNTERCHANGEABLE 

277 at the frequeney in question, a certain 
ni ex hyperbolic angle 0, and also a certain single 
FE Impedance or “characteristic impedance," zo 
— 8 Consequently, both the symmetrical T con- 
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ducting system of Fig. 1, and the symmetrical z of 
Fig. 2, may be regarded as possessing the same complex 
angle 0 and surge impedance zə, by virtue of equiva- 
lence. In any network of conductors, carrying alter- 
nating currents of the assigned frequency in the steady 
state, any of the three equivalent systems of Figs. 1, 


“ 6, 0, 0 3 [^ 
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Fics. 4, 5 and 6 一 COMPOSITE LINE OF THREE SECTIONS WITH 
ITS EQUIVALENT DISSYMMETRICAL T' AND T 


2 and 3 may be interchanged, without disturbing the 
distributions of potential, current and power, at and 
outside the terminals. of the system. 

Moreover, if two or more simple smooth line con- 
ductors, of the type indicated in Fig. 3, are connected 
in series to form a composite line as in Fig. 4, it is known 
that the composite line system a” b” g” b” (Fig. 4) 
can be replaced by its equivalent T of Fig. 5,1 or by 
its equivalent + of Fig. 6, and that these equivalent. 
conductors are, in general, dissymmetrical. Fora given 
composite-line system of Fig. 4 operated in the steady 
state at a single frequency, there can be only one equiv- 
alent T and also only one equivalent 7. On the other 
hand, neither the equivalent T of Fig. 5, nor the equiva- 

* Bibliography (1). 
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lent r of Fig. 6 can determine the elements of smooth 
line section in the equivalent composite-line system 
of Fig. 4. In the absence of specific information con- 
cerning the number and characters of the component- 
line sections, an infinite number of different composite- 
line systems like that of Fig. 4, could be found as 
equivalents for the dissymmetrical T or r of Figs. 5 
and 6. The question thus arises as to what is the sim- 
plest smooth-line system which can replace a given 


a 


R m 


Fig. 7—NETWORK OF IMPEDANCES REPLACEABLE BY One 
DissYMMETRICAL T OR 7 WITH RESPECT TO THE Two PAIRS OF 
TERMINALS, G, g, AND b, h. 


dissymmetrical T' and its corresponding dissymmetrical 
m? 
Again, if any network of conductors, such as that of 


- Fig. 7, carrying alternating currents of a single fre- 


queney, and two pairs of terminals are selected arbi- 
trarily from the system, such as a, g, and b, h; then the 
systema behaves with respect to these pairs of terminals 
like a certain equivalent T of Fig. 5, and also like the 
corresponding equivalent 7 of Fig. 6, these two systems 
being mutually equivalent and, in general, dissym- 
metrical. The question then also arises as to what is 
the simplest smooth-line conductor system which can 
be regarded as replacing the network, with respect to 
the pairs of terminals a, g and b, h? 

A symmetrical network of conductors, with respect 
to two pairs of terminals, may be defined as one whose 
equivalent 7 and 7 are symmetrical, like Figs. 1 and 2. 
A dissymmetrical network is then one which, with 
respect to the two pairs of terminals, has a dissym- 
metrical equivalent T and r, like Figs. 5 and 6. In 
general, networks are dissymmetrical, and symmetrical 
networks are particular cases. 

Although we always discuss a network with reference 
to two pairs of terminals, such as a, g and b, h in Fig. 
7; yet the number of terminals involved may be re- 
duced, in particular cases, from four to three, by the 
merging of two terminals. Thus, we might select the 
two pairs of terminals, a, g and b, g, whieh employ the 
terminal g in common. 

As has already been pointed out, a network which is 
symmetrical with respect to two pairs of terminals 
becomes defined either by a certain symmetrical T 
and r, or by a certain hyperbolic angle 0 and ac- 
companying surge impedance zo. We can then readily 
compute what changes in alternating-current potential 


difference, eurrent and power will be brought about at 
one of the pairs of terminals by any assigned electric 
change steadily impressed on the other pair. When, 
however, the network is dissymmetrical, we shall find 
that it requires for its definition one additional charac- 
teristic. A convenient additional characteristic is the 
"inequality ratio" q; so that the network acquires 
three characteristics 6, ومع‎ and q, which are capable of 
being measured in any given case; by following a suit- 
able technique at each pair of terminals. Having 
ascertained these three determining characteristics 
for the network, we can either draw immediately its 
equivalent T and 7; or we can compute what change 
in potential difference, current and power will be pro- 
duced at one pair of terminals, say o, g, by a given 
change established steadily at the other pair b, h. 
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Fras. 8, 9 AND 10—DissyMMETRICAL T REDUCED TO SYM- 
METRY BY THE RETENTION OF A PERMANENT SERIES LOAD AT THE 
b TERMINAL WITH THE CORRESPONDINGLY LOADED 7F AND 
SMOOTH LINE. 


The computation is effected by following the same pro- 
cess as in the symmetrical case; but with a slight modi- 
fication. 

The subject can best be studied with reference to 
continuous-current lines and networks; 7. e., with 
reference to the frequency of zero. The results so 
obtained are immediately applicable to alternating- 


` millimho. 
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eurrent systems, by the substitution of complex num- 
bers for real numbers, in the well known way. 

Fig. 8 shows a simple dissymmetrical T. In this 
particular case the branch a o has 1000 ohms and the 
branch b o 2000 ohms. Two methods suggest them- 
selves for dealing with it; namely, (1) by reducing it 
to a terminally loaded symmetrical T, and (2) by re- 
taining the dissymmetry, but making measurements 
from each end in turn. We may consider these two 
methods suecessively. 

(1) Method of Reduction to Symmetry. The simplest 
way to reduce the T a, o, b of Fig. 8 to symmetry, is 
perhaps to eut off the portion B, b from the preponderat- 
ing arm, so as to leave the remainder o, B equal to the 
branch a, o. We thus obtain symmetrieal T a, o, B 
with a load B, b, in this case 1000 ohms, permanently 
attached to the terminal B. The corresponding loaded 
symmetrical r is shown in Fig. 9; where the same load 
B',b', is permanently attached to the B' terminal. 
Fig. 10 gives the terminally loaded smooth line a" B", 
corresponding to the combinations of Figs. 8 and 9. 
It has an angle of 0.622362 hyp., a surge resistance of 
Zo = 3316.63 ohms, a total conductor resistance R of 
2064.14 ohms and a total leakance G of 0.187649 
The same load of B” b” = 1000 ohms is 
permanently attached to the 5” terminal. 

It is readily possible to substitute for the dissym- 
metrical T a,o,b, of Fig. 8, the terminally loaded 
smooth line of Fig. 10, in any -steady-state system. 
The smooth line is readily dealt with in computation, 
and the permanent load B” 5” merely modifies the 
conditions existing at B, in a manner deducible directly 
from Ohm’s law. If preferred, the load may be dealt 
with by assigning to it an auxiliary angle 6”, so that 
the position angle at a, becomes 0 + 6”. 

Although the reduction to symmetry effected by 
reducing the higher-impedance branch of the T to 
equality with the lower-impedance branch, is perhaps 
the simplest, especially in a continuous-current case 
like that of Fig. 8; yet it is by no means the only method 
of procedure open to use. Thus, the dissymmetrical 
T of Fig. 8, might be replaced by its corresponding 
dissymmetrical equivalent r, as in Figs. 11 and 12. 
The v, a’ b” g” h' of Fig. 12 may be reduced to sym- 
metry, by cutting off part of the leak a’ g’ as a terminal 
load, and leaving the remainder at a” G, Fig. 13, equal 
to the leak b” h”. This leaves the symmetrical 7, a” 
b” Gh", of angle @ = 0.622362 and surge impedance 
Zo = 5125.67 ohms, but with the leak a” g” of 0.058824 
millimho permanently attached to a”. This symmet- 
rical x might now be replaced either by a 1+0 08 
symmetrical 7, or by a smooth line with the same ter- 
minal leak load. 


Again it would be possible to reduce the 7 of Figs. 
$ and 11 to symmetry, by employing permanent pairs 
of unequal terminal series loads, as for instance 100 
ohms at a, and 1100 ohms at b, so as to produce a 
symmetrical T with 900 ohms in each branch. This 


could be done in an infinite number of ways, each giving 
a different 6, a different zə, and a different pair of per- 
manent terminal loads. It would also be possible 
to reduce the + of Fig. 12 to symmetry, by employing 
pairs of unequal terminal leak loads. There would be 
an infinite number of such possible pairs, with cor- 
responding values of 0 and z) Consequently, a dis- 
symmetrical network admits of being reduced to sym- 
metry for computation, in one way with a single ter- 
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CORRESPONDINGLY DissyMMETRICAL 7, WHICH 1s REDUCED 
TO SYMMETRY BY THE RETENTION OF A PERMANENT LEAK AT 
THE a” 'TERMINAL. 


DISSYMMETRICAL 7 REPLACED BY İTS 


minal series load, in one way with a single terminal 
leak load, in an infinite number of ways with a pair of 
opposite terminal series loads, and also in an infinite 
number of ways with a pair of opposite terminal leak 
loads. In practise, the two single-load methods are 
preferable and easter to use. 

(2) The Method of Applying the Inequality Ratio ۰ 
In this method, the system is left in its dissymmetrical 
state; but measurements are made from each pair of 
terminals in turn, to determine the apparent surge 
impedance from each end. 

Returning to the dissymmetrical T of Figs. 8 and 11, 
we can measure the resistance Ras offered between the 
terminal a and the return conductor or ground g, when 
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the system is freed or opened at the opposite terminals 
b, h. We may also measure the resistance Ra, offered 
between the same terminals at a, when the end b is 
grounded, or the terminals b, 7 are short-circuited. 
Following the usual rule for finding the angle of a 
symmetrical smooth line, the apparent angle of the 
system is expressed by: 

tanh 0 = V و‎ numeric Z . (1) 
In the particular case of Fig. 8, the resistance Rag 
would be 2428.57 ohms, and Ras = 6000 ohms, from 
which tanh 6 = 0.636209 and 0 = 0.751779 hyp. ra- 
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Aum $123.11” 
Fic. 14—Network or Fic. 11 WITH SINGLE ANGLE 6, BUT 
WITH SURGE IMPEDANCES Zoa AND Zob AS DETERMINED FROM 
EACH ۰ ۱ 


dians. Similarly, according to the regular rule for a 
symmetrical line, the apparent surge impedance Zoa 
from the end a, is | 

= v Ry . Ra, ohms Z (2) 
which in the case considered, is 38817.26 ohms. If now 
we repeat the measurements from the.b end, succes- 
sively freeing and grounding the a terminal, we find 
Ras and R,,, which would be 7000 and 2838.33 ohms 
respectively. As ۱5 before, 


tanh 0 = V R,,/Ry = 0. 636209 


whence 0 = 0.751779. This is the same value as was 
found by (1), from the measurements at a. We may 
express this relation as a general rule thus: If in any 
dissymmetrical network of conductors, such as that of 
Fig. 7, we take two pairs of terminals like a, g, and b, h, 
and measure, at an assigned frequency, the impedance 
offered by the network between each pair, with the 
opposite pair first opened and then shorted, we may 
denote these impedances by Ras, R,, and Ras, E, 
ohms Z respectively. Then the ratios R,,/R,; and 
R.,/Rs; will be equal to each other, and to the square 
of the tangent of the angle subtended by the network 
with respect to those four terminals. 

On the other hand, if the network is dissymmetrical, 
the apparent surge impedance from b, or 

Zo = V Roy š Roo ohms Z (4) 
will not be the same as Zo. İn the case considered, 
Zo = 4453.46 ohms. The geometrical mean Za» of the 
two opposite apparent surge impedances is defined by 
Zab = V Zoa . Zob ohms Z (5) 

is a characteristic property of the network, in this case 
4123.11 ohms. It may be called the “geomean surge 
impedance” of the network. 

The ratio of the a-end surge impedance Zoa to the 
geomean surge impedance z¿ may be called q, the 
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“inequality ratio” of the system. It will be defined by: 
q = 2a/%an = V Rag/Rog = v Raj/ Roy = V 2 0/20 ۱ 


۳ s e m ا‎ 
p: + @ g. + سر‎ 

In this case the inequality ratio is q = 0.92582. 

We may thus represent either the dissymmetrical 
T of Fig. 11, or its eorresponding z of Fig. 12, by the 
dissymmetrical line system of Fig. 14. Here the line 
a, b, is to be regarded as a composite line, offering the 
same angle from each end, but different surge im- 
pedances. 

Terminal Potentials and Currents in a Dissymmetrical 
System. If we have any physically consistent distribu- 
tions of potential, current and power at the a and b 


numeric Z (6) 


terminals of a known dissymmetrical system, and we 


know the values of potential and current at one end, 
we can readily ascertain the corresponding values at the 
other, for one and the same frequency. 

In the case of Fig. 15, we have the same dissym- 
metrical system as in Figs. 11, 12 and 14, for which 6, 
Zo, and zo are given. A potential of 1.0 volt and an 
entering current of I, = 0.347826 milliampere, are 
found at the generating end a. Required the cor- 
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Fras. 15 AND 16—NETWORK LOADED WITH A TERMINAL İM- 
PEDANCE 0 AT HACH END ALTERNATELY 


responding values at the motor end b. The value of 
the load o at b need not be known. VVe have then: 


E, = Û (FE, cosh 6 — I, za sinh 0) volts Z (7) 
q 


and I, = q (1, cosh 0 — E, yo, sinh 6) amperes Z (8) 
where yo, = 1/20a, is the surge admittance, as measured 
froma. If و‎ = 1, so that the dissymmetry is made to 
disappear, the above formulas are those which apply 
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to any ordinary smooth line of 6 and zo. İn the case 
considered, E, = 0.21739 volt, and I, = 0.21789 mil- 
liampere. 
. If on the contrary, we know E, and I, only, without 
necessarily knowing the physieal conditions at the 
motor terminal b, we may proceed to find the cor- 
responding values at the generator end a. We obtain: 

E, = q (E, cosh 0 + و‎ zə sinh 0) volts Z (9) 
and T, = 5 (I, cosh 6 + E, yo, sinh 0) 

amperes Z (10) 

where Yon = 1/20: is the surge admittance as measured 
at b. If we make q = 1, these formulas become iden- 
tical with those pertaining to a simple smooth line and 
symmetrical system. In any such smooth-line case, 
formulas like (7), (8), (9) and (10) can be employed 
for any point along the line; but in the case of a dis- 
symmetrical system, they apply only to the terminals 
aand b. 

Again, if b is the generator end, and a the motor end 
of the dissymmetrical system, (see Fig. 16), and we 
know the potential E, and current I., at a, we can find 
the corresponding values at 6 as follows: 


E, = = (E, cosh 6 + I, Zoa sinh 0) volts Z (11) 
I, = q (Ia cosh 0 + Ea Yoa sinh 0) amperes Z (12) 
Similarly, if we know the generator-end values, E, and 
L, the corresponding motor-end values E, and I, 
become: 

E, = q (E, cosh 0 — I, zo» sinh 0) 


volts Z (13) 


J < r (I, cosh 0 — EP yo» sinh 0) amperes Z (14) 
Here again, if q = 1, the formulas are those which apply 
at the terminals of a smooth line and symmetrical 
system. Hence we may formulate the following 
rule for any dissymmetrical case. Having given the 
potential and current at one end, write the correspond- 


ing values for the other, with regard to the direction: 


of power transmission, as though the line were smooth 
and the system symmetrical. Use the surge impedance 
or admittance belonging to the end at which the 
values are known. Then an unknown potential E, 
at a, will be q = V Zoa/Zo times that found by the 
regular formula, and likewise an unknown current I, 
.at b. On the other hand, an unknown current I, 
at a, or a potential E, at b, will be 1/q times that found 
by the regular smooth-line formula. The power 
products E, Ta and E, L, will therefore not differ from 
those given by the corresponding smooth-line formulas. 


Position Angles at the Terminals of a Dissymmetrical 


System. If weload the motor end b of a dissymmetrical 
system, like that of Fig. 15, with an impedance c 
ohms Z, we may find the position angle of that end in 
the same manner as though the load were applied to 
the motor end of a simple smooth line; z. e. 


tanh ó, = 00 numeric Z (15) 


Thus in Fig. 15, the position angle ó, at b is 0.228437 
hyp. The position angle 6, at the generator end a 
is then obtained by the usual rule: 

Ôa = 0 + ó, hyps Z (16) 
which in this case is 0.980217 hyp. The potentials and 
currents at the ends of the dissymmetrical system are 
then related as follow: 

With E, and I, known, 


h ó, 

Ea = q ( E, = 5, ) volts Z (17) 
1 cosh 6, 

p rit I ) amperes Z (18) 


The impedance beyond b is c. The impedance at a is: 
Za = Zo, tanh ó, ohms ۶ (19) 
The power at a is the local product of E, and J, in 
vector watts, according to the regular rule. 
With E, and J, known, 


_ 1 sinh ó, : 

E, -— (E. E) volts Z (20) 
5 cosh ó, ; 

L, = q (1 p^r ) | amperes Z (21) 


Formulas (17) to (21) apply to ordinary smooth- 
line symmetrical systems, if q = 1. Whereas, however, 
in a smooth-line system, a position angle can be as- 
signed to each and every point of the line, we are only 
justified in assigning position angles to two points— 
the a and b ends—of a dissymmetrical system. 

Similarly, if 6 is the generator end (Fig. 16) and a 
the motor end of a dissymmetrical system, to which a 
load c is applied, the position angle at a is defined by: 

tanh 6. = 00 numeric Z (22) 
In the case considered, ó, = 0.268221 hyp. The posi- 
tion angle at b is: 
و۵‎ = 0 + 6, hyps Z (23) 
The potentials and eurrents at the two ends are thus 
related: 
With E, and I, known, 


sinh ó, 
E, = 0 ( Es- inh, ) volts Z (24) 
1 cosh ó, 
T; sS zn r1 amperes Z (25) 
With E, and I, known, . 
_ 1 sinh ó, 
E, = = (2. i. ) volts Z (26) 
5 cosh by ۱ 
I, = q ( I, ER) amperes Z (27) 
The impedance at the generator end is: 
Z, = zo tanh 6, ohms Z (28) 


Formulas (24) to (27) are identical with (17), (18) and 
(20), (21). 

We may formulate the deductions from the last 13 
formulas as follow: Im any dissymmetrical system for 
which 6, Zoa and zo are given, and which is loaded at the 
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motor end with an impedanee o, find the position angle 
at the motor end, with reference to that end”s surge 
impedance. If the potential and current are given at 
one end of the system, the values at the other end are 
found by using the ordinary position-angle formulas 
for a smooth line. An unknown E,, or L,, will then be 


q = V Zo. /Zo times the value so obtained. On the 
other hand, an unknown E, or Ia will be 1/q times the 
value obtained by the ordinary formula. This rule 
agrees with that given above for use with formulas 
(7) to (14). 

In view of the above mentioned reciprocal factors 
for potential and current at each terminal, leaving the 
power products unchanged, it becomes possible to 
reduce a dissymmetrical system to a, corresponding 
symmetrical system plus a terminal ideal transformer, 
without losses, and with a transformation ratio of q. 
It is doubtful, however, there would be any advantage 
in this plan over the method above outlined. 

In the case of an actual composite line, like that of 
Fig. 4, with definite sections, having respective angles 
6, Û», 03 . . . and surge impedances Zoi, 2o2, Zos, . . ., 
it is advantageous to work up the position angles at 
the successive terminals commencing with the motor 
end. The terminal position angle at a would thus 
become ö,, with a terminal surge impedance zo; while 
working from the other end, the terminal position angle 
would be ö,, with a terminal surge impedance of ۰ 
The impedance R., as measured at a would be . 

Ra, = 2o, tanh ó, ohms Z (29) 
while that measured at b would be 

Ro, = Zon tanh ó, ohms Z (30) 
These values are not inconsistent with (19) and (28), 
which pertain to a general network; whereas (29) and 
(30) pertain to a particular composite line with detailed 
sections. 


Recewing-End Impedance of a Dissymmetrical Net- 
work. If we define this quantity as the ratio of the 
e.m.f. impressed at the generator terminals to the 
current received through a load connecting the motor 
terminals, we find that for a load c at b, (Fig. 15), the 
receiving end impedance is: 

Zw = q (ze sinh 0 + c cosh 0) ohms Z (31) 
In the case represented in Fig. 15, Zi, = 4600.0ohms. 
When q = 1, this expression reduces to the ordinary 
value و2‎ sinh 6 + o cosh 0, for a simple smooth line 
or symmetrical system. 

When c = 0, or the b terminals are short-circuited, 
(31) becomes: 

Zw = q zo» Sinh 0 = Za sinh 6 ohms Z (32) 
which is the arehitrave impedance in Fig. 12, or 3400.0 
ohms. 

Similarly, for a load of c ohms Z at the motor end 
a, the receiving-end impedance becomes (Fig. 16.) 


Lu = - (Zoa sinh 0 + c cosh 6) ohms Z (33) 


İn the case of Fig. 16, this becomes 4800.0 ohms. 


When q = 1, the expression (33) again coincides with 
the ordinary smooth-line formula. Moreover, vvhen 
6 = 0, 


Zu = ” (Zo, sinh 0) = 2,, sinh 0 ohms Z (34) 


or 3400.0 ohms, the architrave impedance of the 
equivalent dissymmetrieal 7. 

We may formulate these results as follows. The 
receiving-end impedance of a dissymmetrical network 
loaded at the receiving end, is found by the usual rule 
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Fic. 17 —DISSYMMETDİCAL T OF IMPEDANCE REPRESENTING A 
NETWORK OF ALTERNATING-CURRENT CONDUCTORS 


for a terminally loaded smooth line, using the surge 
impedance of the motor end. The result must then 
be multiplied by q, when the generator is at a, and by 
1/q when the generator is at b. In case the load is 
short-circuited, the receiving-end impedance reduces to 
p” = Za sinh 6, the architrave impedance of the equiv- 
alent dissymmetrical 7. 

Technique of Measurements. We have assumed that 
the measurements made on the network are the four 
terminal impedances Ras, Rags Ros, and Reo. In cases 
where there is but little difference between Ras and Rag, 
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Fig. 18 —DıssYMMETRICAL EQUIVALENT T OF SAME 


or between R and R,,, the value of 0 may not admit 
of being determined with satisfaetory precision; al- 
though the value of Zab will usually be acceptable. In 
such a case, if the two ends of the network can be 
brought to the same testing table, it may be preferable 
to adopt the following technique: 

1. Measure, with b open, the impedance R. ohms 
at a, and the entering current J,ata. At the same time 
measure the vector p. d. e, volts thereby produced at 
the b terminals. From these data, رم‎ and ® can be 
deduced. 
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“2, Measure Ra ohms > at a. With this value 
and Ras, find Zoa, by (2). 

3. Repeat the measurements (1) and (2) reciprocally 
with the ends of the network reversed, so as to obtain 
Ry, ea and Ra. These will give وم‎ and zo, together 
with check values of 6 and ۰ 

Example of a Dissymmetrical Alternating-Current 


System. Fig. 17 offers a particular dissymmetrical 


R, z ۵395-07۲ 21°48'S "= اس‎ 2000”. Reg = 6946-22 30°1S'23'= 6000-j3500" 
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Fic. 19—NErvoRK or Fic. 17 WITH SINGLE ANGLE 6, BUT 


WITH SURGE IMPEDANCES Zoa AND Zob AS DETERMINED FROM 
12 ۸ 0۳ EN». 


T system, whieh might be the equivalent of a certain 
network between the two pairs of terminals a, g and 
b, h. The impedance a o, is 1000 + 7 1000 ohms, such 
as might be produced by an inductor of 0.1 henry and 
1000 ohms, at an angular velocity of w = 10000 
radians per second, (1591.5 ~). The impedance 
ob is 2000 — 7 500 ohms, such as might be produced 
at the same frequency by a condenser of 0.2 micro- 


farad, in series with 2000 ohms. The impedance 
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Fig. 20—SAME WITH TERMINAL LOAD AT b 


og is 4000 — / 3000 ohms. The equivalent m of this 
system is given in Fig. 18, and the equivalent composite 
line network in Fig. 19. The angle subtended by the 
network is 0 = 0.705768 + j 0.305124 hyps; or quad- 
ranting the imaginary part 0 = 0.705768 + j 0.194248 
hyps. The functions of this angle are: 

sinh 0 = 0.822652 Z 27° 237 05” 

cosh 6 = 1.228218 Z 10° 50’ 28” 

tanh 0 = 0.672581 Z 16° 32’ 37” 


Many of the formulas in Appendix A can be checked 
from Figs. (17), (18) and (19) by an inspection of the 
slopes of the various complex quantities presented. 

In Fig. 20, the system of Figs. (17), (18) and (19) 
is represented as being loaded at the b end with c 
= 8516.08 + 7 14409 .4 ohms, and supplied at the 
a end with an e. m. f. of 1.0 Z 0° volt at the frequency 
of reference. The terminal position angles, potentials 
and currents are marked on the Figure, in accordance 
with formulas (15) to (21). 

In Appendix A, a number of formulas have been 
collected relating to dissymmetrical networks and 
their equivalent circuits. The notation is in accord- 
ance with Figs. 21 and 22, and the numerical values 
pertain to the particular case of Figs. 14, 15 and 16. 
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Fics. 21 AND 22-——DISSYMMETRICAL T AND 7 WITH THEIR 
ELEMENTS SYMBOIIZED, ALSO NUMERICAL VALUES CONFORMING 
WITH THE PARTICULAR CASE or Fics. 11, 12 AND 14. 


In Appendix B, a corresponding series of formulas 
has been collected, relating to symmetrical networks 
and their equivalent circuits; 7. e., for the case q = 1. 
The notation is in accordance with Figs. 28, 24 and 
25. The numerical values pertain to the system having 
the same 6 as in Figs. 14, 15 and 16 and with its surge 
impedance 2, equal to the geomean surge impedance 
Zap Of the dissymmetrical system. The values of 
R, and R, in this case are respectively the geomean 
values of Rar, Ro, and Rag, Rog. 

We have considered, in all of the cases represented 
by the figures, that in any dissymmetrical T, all of the 
series impedance is lodged in the conductor a, o, b, 
and that the return conductor g,h, has negligible 
impedance. According to the known theory of T 
conductors, however, any desired share of series im- 
pedance can be transferred from a, o, b, to g, h, without 
disturbing the e. m. fs., currents and powers at and 
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beyond the system terminals, provided that the total 


series impedance on the a side of the leak o q is constant, 


and likewise the total series impedance on the b side. 
In other words a T-system may be converted into a 
H-system. 

Similarly, although in all the z-systems indicated 
in the Figures, all the series impedance is lodged in the 
conductor a’, 5”, yet any desired share of this imped- 
ance may be transferred to g’, b”, without disturbing 
the e. m. fs., currents and powers at and beyond the 
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Figs. 23, 24 AND 25—SYMMETRICAL T, T AND SMOOTH LINE 
HAVING THE SAME ANGLE FEEDER AS IN Fics. 14, 21 AND 22, 
ALSO WITH Zo = Zab OF THE DISSYMMETRICAL CASE. 


system terminals. In other words, the dissymmetrical 
m s may be converted into corresponding O-systems. 


SUMMARY OF CONCLUSIONS 


1. Any dissymmetrical conducting network, each 
element of which obeys Ohm’s law generalized for 
alternating currents, may be reduced to symmetry, 
with respect to two pairs of terminals, and at an 
assigned frequency, either by finding its equivalent 
T and equating the two arms, through the medium of 
a permanent series terminal load; or, by finding its 
equivalent ع‎ and equating the two pillars, through the 
medium of a permanent terminal leak load. The 
system then offers a certain angle 0, a surge impedance 


` 2, and a single terminal load, which may be treated 


as having an auxiliary angle 6”. 

2. Any such dissymmetrical network may also be 
regarded as possessing and being defined by (a) an 
angle 0, (b) a geomean surge impedance Za», and (c) an 
inequality ratio g. The electric behavior of the net- 
work, at the assigned frequency, with respect to the 
two pairs of terminals, may then be computed by the 
usual formulas applying to a simple smooth line; but 


with the further application of the coefficient q. 

3. The theory of simple alternating-current lines 
with smoothly distributed constants may thus be 
regarded as being a particular case (q = 1) of the 
general theory of . dissymmetrical networks having 
any degree of complexity. 


LIST OF SYMBOLS EMPLOYED 


ba, Oo; . position angles at a and b ends of a diss. 
network (hyps Z). 

61, On position angles at ends of a detailed 
composite line (hyps Z). 

Ea E, e. m. fs. applied to a and b ends of a diss. 
network (volts Z). 

£a, Cb, p. d. s at open a and b ends of diss. net- 
work by application of a testing 
e. m. f. at the opposite end (volts 2). 

f frequeney applied to netvvork (eycles per 

` second). 

G total conductance of a smooth sym. line 
(mhos 2). 

Gas Go, admittance of a sym. system, when 
respectively freed and grounded at the 
other end (mhos <). 

Gar, Gags admittance of a diss. system, at the © 
end, when respectively freed and 
grounded at the other end (mhos 2). 

Göy, Gig; admittance of a diss. system, at theb 
end, when respectively freed and 
grounded at the other end (mhos 2). 

g admittance of leak in equiv. 7 of a diss. 

۱ system (mhos Z). 
gu 02) admittances of leaks at a and b ends in 
equiv. r of a diss. system (mhos <). 
g' admittance of leak in equiv. T' of a sym. 
system (mhos Z ). 
g" admittance of each leak in equiv. r of a 
sym. system (mhos Z). 
0 angle subtended by a line or system 
(hyps Z). 2 
0! angle subtended by a terminal load 
(hyps Z ). ۱ 
0i, ,5و6‎ ... 0, angles subtended by successive sections 
| 01 a composite line (hyps 2). 

I, I current strengths at ends of a network 
_ (amperes Z) 

j=v-1 | 

y ۳۳ series adm. of a sym. or diss. z (mhos 2). 


V /رووج‎ Zoo inequality ratio of a diss. network 


(numeric Z). 


R 
Ry, R; 


Fan Ras 


Ry, Rog, 


Al و02‎ 


Yo 


Job,‏ رہہ لا 
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total conductor resistance of a smooth 
line (ohms Z ). 

impedance at one end of a sym. system, 
when the other end is freed and 
grounded respectively (ohms Z ).. 

impedance at the a end of a diss. system 
when the other end is freed and 
grounded respectively (ohms 2). 

impedance at the b end of a diss. system, 
when the other end is respectively 
freed and grounded (ohms 2.). 

impedance of leak in equiv. T' of a diss. 
system (ohms <). 

impedances of a and b leaks in the equiv. 
T of a diss. system (ohms <). 

impedance of leak in equiv. 7 of a sym. 
system (ohms 2 ). 

impedance of each leak in equiv. r of a 
sym. system (ohms 2). 

impedance of architrave in equiv. 7 
of a diss. system (ohms 2). 

impedance of a and b arms in equiv. T' of 
a diss. system (ohms Z ). 

impedance of each arm in equiv. T of a 
sym. system (ohms Z). | 

impedance of architrave in equiv. 7 
of a sym. system (ohms Z ). 

impedance of a load at motor end of a 
system (ohms Z ). 

surge admittance of a sym. system 

(mhos 2). ۱ 

surge admittance at a and b ends of a 
diss. system (mhos Z ). 


Yas = 252001 Yor geometrical mean surge admittance of a 


۱ Rə, /Re = Bö / Roy, = G. /G., = Gyu/Gy = 


pi/p° 


1/q 


1 
9 (p. + pə) 


R 


diss. system mhos Z ). 


generator-end impedances of a, diss. 


Zo Z by 
system (ohms 2 ). 

Z tas Z u reeeiving-end impedances of a diss. sys- . 
tem, including load (ohms Z). 

surge impedance of a sym. system‏ و20 
(ohms 2).‏ 

Zoas Zob surge impedance at a and b ends of a 


diss. system (ohms 2). 

Zab = V Zoa Zob geomean surge impedance of a diss. 
system (ohms <). . 

291, 202, . . . Zon Surge impedances of successive sections 
in composite line (ohms 2). 


angular velocity of frequency impressed‏ رج 2 ح ن 
on a system (rad. per sec.)‏ 

~ sign for “cycles per second.” 

Z sign for a complex numerical quantity. 

hyp. abbreviation for “hyperbolic radian.” 
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APPENDIX A 


List of formulas applicable to a dissymmetrical 7, ع‎ or Network. The numerical values pertain to the case 
ef Figs. 17, 18, 21 and 22. 


V 0 
R 


R (bi + p2) + pip. _ 
(R + pi) (R + pə) 


v (gı +92) +9192 _ ; 
(v + gı) (? + gə) tonne 


= 0.404762 numeric Z (35) 
= 92/9; = 66 _ 909-1 = 0.5 numeric Z (36) 
q 09b 0—1 
—--—— —— EAS — — +R ga + v 
nu M Ra / Ë, 一 V Ra / Ror = V Gog/ Gag V 6/7.۸ V 0/0۵ - أو‎ a 1 m 
= Za/Zay = 0.92582 numeric Z (37) 
a. 6 +? 
= M Ey / Rag = V Ry/ Rar = V Gag/ Gog 到 M Gas/Goy V 200/200 ” — E P x 
= 2o/2.5”— 1.08012 numeric Z (38) 
i (gy +92) 
© 1 2 
= eee = ( 45,13. ) cosh 0-1 = 0.3 numeric Z (39) 
V 2 
مم‎ = Yah - 5 (q cosh ۵ — 1) ) cosh 0—1) = 0.28284 numeric Z (40) 


10 


۰م 


gi/ p» 


“0a 
Zob 


sinh 0 


cosh 6 


tanh 0 


01 
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= pə/gı = وم = 6 رم‎ Ri = pR 


= p/g us G/v = Y ۷ ظط‎ = Zoa Zoo = Zat = Ras E 


01 9° 
= Rag Ry = 0: وم‎ +R (pi + pz) = v/ Ras Ras Roy Rog = V رم‎ p» Ki Re 
= V | رم‎ (pə + GO + وم | وم‎ (ps + R) + م‎ © ] = 17,000,000 = (4123.106) 
۱ ohms° ۶ (41) 
= g:/pi = vg = g/p = v/@ = N gm = Yoa Yor = Yar? = Gas Gog = Gag Gos | 


= gi g: + v (gı + g:) = VGas Ga Gu Gro = V {91 (gs + V) + gs v] İg: (gi + v) + ,و‎ | 
= 0.058824 x 10-6 = (0.2425355 x 10-3)? . mhos?Z (4) 


= y Rar Rag = X pi (Pı + 8) don azdı. = Q xa = 02 رھ‎ = 3811.26 ohms Z (43) 


- کب‎ Ras — Rag - te 


+R‏ وم 


falsis p» (p: + R) + رم‎ ate = üə = FI — 4453.46 ohms Z (44) 


pi + & 


V G(pm + p ۲ pip. _ _V )9۰ + 9) +9192 
q 


V 
= /p/R= V م‎ 9 = 24/0 = Yar/v = 0.82462 numeric Z (45) 
| Re, _ x (& + p1)(& + pə) ےم‎ 
E EN وس‎ RR - EM ə 
— Vv + nv + g: po tg). = -V Ga Gro Ges Cw _ 1.296148 numeric Z (46) 


= V R.,/R.; = V R, /R.; alban 1 


— — ——— —ÓÀ— —MÀ — 


6 (bi + p) + pip: . 
R? + 6 (gi + وم رم + (وم‎ 


_ |  »(n39) + 9۱ دق‎ Zə, Yas | 
د سد مو جال‎ C EET MET G.. Ga = 0.636209 numeric Z (47) 
= ورم‎ = gs/v = g. p = PİR: = و‎ cosh 0 — 1 = 0.2 numeric Z (48) 
= ووم‎ = gi/v = و‎ P = ۵/۱ = X cosh 0—1 = 0 4 numeric Z (49) 
= Za Í ت83‎ ) = 1000 ohms > (50) 
| + cosh ۵-1 ۱ ۱ 
= مك )مو‎ ) = 0 ۱ ohms ۶ (51) 
= z,,cosech 0 = V Rs; (Ra — — Ru) = کیہ‎ Ra (Ro, — Ry) = 5000 ohms ۶ (52) 
= Ya sinh 0 = 0.2 x 1073 mhos ۶ (53) 


Gas — Gas 


= Yab cosech 0 = V Goo (Gag 一 Gu) = V Gag (Gog =R Gry) ss 


Zap Sinh 0 = N Ru __ = 


um Fo, —— = Ry J 
“6 Goja = Rao V Ra و سی ےد‎ = Ew BL P "E. 
‘Rar Ro, Rag Ras 


5 = — .ا‎ = 8400 ohms ۶ (54) 


Gas Goo _ _ Gag Gy 


g 
= 0.294118 x 10-3 mhos > (55) 


— cosh 0 — 1 
= Wad جس سب اس‎ — 0.117647 2 10-3 mhos 4 (96) 


sinh 6 
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Yar (seem 2-1 ) = 0.058824 x 10 


sinh 0 


Zo, tanh 0 = q zə tanh 0 = q? Zo tanh 0 = 2428.57 
zə tanh 0 = rZ tanh 0 = a zə tanh 0 = 2833.3 
Zo, coth 0 = q Za» coth $ = 0? Zo coth 0 = 6000 


= Zə coth 0 = = Za coth 0 = x Zo, Coth 0 = 7000 
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mhos Z 


ohms Z 
ohms Z 
ohms Z 


ohms Z 


11 


(97) 


(58) 
(59) 
(60) 
(61) 


List of formulas applicable to a symmetrical T, v, or network, corresponding to the formulas of Appendix A. 
The numerical values pertain to the case of Figs. 23, 24 and 25. 


R,/Ry 
p'/&' 
p" /g' 
g'/ "م‎ 


sinh 0 


cosh 0 


tanh 0 


_ p! (p' + 2 @” B g" (g" + 2 v) _ 2 
G,/G, = — TL wy - ay 7 gx» tanh? 0 = 0.404762 
= g"/v = p'g' = p" g" = cosh 0 — 1 = versh 0 = 


mho? Z (42a) 
. IRA R” ” ” 2 
- zr _ ve 2 _ vg کا‎ 二 
۱ = 0.82462 numeric Z (45a) 
- Ax R R -— = = S EY 0۲/ج ۔‎ = G,/v= IFIP” = 1.296148 numeric Z (46a) 
c= 
= VR R; = P تست اود = م‎ = 20/10۲, = 17/0, = 0.636209 
| ; ۱ numeric Z (47a) 
= a (OES) = zy tanh (6/2) = 1480. 744 ohms Z (80a) 
"32 
= zü cosech 0 = رھ‎ sech 0 = V Rr رق‎ R.) = a = مه‎ = — LI = 5000 
ohms Z (52a, 
= yosinh 6 = G; cosh 0 = == = وا‎ x5») = yè/v = کته‎ 
| = 0.2 X 109 mhos Z (53a) 
= zə sinh 6 = R,cosh:0 = R, Vee e (e FER) = 20/6۲ = 3400 ۱ 
ohms Z (54a) 
= ون‎ cosech 0 = G, sech 0 = V G, (G, — Gy) = 5 x yg = Ee 
= 0.294118 × 10-3 mhos > (55a) 
= yo (2>1) = yo tanh (6/2) = 0.0871024 x 10“ ohms Z (56a) 
= z tanh 0 = BU IM = 2628 57 ohms > (58a) 
= z,coth 0 = p' + &' = 6480.741 ohms Z (60a) 
= go coth 0 = g” + v = 0.38122 x 10° mhos ۶ (62) 
= gə tanh 0 = TQ = 0.154149 x 10^ mhos Z (63) 


yg’ = 7/6۲ = "و‎ (g" + 2v) 


numeric Z 


2 sinh? (0/2) = 0.296148 numeric ۸ 
"م‎ QR! = R'/v = ام ام‎ +28") = R; R, = ۶و2‎ = (4123.106)? = 17,000,000 ohms? Z 
= G; G, = yo? = (0.2425355 x 10-39)? 


= 0.058824 x 10-6 


(35a) 


(39a) 
(41a) 
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Transient Conditions in Electric Machinery 
BY W. V. LYON 


Associate, A. I. E. E. 
Assoc. Prof. of Electrical Machinery, Mass. Inst. of Tech., Cambridge, Mass. 
€ 


Review of the Subject: — The vector method 4s just as useful 
in solving problems involving transient conditions in electric 
circuits as it has proved to be when the currents and potentials are 
steady sinusoids. As far as the writer is aware, the vector method 
for determining transients in rotating electric machines was first 
used by L. Dreyfus. Previously the method had been applied to 
fixed combinations of resistances, inductances and capacitances 
by Kennelly and others. 

By making certain assumptions that are, however, quite reasonable 
in many cases, the transient currents in nearly all of the common 
types of electric machinery are damped sinusoids. Fortunately the 
damping is exponential and is thus readily accounted for. It is 
interesting to trace the development of the method. In the solution 
of all problems in direct currents the potentials, currents, and 
circuit constants are real numbers. In the corresponding problem 
in which the applied potentials are steady sinusoids, these quantities 
are all represented by complex numbers. In.all other respects the 
working out of the solution is identical with that followed in the 
direct-current case. When the currents are damped sinusoids, they 
and the potentials and the circuit constants can still be represented 
There ts this difference, however; the vectors 
which represent the currents and potentials shrink exponentially as 
they rotate and the values of the circuit constants depend not only 
upon the frequency of the current, but also upon its rate of shrinking. 
Again the solution of any problem follows the same procedure that 


HE problem of calculating the transient currents in 
alternating-current apparatus has rightly been 
considered a difficult one. Fortunately, the 

difficulty lies not so much in the problem itself as in the 
inadequacy of the mathematical treatment that has been 
employed. If a vector method of analysis is used, very 
great difficulties are removed. This was first appre- 
ciated by the author early in 1921 while studying the 


‘paper “Short-eireuit Current of Induction Motors and 


Generators” by Doherty and Williamson. Later it 
was found that the vector method had already been 
developed by L. Dreyfus. Notwithstanding this 
previous publication, the subject was considered of 
sufficient importance to warrant the presentation of 
this paper.. 

As is the case with many other complex problems, the 
method of attack is presented most clearly when certain 
Ideal conditions are assumed. These conditions are 
that the resistances and inductances are constant and, 
In the case of rotating niachines, that the speed is con- 
stant and that the space distribution of air-gap flux is 
sinusoidal. After one is familiar with the method in 
its simplest form, it may be possible to introduce re- 

ements which will account for iron losses, variable 


permeability and harmonics. 


The form of the transient current in any system of 
mutually Inductive circuitsis determined by the differen- 
dal equations of the fall in electric pressure through 


: See Bibliography. | 
7 resented at the Midwinter Convention of the A. I. E. E., 
w York, N. Y., February 14-17, 1928. 


of the vector. 


would the corresponding one in which the currents are steady 
sinusoids. In both the steady and damped sinusoidal cases the 
circuit. constants depend upon the angular velocity of the vectors 
which represent the currents. Im the former, the angular velocity is 
purely imaginary while ün the latter it is complex, the real part being 
the rate at which the current vector shrinks and the imaginary 
portion, tts angular velocity. In electric machinery tn which 
rotating magnetic fields are produced, these fields shrink exponen- 
tially as they rotate when the currents are damped sinusoida. If 
these rotating magnetic fields are represented by vectors, the vectors 
will have a complez angular velocity just as do the currents. The 
e. m. f. which is produced by a steady sinusoidal variation of flux 
lags the flux by 90 degrees, whereas if the flux variation is a damped 
sinusoid, the angle of lag is less than 90 degrees, depending 01 
the damping. The mathematical relation, however, is the same, viż., 
the e. m. f. 28 proportional to the negative of the product of the flux 
and its angular velocity. It ts then readily appreciated that the 
form of the solution for the transient state is the same as that which 
18 used for the steady state. Before the method can be expected to” 
give as accurate results as are obtained when predicting the steady 
operation, considerable experimental data must be obtained in order 
to determine the best methods of measuring the necessary constants, 
for these may be somewhat different during the transient period than 


during steady operation. 
* * * * * * * * 


them. With the foregoing assumptions in regard to 
the ideal character -of the electric circuits, the differen- 
tial equations which apply to most of the eommon 
types of direct and alternating-current machinery may . 
be reduced to the form which is described as “linear 
with constant coefficients.” This being true, the 
transient currents consist of one or more terms of the 
general form: 
i = I e“ sin (wt + 6) 

The way in which this eurrent varies depends upon 
the values of œ and w. There are four types of varta- 
tions, only two of which occur as transients. If o and 
w are both zero, it is a steady direct. current of value, 
I sin 0. If ais a positive quantity, as it always is, and 
o) is zero, it is a diminishing direct current of value, 
I €” sin 0. If a is zero and o is a positive quantity, 
it is a steady sinusoidal eurrent of value, I sin 
(wt + 0). And finally, if neither o nor w are zero, it 
is a diminishing sinusoidal current. In each of these 
cases the current, whether steady or transient, may be 
represented by a vector which in general has a length of 
I €” and which makes an angle of (wt + 0) with the 
horizontal axis. The projection of this vector on the 
vertical axis is the instantaneous value of the current. 
The four cases are illustrated in Fig. 1, a, 6, c and d. 

In each of the diagrams I marks the initial position 
In the first case (a) the vector is fixed 
both in magnitude and in angular position. In the 
second ease (b) the magnitude of the vector diminishes 
exponentially although its angular position does not 
change. In the third case (c) the vector has a constant 
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magnitude but rotates with an angular velocity of w. 
In the fourth and general case (d) the magnitude of the 
vector diminishes exponentially while it rotates at a 
constant angular velocity of w. Thus the general form 
of transient current may be represented by a vector 
which rotates at a constant angular velocity and whose 
terminus travels along a logarithmic spiral. 

The velocity of the terminus of the vector has two 
components, one along the vector and one perpendi- 
cular toit. It is convenient to measure these component 
velocities in multiples of the instantaneous length of the 
vector. The term angular velocity is thus appropriate. 
This is a broader meaning of angular velocity than is 
ordinarily assigned to it, and for this reason, it is spoken 
of as the generalized angular velocity. Differentiation 
shows that the generalized angular velocity of this 
vector is, m=(—a+jw). See Fig. 2. Notice 
that the usual significance is attached to the symbol j. 
It indicates that the vector 7 w, lags the vector, — a, 
by 90 degrees. 


Projection 


7 
5° 
x 


I sin اه‎ 5 


Projection 


Zero Axis 
c d 


Fic. 1-127 Enp or EAcg VECTOR WHICH IS Drawn FROM 
THE ORIGIN 18 INDICATED BY A LARGE PERIOD; THE END OF 
SMALL VECTORS DRAWN FROM THESE PERIODS IS INDICATED 
BY ARROWHEADS. 


The form of the current variation is determined by 
the generalized angular velocity; the magnitude is 
determined by the coefficient J, and the phase, by the 
angle, 0. 

It is interesting to note that the current is passing 
through a maximum value at the moment that the 
current vector makes an angle are tan a/w with the 
projection axis. This occurs when the vector is in 
such a position that the angular velocity m, drawn 
from its terminus, is perpendicular to the projection 
axis. See Fig. 3. A steady sinusoidal current is 
maximum at the moment that the vector representing 
it lies in the projection axis. 

The behavior of damped sinusoidal currents in 
electric circuits, having resistance inductance and 
capacitance, has been the subject of considerable 
investigation, and is familiar, especially to those who 
have been interested in radio communication. The 
field of application of this theory, however, is much 
broader than this, and may easily be extended to include 


practically all electric machinery, both alternating and 
direct-current. 

Currents and pressures that have the same general- 
ized arİgular velocity can, of course, be combined by the 
customary principles of vector addition. The resulting 
vector diagrams will then be similar to those which 
apply to the steady state, save that the electromotive 
forces due to induction, both self and mutual, are not 
in time quadrature with the current, as we shall 
presently see. It is further understood that, as they 
rotate, all of the vectors shrink exponentially. 


P Axis 
P Axis 


0 Axis 


Fic. 2 


0 Axis 
Fic. 3 


Inasmuch as those who are interested in the opera- 
tion of electric machinery have had little contact with 
the theory of damped sinusoidal currents, it may be 
well to show how these currents and the electric pres- 
sures they produce can be represented by vectors. 
In all of the following analyses any capacitance that 
the circuit may possess will be disregarded. This, of 
course, places a definite limitation upon the scope of 
the analysis and prevents the discussion of certain 
transient conditions, particularly those of relatively 
high frequency, which are very important in trans- 
mission circuits. 


Fra. 4 


When a current of the general form flows through a 
non-reactive resistance the fall in potential is r I €” 
sin (wt + 0). This may be represented vectorially as 
in Fig. 4. The fall in potential through an inductance 
is: ۱ 

L Z = — a L I c sin (wt + 0) 
+ w L I e“ t cos (wt 0) 
The first of these components may be represented by a 
vector in opposition to the current vector, and the 
second by a vector which leads the current by 90 degrees. 
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See Fig. 4. In this figure and hereafter a line drawn 
above a letter Indieates that it represents a current or 
pressure of the genera] form. 

Thus: 7 w) LI 

or = (r + m L) Í 
where m — the generalized angular velocity of 
the current vector; the generalized impedance is 
(r + m L). 

In Fig. 4 notice, that there is a rise in potential 
through the circuit, numerically equal to a L I, which 
is in time phase with the current. This indicates that 
eleetrie energy is being generated within the circuit, due 
to the diminution of the magnetic field. The power so 
generated, œ L 12, partially supplies the copper loss, 
P r, in the circuit. The vector representing the fall in 
pressure due to self inductance leads the vector repre- 
senting the current by arc tan a/w more than 90 deg. 
It is in this respect that the vector diagrams repre- 
senting the transient state differ from those that 
represent the steady state. 

If the electric circuit is linked with a magnetic circuit 
in which eddy current and hysteresis losses may occur, 
it is customary, though strictly speaking it is incorrect, 
to represent both the current and the magnetic flux by 
vectors when they have a steady sinusoidal variation. 


Flux 


Current 


Fie. 5 


Fia. 6 


In this case the flux is represented as lagging the current 
which produces it by a constant angle 8. See Fig. 5. 
This is equivalent to assuming that the hysteresis loop 
is an ellipse and that the eddy currents are all in time 
quadrature with the flux in the core. The fall in 
reactive pressure J w L I, leads the flux, o, by 90 deg., 
but the current by an angle which is somewhat less. 
This being the case, it is necessary to represent the 
inductance by a complex number, i. e. L = L,—j Lə, 
where L/L, = tang. Thus] w L I = j o (L, — j L.) I 

= ور ده‎ I + jo L, I 
Similarly if the eurrent is of the general form,—and we 
make the same assumption in regard to the constaney 
of the phase angle between the flux and the current, — 
the fall in pressure due to inductance is (— o + w) 
(Lı — j Lə) I, which may still be written m L I, under- 
standing now that both m and L are complex numbers. 
On expansion, this becomes [(— a L, + o Lə) + j 
(w Lı + a Lə) ] I. If a L, equals w L, i. e. if a/w = 
tan 8, the fall in pressure due to inductance is in exaet 
quadrature with the current and the energy given up 


by the magnetic field due to its collapse is just sufficient 
to supply the eore losses. The power absorbed in the 
coil is then equal to the copper loss. If, with the same 
frequency, the rate of collapse of the magnetic field is 
greater than this, the coil absorbs less power than the 
copper loss. 

If there is another electric circuit linking the same 
magnetic core, the fall in electric pressure through it 
produced by the action of the current in the first circuit 
is j w M I, Notice that this is a vector which leads 
the current vector by more than 90 deg. When there 
are core losses in the mutual magnetic circuit the 
mutual inductance, M,, might likewise be represented 
by a complex number of the form, M, — 7 M.. 

Only considerable research can determine the value 
of this suggestion that the self and mutual inductance 
be represented by complex numbers when the currents 
are decaying sinusoids. It should be noted that 
there is a constant lag angle of flux with respect to 
magnetizing current only if the hysteresis “loop” is a 
diminishing ellipse such as shown in Fig. 6. This is, 
of course, not true and whether it is justifiable can 
only be determined by experiment. 

In any inductive circuit the current cannot immedi- 
ately respond to an abrupt change in the magnitude 
or form of the applied electromotive forces. There is 
always a transient condition of current flow, the dura- 
tion of which is determined by the inductance and 
resistance of the circuit. Moreover, the form of this 
transient flow is determined solely by these constants 
and is wholly independent of the character of the 
applied electromotive forces. The magnitude and 
phase of the transient, however, are determined by 
the initial and final values of the current in the circuit, 
and these are fixed by the applied electromotive forces. 


Two CIRCUITS HAVING MUTUAL INDUCTANCE BUT NO 
MOTION WITH RESPECT TO EACH OTHER 


The simplest case we shall consider is that of a single- _ 
phase transformer which is simultaneously short-cir- 
euited on both primary and secondary sides. 

With the primary winding short-circuited the fall 
in pressure through it is zero. 

Lnthmb,+i.mM = O (1) 
r, and L, represent the resistance and self inductance, 
not leakage inductance, of the primary. M is the 
mutual inductance between primary and secondary. 
Likewise: 

Irs + mL, + Í, m M ع‎ 0 (2) 
The primary and secondary currents have the same 
form and thus the same generalized angular velocity. 


Eliminating Î, and Í, from equations (1) and (2) gives: 


(ri + m Lı) (r, + m L,) = m? M? 
Expand and divide each term by L, L: giving, 
g m? + (k, + kə) m + ki ke = O 


M? 
where c = 1 — ps 


(3) 
7 = ri/Li and ke = ٣/ L, 
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By some writers c is called the leakage coefficient, 

whereas others define the leakage coefficient as 
L, L, 
M? 


L 
— 1. Inthe latter case ( 


, 


T o ۱ 
اوس[‎ 0۲ = “Tag Š In a transformer in which 
the resistance of the windings may be neglected this 
leakage coefficient, ø, is the ratio of the current on 
open circuit to that on short circuit with the same 
applied potential. kı and k. are the reciprocals of the 
time constants of the two windings on open circuit. 


The solution for the generalized angular velocity is 


ki + k: v (k, + k.)2— Ao k, ke 
20 2g 


Since ø is less than unity, the values of m are both real, 
and the transient current consists of two direct-current 
components one of which diminishes much faster than 
the other. These currents may be represented by 
vectors as in Fig. 1B. 

In this case, in which both windings are short- 
circuited, there will be no current in either winding 
after the transient has disappeared. Thus the only 
current after short-circuit occurs is the transient. Let 
the vectors I, and Jo represent the primary and 
secondary currents before short circuit occurs. Also 
let I,’ and I,” represent the components of the transient 
eurrent in the primary and similarly for the secondary. 
'Then at the moment of short circuit, the sum of the 
instantaneous values of the two transient components 
must equal the instantaneous value of the initial 
current in both primary and secondary. The transient 
components might thus be represented by an indefi- 
nitely great number of vectors since it is only necessary 
that their projection on the given axis have specified 
values. It is more convenient, however, if they are 
chosen so that their vector sum is equal to the vector 
representing the initial current before short circuit. 
Thus we will write: 


I| = İy + ۳ (5) 

20 — Ls’ 十 I," (6) 
Remember that the vectors J,’, ete., are diminishing 
exponentially and that this equality holds only at the 
moment of short circuit. Further, notice that the 
vectors ور‎ and Jo) are drawn to represent the maximum 


values of the currents before short circuit. Each of 
these transient components must satisfy equations (1) 


m = — 


(4) 


and 


and (2). Substitute in equation (1) and write in the 
following form: 

- (rı +m’ Dy) R 

İV MT FRE I! 

- (r, +m” L.) 2 
also 1 — مک پک‎ pP 


where m“ and m” are the two values of m as indicated 
by equation (4). If these equations are added and 
(Jio — Ii!) is substituted for I,” the solution for J,’ is: 


I,’ = 


kı (m” m”) { (ki/m" zs 1) Iio 
+v 1 - o V L:/ Lı In | (7) 
m" m" 
also I,” = 7 (m” = m’) { (ki/m' + 1) Ii 
+ V1— © V L,/Liləo) (8) 


The product of the two roots divided by their difference 
15 
m" m" m ki kə 
10 m 7 V (ki + kə) — 4 ©» ki ke 


In a transformer the ratio of the self-inductances of 
the primary and secondary 15 practically equal to the 
square of the ratio of the turns. Thus as is commonly 
said, “بد‎ L;/L; I2 is the secondary current referred to 
the primary. In order to reduce these expressions to 
a more understandable form, we will assume that k, 
equals .وم‎ This is very nearly true. If this is done the 
two components of the primary current become: 


77 “ə 1/2 (Lio + Iz) (9) 
IL E 1/2 (Ji ə Tə) (10) 


where I+, is now the secondary current referred to the 
primary. The values of m reduce to: 


NEN T 

m =- TM 
mn _ EN RM 
m = L-M 


where r and L are the resistance and self-inductance 
of the primary winding, and M is the mutual inductance 
on the assumption that both windings have the same 
number of turns. In this case the difference of the 
self and mutual inductances is the leakage inductance 
of one winding. The two components of the transient 
secondary current are similar to those of the primary. 
Thus at the moment of short circuit each of the currents 
in the two windings breaks up into two components, 
one of which is one-half of the no-load current; (I, 
+ 120) is the no-load current. This component decays 
at a rate determined by the resistance and the sum of 
the self and mutual inductances of the two windings. 
It is due to the energy stored in the magnetic core of 
the transformer. Its rate of decay is much slower than 
that of the other component. The latter’s rate of 
decay is determined by the resistance and leakage 
inductance of the windings, and this component of 
the currents is due to the energy stored in the leakage 
field. The vector diagram and oscillograph record 
are given in Fig. 7. In this figure the vectors J, and 
I», representing the primary and secondary currents, 
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are rotating at an angular velocity of w. At the 
moment of short circuit each of these vectors abruptly 
stops and breaks up into two component vectors which 
diminish at different rates. 

The more usual transient condition would be one in 
which the short circuit occurred on the secondary side 
alone. In this case the current immediately after 
short circuit consists of the transient current together 


with the steady value of the short-circuit current. The 
conditions that must be satisfied are now: 

Lio = ey + 77 + [se1 (11) 
and Izo = Í; -+ 1” + 5: (12) 


I, 


j^ Moment of 
Short Circuit 
7 
— 
p 
L-I 一 -一 A Secondary 
1 42 سے‎ m” 2 Current 
Pd \ A 
m 
کے ———— سس کر ےے‎ —”-Primary 
ly \ m” n Current 
1 | 
J 
Fia. 7 


where I,.ı and I,. represent the steady values of the 
short-cireuit currents in the primary and secondary 
windings. As before this is a vector equation which 
holds only at the moment of short circuit. If these 


steady short-circuit currents are transferred to the. 


left of the equality sign it is readily seen that the 
values of the components of the transient currents 
become: 


= 1/2 (Lio T Izo) zü 1/2 (Ee + İ.cə) (13) 
I” = 1/2 (Lio 一 120) — 1/2 (He — I (14) 


As before the secondary components are similar to the 
primary components. It is now necessary to determine 
the phase relation between the steady short-circuit 
current and the current in the transformer before the 
short cireuit occurs. 

In order to keep the adamı as simple as possible, 
we will make the customary assumption that the 
primary and secondary currents are equal and in phase 
opposition. In that case the first component of the 
transient. current is zero, and the second component 
is equal to the current before short circuit minus the 
steady short-circuit current. [See equations (13) 
and (14) ] with 10 = — ET and İle = 一 | re The 
error thus introduced is usually very slight. The 
vector diagram at the moment of short circuit is given 
in Fig. 8. 

Before short circuit occurs only the four vectors 
representing the primary and seeondary potentials and 
current appear. They are rotating at a uniform 
angular velocity, w. At the instant after short circuit 


se2) 


all of these except the primary potential disappear and 
in their place spring the steady short-circuit currents. 
I,a and I... and the transient currents I,” and ۰ 
The former continue to rotate at an angular velocity of 
w without change in magnitude, while the latter remain 
fixed in angular position but diminish exponentially at a 
rate, a, depending upon the resistance and leakage 
inductance of the windings. The arrows at the ends of 
the vectors indicate the direction of their motion. 
Observe that, at the moment of short circuit, the current 
before short circuit is the vector sum of the steady 
short-circuit current and the transient current. Since 
the steady short-circuit current is from 10 to 25 times 
the full-load current, it makes little difference whether 
the short circuit comes at a time of full load or no load. 
It may be a little more severe at no load. On account 
of the excessive value of the steady short-circuit current, 
the vectors representing it and the transient current are 
essentially in phase opposition at the moment of short 
circuit. Half a cycle after short circuit occurs the 
vector بر‎ will thus have overtaken the stationary 
vector I,” and their vector sum will be a maximum. It 
at once follows that short circuit is most severe when it 
occurs at a moment that the steady short-circuit current 
would be a maximum.” Since the steady short-circuit 
current and applied potential are nearly in quadrature, 
for example, the power factor of a certain 1000-kv-a., 
60-cycle transformer is 0.15 at short circuit,—the worst 
moment for a short circuit to occur is at about the time 
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that the applied potential is passing through zero. In 
the transformer just cited the ratio of the resistance to 
the leakage inductance is 56 and thus in the time of one 


half cycle the transient current will have fallen to 

56 
e 120 or 63 per cent of its initial value. Thus the 
first current rush would be to a value which is 63 per 
cent greater than the maximum value of the steady 
short-circuit current which is itself more than 20 times 


the maximum value of the full-load current. The 


2. To determine the exact moment at which the short-cir- 
cuit current would be a maximum involves the solution of a 
mixed trigonometric and exponential equation, which is decidedly 
not worth while in the present problem. 


6 LYON: ELECTRIC MACHINERY 


forces between turns vary as the square of. the current, 
and thus the first shock on this transformer may be 265 
per cent of its ultimate maximum value. 

When these forces have their maximum value they 
occur one-one-hundred-and-twentieth (1/120) of a 
second after short circuit in a sixty-cycle transformer 
and one-fiftieth (1/50) of a second after short-circuit in 
a twenty-five-cycle transformer. Nevertheless with 
transformers having the same short-circuit characteris- 
ties (volts, amperes, and watts on short-circuit test), the 
maximum forces that the windings must sustain are the 
same for all frequencies. 


If the transformer is but a portion of the electric 
circuit in which the transient takes place the character- 
istic equations (1) and (2) must be written: 


ln + 7, m Li + L, m M = V. (15) 
I, و۲‎ + I, m L, + I, m M = و7‎ (16) 


V; and V; are the transient potentials which appear at 
the terminals of the transformer. If the transformer 
has a relatively small exciting current little error will 
be made if it is treated like a simple inductive reactance. 
In this ease the relation between the primary and 
secondary pressures is much simpler. 


Vi— Il, ) + m L) = a V, (17) 


R and L are the equivalent resistance and leakage 
inductance of the two windings and g is the ratio of 
transformation at no load. 


Two CIRCUITS HAVING MUTUAL İNDUCTANCE, IN ONE 
OF WHICH THERE IS A COMMUTATED ROTATIONAL 


Before turning our attention to polyphase rotating 
machines, it might be well to consider briefly the direct- 
current dynamo of either the shunt or compound type. 
This apparatus is in some respects quite similar to the 
transformer. It has two circuits which are mutually 
inductive, the armature and series field forming one and 
the shunt field the other. There is this difference, 
however: an electromotive force is generated in the 
armature circuit by its rotation in a magnetic field that 
is due to the joint action of the current in both circuits. 
With constant speed this generated e.m.f., is pro- 
portional to the magnetic field. The relation between 
the magnetic field and the currents in the armature and 
shunt field circuits which produce it is similar to the 
well-known hysteresis loop, for which there is no simple 
mathematical expression. In some problems in which 
the exciting current is a steady alternating one, it has 
been found convenient, however, to assume that this 
hysteresis loop is an ellipse. If this assumption should 
prove sufficiently accurate in the present case also, it 
will simplify the problem greatly, for as has been 
already suggested, we may then represent the magnetic 
field and the rotational e. m. f. which it produces, to- 
gether with the net exciting ampere-turns by vectors 
even when these quantities are damped sinusoids, 9. e., 


of the general form. This vector relation is shown in 
Fig. 9, which is drawn for the steady state. 

The flux, ó, and the rotational e. m. f. E which is 
proportional to it, lag the exciting current I by a con- 
stant angle 8. Thus we may write 

or E = KI 
where K is a complex number, K, — 7 Ke, having the 
dimensions of resistance. The vector equation shows 
that when the instantaneous value of the current is zero 


(18) 
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the generated e.m.f. is KI. This is generally 
referred to as the residual voltage. At the moment . 
that the current has its greatest value the rotational 
e.m.f.is K, Z. If the exciting current is going through 
a steady cyclic variation, it is a simple matter to 
determine these two components of the e.m.f. See 
Fig. 10. When the current is a damped sinusoid, how- 
ever, the vector representing it is shorter at the moment 
that it is perpendicular to the time axis than it was 
earlier in the cycle when the current was passing through 
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its maximum value. See Fig. 11. In this case the best 
method of determining the constants K, and Kə, in 
order to assure the greatest accuracy in the calculations 
can only be determined by experiment. 

Let I, and I, be the currents in the armature and 
shunt field circuits and similarly, rı and Li, the resist- 
ance and self-induction of the one, and r: and Lə, the 
same constants of the other. Also let M be the mutual 
inductance between the two circuits. The e. m. f., E, 
generated in the armature circuit by rotation is: 

E = K (Ni/Ne L + I»); 
where N, and N; are the effective turns in each circuit 


. (Ka— rı— m Lı) (r: +mL)— (K—mM)mM 
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acting on the main magnetic axis. Let N;/N; = a. 


Also let V be the rise in terminal pressure across both 


the armature and shunt field circuits. 
Then we have: | 
V = Ë— Lr,— L m L, — وآ‎ M 
(armature circuit) 
= Îl, (K a — r, — m Li) + L, (K — mM) (19) 
also V = — İ,r, — Í, m L, — I, m M 
(shunt field circuit) 
= — L, (rs + m L.) — 7, m M (20) 
_ If these equations, (19) and (20), are solved for T, and 
I, we have: 
I, = 
T; + m ور‎ + K — m M = 
(K a — ۲۱ — m Lj) (rs + m Lj) — (K — m M) m M ۷ 


and I, = 
_ K a— r, — m L, + m M = 
(Ka — r, — m L.) (rs + m L) — (K —m M) m M ۲ 


_ The total net current delivered, Î = T, + Îs, is 
[= 


۲۱ + rs + K (1 — a) + m (L. + L, — 2 M) V 


The ratio of the fall in terminal potential, — V, to the 
total current is the generalized impedance, Z, of the 
dynamo. Thus we have: 

Z = 


(K — m M) m M — (K a — r: — m Lı) (r; + m L) 
rı + r; + K (1 — a) + m (Lı + ور‎ — 2 M) 


For further discussion of this problem the reader is 
referred to the paper by Doherty on “Exciter Instabil- 
ity” presented at the Pacific Coast Convention, August, 
1922. 


Two SYSTEMS OF BALANCED POLYPHASE CIRCUITS 
HAVING A CONSTANT RELATIVE ANGULAR VELOCITY 


The author believes that the simplest method of 
presenting the analysis of the transient condition in 
polyphase induction and synchronous machines is that 
which considers the rotating magnetic fields in the air 
gap and the reactions which they produce. When these 
polyphase machines are operating under balanced 
conditions, either steady or transient, the currents in 
the several phases of the polyphase winding may be 
represented by vectors of equal length equally spaced 
with respect to each other. For the steady condition 
the vectors are of constant length and rotate at a 
constant angular velocity. This representation has 


been used for many years. The only difference between | 


the steady and the transient condition is that the vec- 
tors which represent the latter state diminish exponen- 
tially as they rotate. If this conception of the transient 
state is borne in mind the following analysis is readily 


followed by anyone who is familiar with the steady 
operation of these machines. In all that follows the 
relative speed of the two windings is assumed to re- 
main constant throughout the time considered. Since 
the maximum current occurs within about a half- 
cycle after the transient begins the angular velocity 
of the moving part has little opportunity to change 
before the first shock is over. 

The parallel lines in Fig. 12 represent the uniform 
air gap of a three-phase induction motor. The phase 
belts of the stator conductors are indicated at the letters 
a, b, c. The pole pitch is from +a to — o. If this 
winding is properly designed, there will be very small 
harmonics in the space distribution of the air-gap flux 
density when the alternating currents are balanced. 
We will neglect whatever harmonics may exist and 
assume that this distribution may be represented by a 
sinusoidal curve such as 8. The line X — X is the 
axis of this curve. If the current in phase b lags that 
in phase a by 120 deg. this magnetic field will move 
toward the right at à constant angular velocity of o 
radians per second; where w equals the frequency 
of the alternating current multiplied by 2 v. On the 
other hand, if the current in phase b leads that in phase 


a the field moves toward the left at the same velocity. 
In all that follows we shall consider that a velocity 
toward the right is positive and that one toward the 
left is negative. Furthermore, the magnitude of the 
field as well as its form, will not change as it moves, 
and thus the field might be represented by a rotating 


veetor whose projection on a fixed axis indicates the 


amount of flux linking any particular phase. The 
angular displacement of this vector from the fixed axis 
is always the angle between the axis of the moving 
field and the axis of the phase considered. See Fig. 
13. One of the most important points in the theory of 
rotating fields that are produced in this manner in a 
uniform air gap is that the flux through any phase is a 
maximum at the moment that the current in that phase 
is also a maximum. For example, when the current 
in phase a is maximum the field is in the position 8.. 
One sixth of a cycle later, the current in — c will be a 
maximum and in the same time interval, the flux 
distribution will have moved one sixth of two full pole 
pitehes toward the right and be a maximum through 
this phase, — c. In a similar manner, it will have 
moved into the position 8, by the time that the current 
in phase b is à maximum. This can be represented 
vectorially as in Fig. 13. The current and flux are 
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represented by the rotating vectors I and ¢. The 
instantaneous value of the current in any phase or the 
flux through that phase is the projection of the cor- 
responding vector on the axis of that phase. In Fig. 
13 the vectors are drawn for the time that the current 
in phase — c has its maximum value. Notice that for 
both the phases a and b the currents in them and the 
fluxes through them are each one half of their maxi- 
mum values. There is this difference, however, 
for while the values are decreasing for phase a they are 
increasing for phase b, in which they will be a maxi- 
mum 1/6 of a cycle later. Ordinarily, we are content 
to know the current and flux in any one phase, a, for 
example, and then the projection axes for the other 
phases are not drawn. 


Axis-phase a 


Axis-phase b Axis-phase ¢ 
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It may be convenient and possibly necessary in some 
problems to consider magnetic fields that rotate to the 
left, that is, with a negative angular velocity. In this 
case the vector, ,رب‎ representing the flux through any 
phase would likewise have a negative angular velocity 
and would thus move in the clockwise direction. If 
the trigonometrie expressions for the instantaneous 
currents in the phases are written it will be noticed that 
giving a negative value to w is equivalent to reversing 
the phase order of the currents. Such a reversal of 
phase order reverses the direction of rotation of the 
magnetic field produced by the current. If the current 
vector in Fig. 13 is rotated in the negative or clock- 
wise direction it will be noticed that the current in 
phase b will lead that in phase a by one-third of a cycle, 
that is the phase order of the currents is reversed. We 
must be most careful, however, to remember that, if it 
is necessary to combine or operate on these vectors, 
whether they have positive or negative angular velocities 
the operator 7 always indicates a phase rotation of one 
quadrant in the positive or counter clockwise direction. 
If this is not done we may become hopelessly confused. 
If the air gap is not uniform, as is the case with salient 
pole machines, the curve of flux density will alter both 
its shape and magnitude as it rotates, but that is a 
complieation which it is well for the present to disregard, 
leaving its consideration until some future time. 

In general if the currents in the several phases of 


a symmetrical n-phase winding have values equal to the 
projections of a “current” vector on each of a system 
of n equally spaced axes, the magnetic flux linking any 
winding is equal to the projection of a “flux” vector, 
which is drawn in phase with the “current” vector, on 
the axis corresponding to that winding. This assumes 
that the windings are so designed that there are essen- 
tially no harmonies in the space distribution of the 
magnetic field and that each winding will, if acting 
alone, produce exactly the same flux per ampere. 

If the currents are of the general form both the 
“current” and “flux” vectors shrink exponentially as 
they rotate. The space distribution of the magnetic 
field is still sinusoidal at all times but its magnitude 
diminishes as it moves through the air gap. Compare 
Fig. 18 with Fig. 14. In the latter figure let us assume 
that the currents took on the general form at the moment 
that the flux distribution was at the point indicated by 
the curve 8, and that their generalized angular velocity 
is (一 a +j w). The equation of the maximum density 
for any later position of the field suchas 8115: B, = Bee”, 
where t is the time required for the field to move from 
the position 8 to 8), 2. e., through an angular distance X. - 

The flux through any phase as vvell as the current in 
it may still be represented by a vector although this 
vector is now of the general type—that is, it shrinks 
exponentially as it rotates. A current or flux variation 
of the general form reaches its maximum value in any 
phase at the moment that the veetor which represents 
it is still are tan a/w from the projection axis of this 
phase. Because the field is shrinking as it rotates 
the flux linking phase a is greater when the field is 
in the position 8: than when it is in the position رم‎ at 
which moment the axis of the field is coincident 


with the axis of the phase. See Fig. 15. Thus the 
flux through any phase is a maximum at the moment 
that the current in the phase is a maximum, but the flux 
distribution is no longer directly opposite the phase 
when this occurs, as was the case when the currents were 
steady sinusoids. 

It is then evident that even in the transient state the 
flux through any phase and the current in that phase 
may still be represented by coincident vectors. 

The electromotive force generated by this rotating 
field is the negative of the time rate of change of the 
flux linkages. "Thus, in general, when this time rate 
of flux change in any phase is zero the e. m. f. is zero 
and the flux linking the phase is a maximum. In Fig. 
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16 the flux in phase a is a maximum when the vector 
representing it is in the position shown. At this 
moment the e. m. f. is zero and the vector, Ë, repre- 
senting it must, therefore, be perpendicular to the pro- 
jection axis. An additional proof of this is interesting 
even if unnecessary. Consider Fig. 15. When the 
field is in the position 6, the e. m. f. produced in the 
phase a has two components, one due to the motion of 
the field and the other due to its shrinking. The 
motional e. m. f. is proportional to the velocity of the 
field and to the sine of the angle, w t, between the field 
and the phase winding. If at the poimt + a, the 


m 


0 


Axis of Feld لح‎ | — 
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direction of the flux density is upward and the motion of 
the field is to the right the e. m.f. generated in the 
belt, + 0, 15 away from the observer. On the other 
hand if the field within a circuit changes an e. m.f. 
is set up which is proportional to the rate of change and 
is m such a direction as to oppose the change. The 
e. m. f. thus produced is proportional to a cos w t, and 
is toward the observer in the belt of conductors + a. 
If the angular displacement, w t, is so chosen that the 
resulting e. m. f. is Zero, w sin wt equals a cos o f. 
That is to say, the angular displacement of the flux 
distribution is are tan a/w. Observe that the fall in 
potential — E, leads the flux, ¢, by an angle, 2 
+ arc tan a/w, just as does the fall in potential due to 
self induction. "Thus, in a polyphase winding we may 
combine the fall in pressure due to the so-called leakage 
flux in one phase with the fall in pressure in this phase 
due to the mutual air-gap flux which is produced by all 
phases. "This is commonly done in the analysis of 
synchronous machines and the resulting loss in pressure 
is the so-called synchronous inductance, or reactance 
drop. If L represents the synchronous inductance, 
a transient current, Î, having a generalized angular 
velocity m, will produce a fall in potential of m LI. 
This fall in potential is represented by a vector which 
is 7/2 + are tan a/w ahead of the current vector. 
This conception of the action of a rotating field of 
the transient type may be summarized as follows: 
Consider two groups of symmetrical polyphase wind- 
ings, a and b, the first of which is stationary and the 
second, moving in the positive direction at an angular 
velocity of » radians per sec. The projection axis for 


the stationary winding is stationary, while the projec- 
tion axis for the moving winding is moving at the same 
angular velocity, viz. p, as indicated in Fig. 17. If the 
angular velocity of the currents in the stationary, or 
a, winding, is w radians per sec., the relative angular 
velocity of the magnetic field due to these currents 
with respect to the b winding is (w — p) radians per 
sec. in the positive direction. The flux linking any 
phase of the b winding is a maximum are tan 

- = 7 before the flux vector lies in the projection 
axis of this phase. If the field is moving with respect 
to the b winding in the positive direction the flux will 
be a maximum when the flux vector is at ."وه‎ See 
Fig. 17. If the field is moving relatively in the nega- 
tive direction, z. e. if p is greater than w, the flux will 
be a maximum when the flux vector is at 6.”. At the 
moment that the flux is a maximum through any phase, 
the generated e.m.f. in this phase is zero and the 
vector representing it will be at E,' for a positive 
relative angular velocity and at E,” for a negative 
velocity. This generated e. m. f. is equal to the pro- 
duct of the current in one phase of the a windings, the 
relative generalized angular velocity with respect to 
the b windings of the field that these currents produce, 
and the mutual inductance between the a and b wind- 
ings. The mutual inductance M can be measured as 
follows: Fix the a and b windings with respect to each 


Jp 


272 


r 53 
arc tan DÙ 
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other and send balanced polyphase currents through 
either, the a for windings for example. The mutual 
inductance M is the e. m. f. generated in one phase of 
the other, or b winding, divided by the product of the 
current in one phase of the a winding and angular 
velocity of this current. By the time that the flux has 


advanced 7/2 — arc tan with respect to the 


b windings, the generated e. m. f. will have reached its 
maximum value, but the current will usually not reach 
its maximum value until somewhat later depending 
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upon the eonstants of the b windings and their terminal 


pressure. In any case, however, the current will be 
Ë,- V, . 

represented by a vector equal to ———— — where 
Tə + m L, 


V, is the terminal pressure, r: the resistance and L, 
the synehronous inductance of one of the b phases. 
By the time the current in the vinding has reached its 
maximum value the field will have advanced a further 
relative angle of 0. Thus whether the motion of the 
a field is positive or negative, its action on the 6 wind- 
ing will develop a field which will always lag the a 
field in its motion with respect to the b winding. Con- 
sider an induction machine in which the magnetic 
field due to the stator, a windings, rotates in the 
positive direction at w radians per sec. If the rotor 
turns in the same direction at less than this synchronous 
speed the relative motion of the stator field with respect 
to the rotor is positive. The rotor component of 
magnetic field actually and apparently lags the stator 
field. If the rotor turns in the same direction at a 
speed greater than synchronous, the relative motion 
of the stator field with respect to the rotor is negative 
and the rotor component of field apparently leads when 
observed from the fixed stator, though it actually lags, 
as would be observed from any point on the rotor 
itself. It is interesting that these well-known facts 
should apply to the transient state as well as to the 
steady operating condition. 

Having established these general principles in regard 
to rotating fields during the transient state the sqlution 
for the transient currents in an induction motor may 
be written down without further explanation. If the 
generalized angular velocity of the stator currents is 
m the relative generalized angular velocity of the field 
produced by these currents with respect to the rotor 
is (m — ? p). Since m = — e + 7 o and the relative 
angular velocity of the stator magnetic field and the 
rotor winding is (o — 2), the generalized angular 
velocity of the rotor currents is — a + 7 (€ — p), 1. e. 
m — j p. The generalized angular velocity of the rotor 
currents is thus m — j ,م‎ but the relative angular veloc- 
ity with respect to the stator of the field produced by 
them is İm — jp — (— 3 p) ], i.e m. The second 
(— j p) is the relative angular velocity of the stator 
with respect to the rotor. Thus we have: 


Ë, = — (m — jp) M I, (19) 
Ë, = — m M و[‎ (20) 


E, and E, are respectively the generated e. m. f. due 
to the stator currents in the rotor and that due to the 
rotor currents in the stator. The rotor current, İş, is: 

I 2 = E 9 / Z 9 (21) 
where Z; is the generalized self impedance of the rotor 
calculated for a generalized angular velocity of m 
— jp. The stator current J), is: 
E,— Vi 


İ, = 7 (22 


where Z, is similarly the generalized self impedance of 
the stator calculated for a generalized angular velocity 
of m, and V, is the transient terminal pressure of the 
stator. Combining these four simple relations gives 
the characteristic equation for the stator current of 
an inductor motor during the transient period. That 
is: 


p Z; E 
h^-ZZ-mm-jpwm ^ ۰. ۵ 
_ = V. 
Zt 


where — V, is the fall in terminal pressure across one 
phase of the stator. That is, it may be said that the 
transient impedance of an induction machine is: 


m Z, و2‎ — m (m — jp) M” 


Y= Z, (24) 
The steady operating impedance 5 
71 2 — j ) w — j p) ۰ 
z, - 2% بط طم‎ J p) (25) 


Zə 


where Z, = r) + 7 o Lı and Z: = r? + (o) — p) La. 
If the transient occurs when the stator winding is 
symmetrically short-circuited or if it occurs when the 
stator is connected to a circuit whose resistance and 
inductance are many times smaller than the similar 
constants of the motor, the transient terminal pressure 
of the stator is essentially zero. It is thus evident 
that the transient impedance is also zero. Equating 
the transient impedance Zr to zero gives the charac- 
teristic equation for an induction motor and determines 
the generalized angular velocity of the stator current. 
Since Z, = r: + mL, and Z: = r; + (m — j p) Ls 
equating the transient impedance to zero gives 
(rı + m Li) [rs + (m — j p) Le] - m (m — jp) M? =O 
Expanding and collecting terms and dividing by Lı La 
gives the following quadratic in m: 


o mö + (kı + دوم‎ jpo)m + وتا ص 7 ح و رط‎ (26) 
M? 
where kı = r /L,, k, = r,/ Ls ando = 1— 
1 Le 


Notice that this characteristic equation for the induc- 
tion motor is identical with that for a transformer with 
the exception of two terms, — pom, and —j p kı 
These terms, hovvever, completely alter the form of the 
transient. The solution for the generalized angular 
velocity, m, is: 


+ ip 
2 


“m 20 + 
ə Se a 


If k, = k,, the quantity under the radical is real and 
will be negative unless the equivalent resistance is 
somewhat greater than the equivalent reactance as 
measured with the rotor blocked, a most unusual con- 


(27) 
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dition except possibly in very small motors. Fur- 
ki + ke . 
thermore the decrement — E is very nearly equal 


to the resistance of the stator winding divided by the 
equivalent inductance, L,, as determined from locked 
saturation test. Thus under ordinary conditions the 
two values of the generalized angular velocity are 
approximately: ۱ 
m" سے‎ — r,/Lo +- 3 ( (28) 
m" = — r,/Lo + 7 0 (29) 

These are the values given by Doherty and William- 
son. The magnitude of the error is best shown by a 
numerical example. The relative constants of a 25- 
cycle, 2600-h. p. induction motor are rı = r: = L. 
= Lə and ø = 0.075. The equivalent leakage induct- 
ance as calculated from these values is relatively L, 
= 0.0751. The approximate values of the generalized 
angular velocity for a speed of 157 radians per sec. are: 

m" = — 13.31 + 7 157 
m" = — 13.31 + 0 

The values as caleulated from the solution of the 

characteristic equation are: 

m’ = — 13.33 + 7155.8 

m” = — 13.33 + 71.2 
The error in the approximate values in this case at 
least is very small. 

The form of the generalized angular velocity shows 
that the transient current during short circuit consists 
of two components, one having a frequency slightly 
below that corresponding to the speed of the motor, 
and the other a very low frequency. Notice that the 
sum of these frequencies always equals that correspond- 
ing to the speed of the motor. If the ratios of resis- 
tance to synchronous self inductance are not the same 
for the stator and rotor windings the components 
of the transient current will decay at different rates, 
otherwise their rates of diminution are ordinarily equal. 
During this transient short circuit the angular velocity 
of the rotor is greater than the angular velocity of the 
stator currents so that the machine is operating as an 
induction generator. In this case, the swm of the angu- 
lar velocities of the corresponding stator and rotor 
component currents is equal to the angular velocity 
of the rotor. In calculating the performance of an 
induction machine under steady operating conditions 
it is customary to use an effective resistance of the 
stator winding and a true, or ohmic, resistance of the 
rotor winding. Similarly, it may prove more accurate 
to use an effective resistance of stator and a true resist- 
ance of rotor when calculating m’, and, vice versa, 
a true resistance of stator and an effective resistance 
of rotor when calculating m”. If, for example, the ratio 
of effective to true resistance for both stator and rotor 
is 1.5, the rate of decay of the transient currents is 
about 25 per cent greater,? than would be calculated 


3. A.I. E. E. Transactions, page 525. 


by using the true resistances. Again, if the ratios of 
true resistance to synchronous self inductance are the 
same for both stator and rotor, the rates of decay of 
the two transient components would be slightly dif- 
ferent if effective resistances are substituted for true 
resistances as suggested. ۱ 
Having determined the form of the transient it 
remains only to determine its magnitude. The funda- 
mental principle is that the current in any inductive 
elreuit eannot abruptly change its value. Thus, if 
the currents in the stator and rotor windings immedi- 
ately before the transient occurs are represented by the 
vectors I, and I, in Fig. 18 the vector sum of the 
currents in these windings immediately after the 
transient begins must be respectively I, and ۰ 
The case is somewhat different from the single-phase 
transformer inasmuch as there are now two or more 
stator or primary circuits and a corresponding number 
of rotor or secondary circuits in which the instan- 
taneous values of the currents immediately before and 
after the transient begins must be the same. It has 
already been pointed out that for any mstantaneous 
values of balanced polyphase currents in symmetrical 
windings both the currents and the magnetic field 
which they produce may be represented by vectors 
in phase with each other. If the currents in the stator 
and rotor windings do not abruptly change at the 
moment that the transient begins, it is evident that the 
component magnetie fields they produce in the air 
gap do not change in magnitude or in position. That 
is, the veetors which represent the component fields 
produced by the stator and rotor currents do not change 
at the moment that the transient begins. Since the 
eurrent vectors are in phase with the flux vectors they 
also do not change abruptly. The general relation 
that holds at the moment that the transient begins is: 


Iş = İF +I +I, (30) 
2007” I,’ + I + Iz, (31) 


where 7 ور‎ and Izo are the vectors representing the stator 
and rotor currents before the transient begins; I,’ 
and I,” are the generalized vectors representing the 
two transient components of the stator currents; 
similarly I,’ and I,” are the generalized vectors repre- 
senting the transient components of the rotor current; 
I,, and I, are the vectors which represent the steady 
currents in the stator and rotor which will exist alone 
after the transient has disappeared. As it is more 
convenient if the rotor currents are determined on the 
basis of a one to one ratio of transformation this will 
be done throughout. In this case L, = Ln, at least 
approximately. The steady operating conditions be- 
fore and after the transient occurs are sufficient to : 
determine the steady current vectors مور مور‎ etc. 
The following relations must hold for the transient 
currents: _ _ 

I’ + I," 一 1: 一 li, 

i + i T 120 一 Ias 
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If an induction motor or generator is short-circuited 
on all phases J,, and Iz, are both zero since there will 
ultimately be no current in the machine. If an induc- 
tion motor or generator is connected to a line when it 
carries no current, Iio and Iz are both zero. We shall 
work out in detail only the first of these eases since 
nothing is gained in generality by doing more. 

Each of the transient. components must satisfy the 

. vector relation: 


(Ti + m Li) T + m M I, = O 
and thus we may write: 

(rı + m' L) İy + m' M L,' = 0 

(Ti + m" Li) 17 + m” M I." O 
also at the moment of short circuit 

I,' جات‎ I," = Ii 


(32) 
(33) 


(34) 


whether short-circuit occurs at no-load or full load. - 
The first of these components rotates at approximately 
synchronous speed while the second is practically 
stationary. Each shrinks exponentially at the same 
rate. The first maximum value of the transient cur- 
rent in any phase is greatest if it occurs at such a time 
that the transient components lie as nearly as possible 
in the projection axis for that phase at the moment 
that the transient begins. For, in that case, about 
one-half a cycle later the first component will have 
overtaken the second and the resulting projection on 
the axis of the phase will have its greatest value. Dur- 
ing this time, however, each of the components will have 


k x 
shrunk to € °? of their initial values. The first 


and İy, + IL = Lo 


winding gives: 


m m” 


m" 


Ty = و‎ [M /rı Iso + (1/ki 
+ 1/m”) Tal 
and I” = کل‎ [M/r, I» + (1/k 
m — m 


+ l/m/) Ii] 


These expressions are readily reduced to understandable e rg ls 
If kı i 1 


forms if suitable approximations are made. 
= وم‎ = k the values of m’ and m" are approximately 


— 16/6 + jp (36) 
— k/o + 30 (37) 
4/1-— oL 
= (1 — c /2) L approximately 

where L is the synehronous self inductance of either 
the stator or rotor winding. The approximate values 
of the transient components become: 

J^ = 110/2 一 l/o I, 

I” = [10/2 + 1/61, 
The error in these values is well within 10 per cent 
ordinarily. Notice that the transient currents depend 
almost entirely upon the no-load current, 7,, and the 
leakage coefficient, o. It makes little difference 


4. More exact solutions for the component currents, based 
on the approximate values of m are: 


I, -(1 til | (12 - SEE 1.) 


ll 


Also 


. R 
k zı: c p 
an “= (on) k I 10/2 
1 )— 
T3 2 
1 — 
m yV "op 


(35) 1 ii 
Solving these four simultaneous equations for the +. 
components of the transient currents in the stator 


jp 


maximum current rush may be to an approximate value 
of 


(38) 


where 7, is the maximum value of the exciting current 
before the transient begins. This approximate solu- 
tion checks with that given by Doherty and William- 
son although the latter is in a slightly different form. 
The vector diagram is given in Fig. 18, which is drawn 
at the moment that the transient begins. The pro- 
jection axes of one phase of the stator winding and of one 
phase of the rotor winding are shown. The former is 
fixed in position while the latter rotates at a speed 
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equal to that of the rotor, viz., p radians per sec. When 
s is the slip, p equals (1 — s) wo. Before the transient 
begins only the vectors representing the stator and 
rotor currents appear. They are rotating at a uniform 
angular velocity, wo. At the moment of short circuit 
both of these vectors vanish and in their places spring 
the vectors which represent the transient currents. 
The arrows at the ends of these vectors represent their 
angular velocities. 


SYMMETRICAL POLYPHASE AND SINGLE-PHASE SYSTEMS 
MOVING WITH RESPECT TO EACH OTHER AT A 
CONSTANT ANGULAR VELOCITY 


The analysis of this problem is somewhat more 
difficult than that of the preceding. When a polyphase 
winding is acted upon by a rotating magnetic field 
it reacts and produces another magnetic field which 
always rotates at the same speed and in the same direc- 
tion as does the first. This is not true of a single- 
phase winding, however, and it is herein that the dif- 
ficulty lies. When a single-phase winding is excited 
by a rotating magnetic field, there is produced in it 
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an alternating current of a frequeney corresponding 
to their relative angular velocities. An alternating 
current in a single-phase winding produces a magnetic 
effect that is equivalent to two magnetic fields which 
rotate in opposite directions with respect to the winding 
at an angular velocity equal to that of the current in 
the winding. This well-known method of analysis 
is ascribed to Ferraris. The magnitude of each of 
these fields is one-half of the maximum value of the 
actual alternating field produced by the single-phase 
current. 

We have already established the relation that 
exists between a rotating field and the field which may 
be produced by its action on a polyphase winding. It 
now becomes necessary to establish a corresponding 
relation when the rotating field acts on a single-phase 
winding. In Fig. 19 let .رف‎ represent a positively 
rotating magnetic field which has a complex angular 
velocity (— a + j w). The flux linking the single- 
phase winding will be a maximum at the moment when 
مرفي‎ is in the position shown. The generated e. m. f., 
E, in the winding will be zero at this instant. By the 


time that ğı, has advanced (7/2 — are tan a/w) this 
generated e. m.f. will be a maximum. The current 
however will not reach its maximum value until a time 
still later that is determined by the generalized self im- 
pedance of the single-phase winding. At the time that 
the current is a maximum the vector representing itis arc 
tan a/w from the axis of the winding as shown in Fig. 
20. At this moment one of the oppositely rotating 
fields موض‎ is coincident with the current and the other, 
$>, makes the same angle with but is on the other side 
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of the axis of the winding. Of course, the axis of the 
field due to a single-phase current always coincides 
with the axis of the winding. 

One of these oppositely rotating fields travels at the 
same speed and in the same direction as does the initial 
field, $1, while the other travels at the same speed 
but in the opposite direction. The condition at a time 
angle wt + 6 after the current in the single-phase 
winding is a maximum is shown in Fig. 21. The op- 
positely rotating fields are represented by the vectors 
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and $z: the subscripts a and b indicate the direc-‏ موقي 
tion of rotation. The constant phase relation between‏ 
will be determined presently.‏ موف and‏ مرق 

If the first magnetic field had been rotating negatively 
with a complex angular velocity of (— a — ? w) the 
flux linking the single-phase winding would have been 
maximum at the moment shown in Fig. 22. By 
the time the current had reached its maximum value 
this field would have reached the position indicated in 


14 LYON: ELECTRIC MACHINERY 


Fig. 28. The oppositely rotating components of the 
single-phase field are now shown at $> are .روفي‎ Ata 
time angle o t + 6 later the vectors shown in Fig. 23 
will have moved to the positions shown in Fig. 24. It 
is important to observe that the oppositely rotating 
components of the single-phase field are exactly the 
same in Fig. 21 and in Fig. 24. That is, as far as the 
reaction produced in the single-phase winding is con- 
cerned, there is no difference between positively and 
negatively rotating fields provided they are always 
conjugate with respect to the axis of the single-phase 
winding. In the case of polyphase windings we have 
seen that positively and negatively rotating magnetic 
fields produce quite different effects. 

A magnetic field that rotates positively with respect 
to a polyphase winding excites the latter so that it in 
turn produces another field which rotates positively 

and which lags behind the first by an angle determined 
` by the relative generalized angular velocity of the 


rotating field and the polyphase winding, and the 


resistance and synchronous self inductance of the 
latter. The reaction of the polyphase winding is 


Fic. 24. 


often analyzed in the following manner: The alter- 
nating current that is generated in each phase of the 
polyphase winding produces two component fields 
that are rotating in opposite directions with respect 
to the winding. Those component fields due to the 
different phases which rotate opposite to the exciting 
field sum up to zero if the phase circuits are sym- 
metrical. The other component fields that rotate 
in the same direction as the exciting field are in space 
phase with each other and their resultant is thus n 
times any one component; where n is the number of 
phases of the polyphase winding. That is, we may con- 
sider that the current in each phase of a polyphase 
winding consists of two equal components, hereafter 
designated by the subscripts a and b, each equal to one- 
half of the maximum phase current. These two compo- 
nent currents produce equal but oppositely rotating mag- 
netic fields. If the winding were single-phase each 
of these magnetic fields would exist, but with a sym- 
metrical polyphase winding only the fields exist which 
rotate in the same direction as does the exciting field. 
Those that rotate in the opposite direction neutralize 
one another, so that their resultant is zero. 


Let M, be the e. m. f. generated in one phase of the 
polyphase winding by one ampere of either the a or b 
component of the single-phase current when varying 
at the rate of one ampere per second. Similarly let 
M, be the e. m. f. generated in the single-phase winding 
by one ampere of the a component of the polyphase 
current, provided that is the one which produces a 
magnetic field, when varying at the rate of one ampere 
per second. Thus defined the n-phase mutual induct- 
ance is 7 times the single-phase mutual inductance. 
If the single-phase winding is supplied with an alter- 
nating current having an angular velocity of w and the 
winding is rotated at this same velocity, cə, with respect 
to the polyphase winding, one of the oppositely rotating 
fields due to the single-phase current will be stationary 
with respect to the polyphase winding and thus will 
generate no e. m. f. in it. The other field will rotate 
at a relative angular velocity of 2 w and generate an 
e.m.f., E, in each phase of the polyphase winding. 
If the single-phase current is I, the single-phase mutual 


E 
inductance is M, = 9 ن‎ 2/9 ` 
/ 


If the single- and 
polyphase windings are fixed with respect to each other 
and one phase of the latter is directly opposite the 
former, the oppositely rotating fields produce equal 
e. m. fs. in time phase with each other in this phase 
of the polyphase winding. If E is still the e. m. f. 
generated in this phase of the polyphase winding the 
e. m. f. due to one of the component fields is E /2 and 
the single-phase mutual inductance is as before M, 


_ Ë . 
^ 2oI/2 ° 


angular velocity in the single-phase winding. If the 
two windings are fixed with respect to each other 
and balanced currents of 7 amperes per phase are sup- 
plied to the polyphase winding the polyphase mutual 
inductance is the e. m. f. generated in the single-phase 


where I and w are the current and its 


winding divided by 4ھ‎ ; where w is the angular 
velocity of the polyphase currents and 7 is the current 
per phase. Thus measured, the n-phase mutual in- 
ductance will be ہہ‎ times the single-phase inductance. 
Let us assume that we have a three-phase winding 
fixed in position and a single-phase winding which is 
rotating in the positive direction at an angular velocity 
of p radians per second. If the single-phase winding 
carries a current having an angular velocity of w radians 
per second there will be produced in the air-gap two 
rotating fields, one traveling in the positive direction 
at an angular velocity of (w + p) and the other in the 
negative direction at a velocity of (w — p) radians per 
second. If p is greater than w both of these fields 
travel in the positive direction. The following rela- 
tions can be worked out on the assumption that cə is 
either greater or less than p. In accordance with the 
established custom when analyzing polyphase induction 
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motors and generators, we will choose the first of these 
assumptions, viz., that w is greater than ۰ 

The alternating current in the single-phase winding 
having an angular velocity, w, produces two oppositely 
rotating fields each having this same angular velocity. 
The single-phase winding, having a positive angular 
velocity of » radians per second, gives to each of these 


fields an added angular velocity, so that one rotates. 


positively with an angular velocity of w + p and the 
other negatively with angular velocity of w — p radians 
per second. Each of these fields produces by its action 
on the polyphase winding a rotating field which has the 
same angular velocity. This is illustrated in Fig. 25. 
In this figure A, represents the fixed axis of one phase 
of the polyphase winding, and A; represents the axis 
of the single-phase winding; the latter having a posi- 
tive angular velocity of p radians per second. The 
oppositely rotating fields due to the single-phase cur- 
rent, which has a complex angular velocity of — به‎ 
+ J w radians per second, are represented by the vec- 


"2 ۱ À 
Pza. jp ı Oxy 
۲ - 06 - (o -p) ۱ 
m Le NS 
-0+j(w+p) @À 
- 02 - J (uo - p) 
-0+jko+p) 
4 j -O6+j(w-p) A -%+)(w+p) 
~-j(w-p) 3 
€ İc 
Tis la 
Fig. 25 


tors موقي‎ and $s. The former has an angular velocity 
of — a + j (w + p) and the latter a velocity of — a 
— j (o — p) radians per second. Since these com- 
ponent fields are proportional to the current in the 
single-phase winding, this current can be divided into 
similar components as indicated in the figure, Is. and 
Ia. Observe that each of these component currents 
is one-half of the maximum value of the single-phase 
current, and that at any moment the actual current 
in the single-phase winding is their vector sum. The 
first of these rotating fields, ,موف‎ produces a cur- 
rent in the short-circuited polyphase winding J,, where 


— Ia Mi (m + j p) 


la = Y + ہیں‎ + Lı 


— Ia M; (m + j p) 


Zi (m +j?) ( 39) 


M, is the single-phase mutual inductance, m = 
— ot + o) is the generalized angular velocity of the 
single-phase current, r, and L, are respectively the 


resistance and synchronous self inductance of the 
polyphase winding. The other rotating field دوق‎ 
produces a current in the short-circuited polyphase 
winding, ,ب7‎ which has a negative angular velocity of 
w — p radians per second. In this case the current is: 


— Iu Mi[- a — j (o— »)] 
rı -[-«—3(o—»lL 


Each of these polyphase currents produces two oppo- 
sitely rotating component currents in the single-phase 
winding. The sum of the components which rotate 
in the same direction must of course be equal to the 
initial components rotating in the corresponding direc- 
tion that were assumed to exist, viz., مور‎ and In. By 
the principle already established the rotating field due 
to İn, produces the same effect in the single-phase 
winding as would a current of the same magnitude 
which rotates positively and is always conjugate to 
I, with respect to the axis of the single-phase winding. 
It is readily seen that this conjugate current, ودر‎ is 


= I, M, [一 ot + 7 (o — 0?) [ 
rı + ]- a + (w— (م‎ lL, 
za I» Mi (m — 7 p) 


41 —— = ہے 
Zi (m-j») ( )‏ 


"nm 一 (40) 


Live = 


The stator current مر‎ and the conjugate current ود7‎ 
which replaces ,ورم‎ both rotate positively and their 
resultant effect on the single-phase winding is thus 
the same as the effect of their vector sum. The posi- 
tively rotating component of the single-phase current 
is thus given by: 


= (Tis + [u 3M, (— 0 +) Q) 
و۲‎ + ) a + 7 o) L. 
5 — (Lie + Tis.) 3/2 Mim 
Jaa لحك‎ 
5 22 


These eharaeteristic equations may be written thus: 


216 = 


(42) 


I. Z, (m4jp) 十 7 M, (m + j?) = O (43) 
Iu Z: (m-jp) + مر‎ Mi (m—jp) =O (44) 
" 5 5 ə M, 

Ix 7 (m) + iuc Ii) 2 m = O (45) 


In these equations it is assumed that any impedance 
outside the armature and field circuits is negligible, 
otherwise the zeros at the right of the equality signs 
would be replaced by the proper fall in terminal pres- 
sure. 

If مر‎ and 7. as determined by the first two equations 
are substituted in the third equation, the latter be- 


comes, on dividing by Tə," 
M. (m +7 Pp) 
Z, (mtJ p) 


Z» رم‎ — 8/2 My m( 


= 0 


M, (m — j 
P (m — j p) ) (46) 


Zi (m-Jp) 
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If the values of the impedances are substituted, this 
reduces to the following cubic in m: 


(ki (1 + 0) + kə) m? 


G 


m? + 


7 (Ki? + 2 6, ks + o p?) m 7 (k? + 9?) ks 


-0 
p c 
(47) 
3 M ° 
where k, = ri/Li, ke = r;/L; and ç = 1— 2L, Le 


r, and r, are the resistances of one phase of the poly- 
phase winding and of the single-phase winding; L, 
and L, are the synchronous self-inductance of the 
polyphase winding, and the self inductance of the 
single-phase winding. At least one root of a cubic 


equation must be real and since the coefficients are . 


all positive this root must be negative. This indicates 
that one of the components of the single-phase current 
is a diminishing direct current. The other roots of the 
cubic equation are conjugate imaginaries. This indi- 
cates two alternating currents, having positive and 
negative angular velocities but the same rate of diminu- 
tion. In the polyphase winding, however, the three 
components of the current will all have different angu- 
lar velocities. The first component of the single-phase 
current will produce a current having an angular veloc- 
ity of p radians per second and the same rate of diminu- 
— tion that it itself has. The second component of the 
single-phase current, having a positive angular velocity 
will produce a current which has an angular velocity 
of approximately 2 p radians per second. The third 
component of the single-phase current having a nega- 
tive angular velocity, produces a current in the poly- 
phase winding of very small angular velocity. The 
sum of the angular velocities of these two latter com- 
ponents of the polyphase current is exactly 2 ,م‎ and 
each diminishes at the same rate as does the alternating 
current in the single-phase winding. 

A literal cubic equation is practically impossible 
to solve but if numerical values are given to kı, k, and 
c, the real root ean be found to any degree of approxi- 
mation by Horner’s method. The equation may then 
be reduced to a quadratic, the roots of which can be 
readily found. Approximate values of the roots may 
be found as follows: The largest term in any of the 
coefficients is p?. With a 60-cycle machine this would 
be 142,122. The ratios of resistance to self-inductance 
kı and ke, are of the order ten, while o is of the order 
one-tenth. For these values of kı, k, and c the real 
root of the characteristic equation is approximately 

m = 一 k, / G (48) 
The degree of approximation of eourse depends upon 
the values of the constants. If kı = k: = 1.0, © 
= 0.1 and p = 877, 


m = — 10.06 
whereas the approximate value would be 
m = 10.0 


In this case the error is about one-half of one per cent; 
which is better than the analysis warrants. 

Since the sum of the roots must be the coefficient of 
m? in the characteristic equation. with its sign changed 
and since the remaining roots are conjugate imaginaries 
of the form — as + j w we have 


k, (1 + k. 
- لظ ے ہی 2 - م/م‎ ٢ (49) 
That is: 
kü + c 
uru (50) 


Purthermore since the product of the roots is equal to 
the last term in the characteristic equation we have: 
ريه‎ (as? + o) = (k + p?) 7 
From this, the approximate value of w is seen to be 
o = V p — (as — kit) 
ay — ky? , 

= p— ^ 2p approximately (51) 
For the assumed values of kı, kə, c and p the calculated 
values of o» are w are: f 

Qo = 5 š 497 

and w = 376.73 
For the same values of the constants the approximate 
values of o, and w are: 


|. 1x1.12 e 
a1 mo 7 
"m 5 5۶ - 1 
@ =o 2x81 
= 376.96 


It is readily appreciated that for such values of the 
constants as were chosen the roots are approximately 


m" =—k/o (52) 
" kı (1 + a) ' 
“x 2g +) 2 (53) 


It may be desirable, however, to extend the analysis 
so that it will apply to the case in which it is necessary 
to consider the constants of circuits outside of the 
machine itself. In that event the approximate values 
of the roots may not be sufficiently accurate. Again 
if the constants are considerably larger, as for example: 
Let k, = 15, k, = 20, o = 0.15 and p = 371 the 
actual values of m are: 
m" = — 143.2 + 0 
m" = — 52.57 + 360.3 
The approximate values of m are: 
m" = — 133.3 
m = — 57.5 2238 © 
In this ease the error in the approximate caleulation is 
quite considerable. 
The general principle that determines the magnitude 
of the various component currents is that the resulting 
distribution of magnetic field in all parts of the machine 
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cannot abruptly change, and what is even more im- 
portant the component fields due to the different cur- 
rents eannot abruptly change. That is, in mathemat- 
ical language, the magnetic field at every point is a 
continuous function of the time, and since the magnetic 
fields are assumed to be proportional to the currents, 
the latter are also continuous functions of the time. 
If it were possible to construct mutually inductive 
circuits which had a leakage coefficient of zero, the 
currents in these circuits might abruptly change, 
1. e., they might be discontinuous functions. 

In order to illustrate the method of determining 
the component currents we will calculate their approxi- 
mate values when the alternator is suddenly short- 
circuited. The general method of treating the single- 
phase current is to represent it by pairs of oppositely 
rotating vectors. In the present problem there are two 
such pairs one having no velocity and the other an 
angular velocity of p radians per second. That is, the 
transient current /,: in the single-phase winding may 
be equated to four components: 

L, = Ine’ + Toy’ + L,” + "رو‎ (54) 
İs." and I,’ are conjugate and have zero velocity with 
respect to the single-phase axis. "مو‎ and Iz,” are con- 
Jugate and have an angular velocity of approximately 
p radians per second with respect to the same axis. 
Each of these transient components in the single-phase 
winding is the cause of a corresponding component 
in the polyphase winding. Thus, similarly: 

I 一 Tha’ T Io "T I" T 1," (55) 
I,' and İn” both have the same angular velocity of 
p radians per second, "مر‎ has an angular velocity of 
approximately 2 and I," an angular velocity of 
approximately zero radian per second. 

The fundamental relation at the moment of short 
cireuit is that: 


lio ləbi. (56) 
Izo = Ter + I». (57) 


where I, is the current in the polyphase winding, before 
` short circuit and Tı, is the steady current after the 
transient I,: occasioned by the short circuit has dis- 
appeared. I, and I», are similar values for the single- 
phase winding. 

If the field excitation is not changed at or after short 
circuit مم‎ = Iz, and thus the transient current in the 
single-phase winding must be zero at the moment of 
short circuit. Since the components of the single- 
phase current are conjugate in pairs it follows that 
at the moment of short circuit: 


Iza’ == EE Ij" (58) 
"وو]‎ = — Ia" (59) 


The corresponding pairs of single- and polyphase 
component currents must satisfy the two equations 
(43) and (44), the a eomponents satisfying (43) and the 


and 


and 


b components, (44). Thus: 
Tie’ (rı + (— ka/ o tJ») Li + Iza’ M, (— 0 


十 72) = O (60) 
İp" [rı + (一 k,/o + ) p) L4] + I, M) (— k,/ o 
+jp) = O (61) 
” io) q 
I. E ےا‎ kı 20 t j2»)1L | 
(1 + o) : 
+ Tne! Mı (= kı م96‎ +j2p) = 0 (62) 
1 + 
T (rı + (— k, 5. L.) 
(1 + o) 
+ Ia" M(- kı و‎ — ) = O (63) 


If each of these equations is divided by L,, it will be 
seen that since k) is much smaller than k,/o, k./o or p, 
the following relations between the polyphase and 
single-phase eomponents are approximately true: 

I! bə M/L, Toa! 
I, EE M,/Ly In' 
I," =— gE M,/L, Toa" 


X (11) - 
İn" = — ,ر‎ Ia” 


If the conditions set forth in equations (55), (58) 


but 


and (59) are applied to these four equations we have: 


Tə = Mit, I Lj (64) 

If it is deemed necessary à more accurate value of 
I, can be calculated by substituting in equations 
(60) to (63) inclusive the actual values of the generalized 
angular velocity as determined by equation (47) and 
by recognizing that I,,' and I,,' are conjugate with 
respect to the single-phase axis. 

Since the current in the polyphase winding before 
short circuit is much smaller than the steady short- 
circuit current: 

I, approximately, at the moment of short‏ — = ہآ[ 
circuit. The steady short-circuit current is calculated‏ 
by dividing the e. m. f. generated in the polyphase‏ 
winding on open circuit by the synchronous impedance.‏ 
The open circuit e. m. f. is‏ 

E = - 7 pM, وآ‎ ۱ )65( 
Since the resistance is much less than the synehronous 
reactance the steady short-circuit current is approxi- 


mately. 
7 p Mı Iz 
I, = 一 E L (66) 
nəl Mı/ Lı 120 
From this it follows that: 
1-— 7 
Ix’ = In 5. (67) 
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The direct-current component of the single-phase 
current has an initial value of Iz’ + Iz’ which is 


(1 — o) 
27 


120 The alternating component of the 
single-phase current is equal and opposite to this 
initially, but one-half cycle later the vectors repre- 
senting all of the component currents are in phase. 
During this time, however, the direct and alternating 
components have shrunk different amounts. The first 
maximum value of the single-phase current will be, 
approximately: 

1— > Jem Ja(l +o) x 
Izo تح‎ €? +e 2 ? (68) 
For ki = ke = 1۳ g = 0.1, this 15 17 Leo. 

The component of the polyphase current which has 
an angular velocity of p radians per second is İT.” 
l— o 
" İr. 


+ İn" which initially equals The com- 


ponent which has an angular velocity of approximately 


_ Ic 
2 pis مر‎ which initially equals 一 EP 


component which has a very small angular velocity is 


1 + c 
20 


of these components shrinks at a rate of k;/o while each 


kı (1 + o) 
20 ° 


tice that initially the first component is opposite in 
phase to each of the others. Thus one-half cycle after 
the transient begins, the three components will be in 
phase. If the transient begins at such a time that the 
steady short-circuit current in one phase would natur- 
ally be passing through its maximum value, this phase 
will carry the maximum possible transient current 
one-half cycle later. This maximum transient current 
` is approximately: 


] 一 € E: ] (1 +o) r 
m [ € 5 b + — c" 2o ? ] 
g 0 


for k, = k, = 1, = 0.1 this is 18 İ,, 7. e. eighteen 
times the maximum value of the steady short circuit 
current. 


Ti The 


Tə" which initially equals — I,. The first 


of the others shrinks at a rate of No- 


(69) 


Two SINGLE-PHASE WINDINGS HAVING A CON- 
STANT RELATIVE ANGULAR VELOCITY 


A current of angular velocity w in one single-phase 
winding produces in the other winding currents having 
angular velocities of w + p and o — p; where p is 
the relative angular velocity of the two windings. Each 
of these currents reacts on the first winding and pro- 
duces additional currents having angular velocities of 


w + 2p and w— 2p. These currents in turn produce 
others in the second winding having angular velocities 
Of w + 3p and w— 3p, and so on ad infinitum. It ` 
is then evident that the solution of this problem cannot 
be reached by the method here developed unless all 
but a few of the component currents are neglected. 
The solution will not be attempted at this time. 

Carson® has suggested a functional form of solution 
of this problem and Dreyfus‘ also gives a solution based 
on certain assumptions. 

In conclusion the writer hardly needs to say that it 
has been his sole purpose to show how certain low- 
frequency transient conditions in electric machinery 
may be analyzed by a vector method. Before accurate 
results can be obtained by this method considerable 
experimental work is probably necessary in order to 
determine the most satisfactory methods of measuring 
the required constants. An interesting development 
would be the solution of the transient conditions in 
certain unsymmetrical arrangements of polyphase 
machines similar to those so ably discussed by Mr. 
Fortescue.’ 


Transient conditions may be the result of a variety 
of causes, one of which is short circuit. The vector 
method, however, is equally applicable in any case in 
which the original premises are reasonably true. 
The actual working out of a numerical solution is 
rather complex and this has been purposely avoided in 
order that the general principles might not be obscured 
a mass of calculation. Furthermore, it seemed best 
at this time not to attempt to generalize the theory 
and present it in such form that it would apply to 
combinations of electric machines that might even 
be unsymmetrical. 


Appendix 

For those who are interested in comparing the vector 
method of solution with the differential equations which 
apply to these problems the following development is 
given. 

Consider a symmetrically wound three-phase induc- 
tion motor having Y-connected stator and rotor phases 
as shown diagrammatically in Fig. 26. The stator and 
rotor currents are represented by the letters x and y 
respectively and their assumed positive directions are 
indicated by the arrows. Fig. 27 shows a develop- 
ment of the air gap. The progression of the phases, 
1, 2 and 8 is toward the right for both stator and rotor, 
and the latter moves at a constant angular velocity 
of » radians per second toward the right. 


Let rı, rs, L'; and L;' represent the resistances and 
self-inductances of one phase of stator and rotor 
respectively. 

Let M, represent the mutual inductance between any 

5. Bibliography No. 7. 


6. Bibliography No. 4. 
2. A.I. E. E. Vol. XXXVII, p. 1027. 
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two phases of the stator. Due to symmetry these 
three mutual inductances are equal. 

Let M, represent the mutual inductances between 
any two phases of the rotor. These mutual inductances 
are likewise equal. 

Let 2/3 M represent the maximum value of the 
mutual inductance between one phase of the stator and 
one phase of the rotor. At a time when the angular 


Fic. 26 


displacement of these phases of stator and rotor is 
pt the momentary value of this mutual inductance 
is 2/3 M cos pt, if the space distribution of the air- 
gap flux density is sinusoidal. This assumption of a 
sinusoidal distribution of flux density is not much in 
error if the motor is well designed. It is also necessary 
to assume that rı, re, Lı’, Le", M., M, and 2/8 M are 
constant. This is commonly done in engineering prac- 
tise. 


Stator 
— مسا‎ 
پر 120 120 ہے‎ 


‘Y2 ۸ 
b pt el” 120°- 2 1207 ¬+ 
Rotor 


Fic. 27 


With the stator phases short-circuited as indi- 
cated in Fig. 26 the following differential equations 
may be written,—considering first phases 1 and 2 and 
then phases 2 and ۰ 


, dx, d Xə _ 0 2: | 
ina + Lı it -M.-x M I 
d, 2 2 ۱ 
+ [yS M cos p t + y» — M cos (pt + 120°) 
"m 3 
+ ys M cos (pt + 240 ([ } 
, d x» d x: dx, 
= (rizə + Lı dt — M. dt — M, dt 
d , 2 RE ۱ 
+ رس‎ (Us چ-‎ M cos pt + ya M cos (pt + 120°) 


+ yı £ M cos (p t + 240?) } = (1) 


There is a negative sign before M, for the reason that 
positive currents in phases 2 and 3, for example, 
produce a flux through phase 1 which is opposite in 
direction to the flux that phase 1 produces through itself. 

The equation which applies to phases 2 and 3 is: 

d Z2 d v3 dx 1 

dt a ون‎ 


d 2 2 ° 
+ -gr Wg M cos pt + yi M cos (p t + 1209) 


I Tı Lo + Ly’ 


yis M eos (pt + 240) [۱ 
= dz _ y, 4: 
t 


d. wu 


— Í rıza + Ly’ — M. 


+ tus j M cos p t +y- M cos (p t + 120°) 


+ yes M cos (pt + 2409) [ ( = 0 (2) 


Since the sums of the stator and rotor eurrents are each 
equal to zero, we may substitute: 

ts = —21— Tə and y; = — و‎ — 
Also replace L,' + M. by رط‎ Multiply the second 
equation by 2 and add it to the first, again, subtract 
the second equation from the first. The results of 
these operations are 77 


mit + مرچ ول‎ + Vr M 一 一 2 [yı cos (pt + 30 ) 
2777770775 (3) 

d x» | 2 نے‎ 
nız: + Li- t Fg Non M 一 一 -9 [y1 cos (p t — 90) 
+ y2 cos (p t — 30°)] =0 (4) 


By a similar procedure differential equations may be 
written for phases 1 and 2 of the rotor. These equa- 
tions will be found to be similar to (3) and (4) except 
that the signs before the phase angles will be reversed. 
Thus we will find " 


۲۵ 7 + Le ous b M —_ |z eos (p t — 30°) 
7777757 (5) 
2 

Ta yə + Le اماد ولا کے‎ e Vs M -Ir [xı eos (p t + 90°) 


+ Za cos a + 80°] = 0 (6) 

The following treatment was suggested to the author 

by Dr. F. L. Hitchcock as being the simplest method 

of reducing these equations to one containing a ə 
variable. 


Let سے نگ رر بت رم‎ 2, zi 
(Z; is an impedance operator, 

Similarly rı وھ‎ + Ly ون کے‎ = Z: 2; 

Also T2 V1 + Le al = الام‎ 

And T» a + Le و کے‎ = Z: ولا‎ 


If D and D- represent the direct and inverse dif- 
ferential operators, we shall have on integrating equa- 
tions (3) and (4): 

Du Z, z, + 2// 3 M [yi cos (p t + 30°) 


+ ون‎ cos (pt + 90°] = 0 (7) 
D^ Z, zx + 2/۰۷ 3 M [yı cos (p t — 90?) 
+ ye cos (pt — 30°] = 0 (8) 


Solve these two equations for z: and z+ 
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— — 90°) DZ 
yı = VM —— [eos (pt + 90°) 1X2 
— eos (pt — 30 ) D-1 Z, >l (9) 


—icos (pt — 90 ( D! Z, >, 


چو کہ 
M‏ 8 يه y» 一‏ 
eos (p t+ 80?) D3 Z, xə)‏ — 

Differentiation of equation (9) gives: 


(10) 


dv: )سح كك‎ 
dt ليه‎ 3M 
— p sin (p t + 90°) D! Zi x: — cos (p t — 30°) Zi 1 
T p sin (p t — 30°) D-1 Z, xı) 
The expression for y» may be differentiated in a like 
manner. 


cos (p t + 90°) Z, zə 


d yi d y» 
If the values of yı, Y2, prs and “dı are now sub- 


stituted in equations (5) and (6), the resulting equations 
may be written as follows: 

In order to simplify the writing of these equations 
let: 


cos (p t — 30”) 

cos (p t + 30”) 

cos (p t — 90”) 

cos (p t + 90°) = d 
Equation (5) becomes: 
) ) - r:a + pLoa) D? Zi — L; a Z, 

| + M? (a D — P 01([ 21 
+Í (r:d - p L, di) DUZ) + Lid 


sin (pt — 30°) = a, 
sin (pt + 30°) = b, 
sin (pt — 90°) = c, 
sin (p t + 90°) = d, 


۱ 
o çən 


+ M? (c D — pe))) zs = 0 (11) 
Equation (6) becomes: | 
{ (roe — L p c) D! Z, + Doc ZA 
(— rs b + L, p b) D? Z, — Lab Z, 
+ M? (dD— pd,)} zi + 
-TM*(bD—pb))x.—-0 ` (12) 


In the process of eliminating z, and x: from equations 
(11) and (12) all of the harmonic terms in the coeffi- 
cients vanish, giving as a final result: 

(r? + Lè p) یں‎ Zi? + Lè Z + (D? + p°) ۰ 

+ 2 ra L, Z, DA 21— 2 (ro D + p? Lə) M? D~ Zi 
— 2⁄2 L, Z, D = 0 (13) 
Differentiate this equation twice, collect the terms 
in descending powers of D and after dividing each 

M? 
L, L. 
.وم = ویو‎ This gives: 

e? Dt + 2 o (ki + kə) Dš + [ ) + k) + 2o ki k, 

+ P o?) D? + [2 ki ke (kı + kə) + 2 ثم‎ k, o] D 
۱ . ky? (ko? + p?) = 0 (14) 
If equation (26) in the body of the paper is written 
for the conjugate roots, z. e., for a negative value of the 


term by L L4, let 1 — =o, T/L = k, and 
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rotor angular velocity it would be: 

o m? + (kı + ks + j po) m + ki ks + j pk = 0 5 
The product of equations (26) and (26;) will contain 
all four roots, and will be found to be identical with the 
differential equation (14). 

This proves that the two methods are identical i in 


‘regard to the results they produce, and that, of the two, 


the vector method is much less laborious. Further- 


more, it is much simpler to determine the constants of 


integration by the vector method. | 


SYNCHRONOUS ALTERNATOR | 


In this case there is only one phase on the rotor, 
viz. the field winding. Making the same assumptions 
in regard to the resistances and inductances as were 
made before, the following differential equations may 
be written: 


dx, 


۲۱ 2۱ + Lı dı + M -ir (y cos p t) = 0 
d xə 

rı و22‎ + L. qi FM Tr ; ly cos (pt — 120) ] = 0 
d y. o 

rey + Le t +M -ir [zi cos (p t — 80°) 


+ 22 cos (p t — 9001 = 0 
Bu these equations are treated by the method already 
outlined the resulting differential equation is: 
c D3 + (kı + kə + o ka) D? + (k + 2k, ke 
تر رہ ےی‎ + p*) = 0 


2 
LiL, 


Again the equation obtained by the vector -— and 
the differential equation are identical. 


In this case g = 1 3/2 —— 
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